Taller

Taller 1: Sistemas de Ecuaciones
Brayan Escola
Resolver los ejercicios planteados:
1. Una caja de 13 Mg contiene una pieza de equipo y está sostenida por tres grúas cuyos cables se unen en
el anillo D, como se muestra en la figura. Determine las tensiones en los cables DA, DB y DC.
Método de Gauss sin pivoteo

clc %permite borrar el area de trabajo
1
clear %permite borrar las variables almacenadas
format long %permite utilizar la maxima capacidad de la maquina
A=[3/5 -2/7 -4/13;4/5 6/7 12/13;0 -3/7 3/13]
A = 3×3
0.600000000000000 -0.285714285714286 -0.307692307692308
0.800000000000000 0.857142857142857 0.923076923076923
0 -0.428571428571429 0.230769230769231
b=[0;127400;0]
b = 3×1
0
127400
0
gaussnpiv%script normal
La Matriz C, que corresponte a la matriz aumentada [Ab] es =

1.0e+05 *
0.000006000000000 -0.000002857142857 -0.000003076923077 0

0.000008000000000 0.000008571428571 0.000009230769231 1.274000000000000
0 -0.000004285714286 0.000002307692308 0
ETAPA 1=
La matriz correspondiente a esta etapa antes del proceso:

1.0e+05 *
0.000006000000000 -0.000002857142857 -0.000003076923077 0

0.000008000000000 0.000008571428571 0.000009230769231 1.274000000000000
0 -0.000004285714286 0.000002307692308 0
Los Multiplicadores correspondientes a esta etapa son:
m(2,1)= 1.333333333333333
m(3,1)= 0
La matriz correspondiente a esta etapa despues del proceso:

1.0e+05 *
0.000006000000000 -0.000002857142857 -0.000003076923077 0

0 0.000012380952381 0.000013333333333 1.274000000000000
0 -0.000004285714286 0.000002307692308 0
ETAPA 2=

1.0e+05 *
0.000006000000000 -0.000002857142857 -0.000003076923077 0

0 0.000012380952381 0.000013333333333 1.274000000000000
0 -0.000004285714286 0.000002307692308 0
2
m(3,2)= -0.346153846153846

1.0e+05 *
0.000006000000000 -0.000002857142857 -0.000003076923077 0

0 0.000012380952381 0.000013333333333 1.274000000000000
0 0.000000000000000 0.000006923076923 0.441000000000000
SOLUCI?N:
La matriz Ab final:
1.0e+05 *
0.000006000000000 -0.000002857142857 -0.000003076923077 0

0 0.000012380952381 0.000013333333333 1.274000000000000
0 0.000000000000000 0.000006923076923 0.441000000000000
La solucion de X1 hasta Xn es:
X1= 49000
X2= 3.429999999999999e+04
X3= 63700
Método de Gauss con pivoteo

A=[3/5 -2/7 -4/13;4/5 6/7 12/13;0 -3/7 3/13]
A = 3×3
0.600000000000000 -0.285714285714286 -0.307692307692308
0.800000000000000 0.857142857142857 0.923076923076923
0 -0.428571428571429 0.230769230769231
b=[0;127400;0]
b = 3×1
0
127400
0
sol = gausspivoteo(A,b) %script con funcion
Ab = 3×4
105 ×
0.000006000000000 -0.000002857142857 -0.000003076923077 0
3
0.000008000000000 0.000008571428571 0.000009230769231 1.274000000000000
0 -0.000004285714286 0.000002307692308 0
1.0e+05 *
0.000008000000000 0.000008571428571 0.000009230769231 1.274000000000000

0 -0.000009285714286 -0.000010000000000 -0.955500000000000
0 0 0.000006923076923 0.441000000000000
sol = 3×3
104 ×
4.900000000000000 0 0
0 3.430000000000000 0
0 0 6.369999999999999
% y=diag(sol)
n=3;
for i=1:n
Xi=X(1,i);
fprintf('\nX%g=',i)
disp(Xi);
end
X1=
49000
X2=
3.429999999999999e+04
X3=
63700
Método de Gauss Seidel

A=[3/5 -2/7 -4/13;4/5 6/7 12/13;0 -3/7 3/13]
A = 3×3
0.600000000000000 -0.285714285714286 -0.307692307692308
0.800000000000000 0.857142857142857 0.923076923076923
0 -0.428571428571429 0.230769230769231
b=[0;127400;0]
b = 3×1
0
127400
0
n=3; %número de ecuaciones

GaussSeidel
K X(1) X(2) X(3) Norma

1 0.0000148633.3333276033.3333313506.4096
2 212333.3333-346811.1111-644077.77781066374.7459
3 -495444.44441304670.37042422959.25933554582.4864
4 1863814.8148-4200267.9012-7800497.530911848608.2882
5 -6000382.716014149526.337426277691.769539495360.9605
6 20213609.0535-47016454.4582-87316272.5652131651203.2017
7 -67166363.5117156870148.1939291330275.2172438837344.0057
8 224100211.7055-522751860.6462-970824884.05731462791146.6856
4
9 -746788372.35181742654835.48753236358980.19104875970488.9521
10 2489506907.8393-5808700818.2916-10787587233.970116253234963.1737
11 -8298144026.130819362484694.305335958900146.567054177449877.2457
12 27660692420.4361-64541467014.3510-119862724455.2233180591499590.8189
13 -92202095734.7872215138372014.5034399542690884.0777601971665302.7297
14 307340531449.2906-717127758081.6781-1331808693580.25902006572217675.7659
15 -1024468225830.96852390426008905.59384439362587967.53126688574058919.2207
16 3414894298436.5625-7968086547718.6465-14797875017191.769522295246863064.0703
17 -11382980782455.209026560288641028.824249326250333339.234474317489543546.8906
18 37943269487184.0312-88534295321462.7500-164420834168430.7812247724965145156.3125
No se alcanzó la convergencia
Ejercicio 2

A=[-5 5 3; 5 6 1; 3 1 7]
A = 3×3
-5 5 3
5 6 1
3 1 7
b=[1;2;3]
b = 3×1
1
2
3

-5 5 3 1
5 6 1 2
3 1 7 3
ETAPA 1=

-5 5 3 1
5 6 1 2
3 1 7 3
5
m(2,1)= -1
m(3,1)= -0.600000000000000

-5.000000000000000 5.000000000000000 3.000000000000000 1.000000000000000
0 11.000000000000000 4.000000000000000 3.000000000000000
0 4.000000000000000 8.800000000000001 3.600000000000000
ETAPA 2=

-5.000000000000000 5.000000000000000 3.000000000000000 1.000000000000000
0 11.000000000000000 4.000000000000000 3.000000000000000
0 4.000000000000000 8.800000000000001 3.600000000000000
m(3,2)= 0.363636363636364

-5.000000000000000 5.000000000000000 3.000000000000000 1.000000000000000
0 11.000000000000000 4.000000000000000 3.000000000000000
0 0 7.345454545454546 2.509090909090909
SOLUCI?N:
La matriz Ab final:
-5.000000000000000 5.000000000000000 3.000000000000000 1.000000000000000
0 11.000000000000000 4.000000000000000 3.000000000000000
0 0 7.345454545454546 2.509090909090909
X1= 0.153465346534653
X2= 0.148514851485149
X3= 0.341584158415842

A=[-5 5 3; 5 6 1; 3 1 7]
A = 3×3
-5 5 3
6
5 6 1
3 1 7
b=[1;2;3]
b = 3×1
1
2
3
Ab = 3×4
-5 5 3 1
5 6 1 2
3 1 7 3
-5.000000000000000 5.000000000000000 3.000000000000000 1.000000000000000
0 11.000000000000000 4.000000000000000 3.000000000000000
0 0 7.345454545454546 2.509090909090909
sol = 3×3
0.153465346534653 0 0
0 0.148514851485149 0
0 0 0.341584158415842
% y=diag(sol)
n=3;
for i=1:n
Xi=X(1,i);
fprintf('\nX%g=',i)
disp(Xi);
end
X1=
0.153465346534653
X2=
0.148514851485149
X3=
0.341584158415842

A=[-5 5 3; 5 6 1; 3 1 7]
A = 3×3
-5 5 3
5 6 1
3 1 7
b=[1;2;3]
b = 3×1
1
2
3
7
GaussSeidel

1 -0.2000 0.5000 0.4429 0.6972
2 0.5657 -0.2119 0.2164 1.0698
3 -0.2821 0.5323 0.4734 1.1570
4 0.6164 -0.2592 0.2014 1.2279
5 -0.3383 0.5817 0.4905 1.3047
6 0.6760 -0.3117 0.1834 1.3862
7 -0.4017 0.6375 0.5097 1.4728
8 0.7433 -0.3710 0.1630 1.5648
9 -0.4732 0.7005 0.5313 1.6625
10 0.8193 -0.4380 0.1400 1.7664
11 -0.5540 0.7716 0.5558 1.8767
12 0.9051 -0.5135 0.1140 1.9939
13 -0.6451 0.8519 0.5833 2.1185
14 1.0019 -0.5988 0.0847 2.2508
15 -0.7480 0.9426 0.6145 2.3914
16 1.1113 -0.6951 0.0516 2.5408
17 -0.8642 1.0449 0.6497 2.6995
18 1.2347 -0.8038 0.0143 2.8682
Ejercicio 3

A=[3.33 4.81 -2.22; 4.81 10 745; -2.22 7.45 15]
A = 3×3
102 ×
0.033300000000000 0.048100000000000 -0.022200000000000
0.048100000000000 0.100000000000000 7.450000000000000
-0.022200000000000 0.074500000000000 0.150000000000000
b=[5;0;2]
b = 3×1
5
0
2
8
1.0e+02 *
0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000

0.048100000000000 0.100000000000000 7.450000000000000 0
-0.022200000000000 0.074500000000000 0.150000000000000 0.020000000000000
ETAPA 1=

1.0e+02 *
0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000

0.048100000000000 0.100000000000000 7.450000000000000 0
-0.022200000000000 0.074500000000000 0.150000000000000 0.020000000000000
m(2,1)= 1.444444444444444
m(3,1)= -0.666666666666667

1.0e+02 *
0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000

0 0.030522222222222 7.482066666666667 -0.072222222222222
0 0.106566666666667 0.135200000000000 0.053333333333333
ETAPA 2=

1.0e+02 *
0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000

0 0.030522222222222 7.482066666666667 -0.072222222222222
0 0.106566666666667 0.135200000000000 0.053333333333333
m(3,2)= 3.491445212959590

1.0e+03 *
0.003330000000000 0.004810000000000 -0.002220000000000 0.005000000000000

0 0.003052222222222 0.748206666666667 -0.007222222222222
0 0 -2.598802584637785 0.030549326538041
SOLUCI?N:
9
La matriz Ab final:
1.0e+03 *
0.003330000000000 0.004810000000000 -0.002220000000000 0.005000000000000

0 0.003052222222222 0.748206666666667 -0.007222222222222
0 0 -2.598802584637785 0.030549326538041
X1= 0.749222866050857
X2= 0.515382830044366
X3= -0.011755154746508

A=[3.33 4.81 -2.22; 4.81 10 745; -2.22 7.45 15]
A = 3×3
102 ×
0.033300000000000 0.048100000000000 -0.022200000000000
0.048100000000000 0.100000000000000 7.450000000000000
-0.022200000000000 0.074500000000000 0.150000000000000
b=[5;0;2]
b = 3×1
5
0
2
Ab = 3×4
102 ×
0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000
0.048100000000000 0.100000000000000 7.450000000000000 0
-0.022200000000000 0.074500000000000 0.150000000000000 0.020000000000000
1.0e+02 *
0.048100000000000 0.100000000000000 7.450000000000000 0

0 0.120653846153846 3.588461538461539 0.020000000000000
0 0 -4.551425374561683 0.053502709595155
sol = 3×3
0.749222866050857 0 0
0 0.515382830044365 0
0 0 -0.011755154746508
% y=diag(sol)
n=3;
for i=1:n
Xi=X(1,i);
fprintf('\nX%g=',i)
10
disp(Xi);
end
X1=
0.749222866050857
X2=
0.515382830044366
X3=
-0.011755154746508

A=[3.33 4.81 -2.22; 4.81 10 745; -2.22 7.45 15]
A = 3×3
102 ×
0.033300000000000 0.048100000000000 -0.022200000000000
0.048100000000000 0.100000000000000 7.450000000000000
-0.022200000000000 0.074500000000000 0.150000000000000
b=[5;0;2]
b = 3×1
5
0
2

GaussSeidel

1 1.5015 -0.7222 0.7143 1.8128
2 3.0209 -54.6654 27.7309 60.3496
3 98.9498-2113.5459 1064.5057 2307.1830
4 3764.0716-81116.193840844.925588528.7145
5 144399.2873-3112403.00801567198.05513396797.2989
6 5540493.4387-119421232.452460132538.6136130333268.0230
7 212585696.3419-4582127846.65562307252847.03095000816705.3991
8 8156797678.0245-175813756786.928588528038593.9888191878620873.6452
9 312971896645.2797-6745878357538.54493396772758281.11237362278467155.7793
10 12008561688633.4746-258835688664175.5625130332325833125.0781282486626082431.0625
11 460761989737005.1250-9931384791631316.00005000780554324630.000010838858414746068.0000
12 17679187290794988.0000-381061840384057280.0000191877233776452768.0000415881110423366080.0000
13 678340814183495680.0000-14621135847967991808.00007362225244989927424.000015957132327853557760.0000
14 26027568610391494656.0000-561005041253347885056.0000282484583976834072576.0000612266500561702092800.0000
15 998663671128280793088.0000-21525458732086842097664.000010838780060263449100288.000023492333083916815564800.0000
16 38318182653189957877760.0000-825920160345811354386432.0000415878104004425108422656.0000901388060949571715465216
17 1470247856502455272472576.0000-31690107967307356804481024.000015957016973191684795924480.0000345857703243054445
18 56412611712682855554351104.0000-1215932230736581161755082752.0000612262074465979053176258560.000013270372226424
Ejercicio 4
11
A=[72 0 0 9 0 0;0 2.88 0 0 0 -4.5; 0 0 18 9 0 0; 0 0 9 12 0 0;0 0 0 0 33 0;0 -4.5 0 0 0 33]
A = 6×6
72.000000000000000 0 0 9.000000000000000
0 2.880000000000000 0 0
0 0 18.000000000000000 9.000000000000000
0 0 9.000000000000000 12.000000000000000
0 0 0 0
0 -4.500000000000000 0 0
% A=[72 0 0 9 0 0;0 2.88 0 0 0 -4.5;...

% 0 0 18 9 0 0; 0 0 9 12 0 0;0 0 0 0 33 0;0 -4.5 0 0 0 33]
b=[2;0.5;1;0;1.2;5]
b = 6×1
2.000000000000000
0.500000000000000
1.000000000000000
0
1.200000000000000
5.000000000000000

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 9.000000000000000 12.000000000000000 0
0 0 0 0 33.000000000000000
0 -4.500000000000000 0 0 0 33.0000000000
ETAPA 1=

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 9.000000000000000 12.000000000000000 0
0 0 0 0 33.000000000000000
0 -4.500000000000000 0 0 0 33.0000000000
12
m(2,1)= 0
m(3,1)= 0
m(4,1)= 0
m(5,1)= 0
m(6,1)= 0

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 9.000000000000000 12.000000000000000 0
0 0 0 0 33.000000000000000
0 -4.500000000000000 0 0 0 33.0000000000
ETAPA 2=

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 9.000000000000000 12.000000000000000 0
0 0 0 0 33.000000000000000
0 -4.500000000000000 0 0 0 33.0000000000
m(3,2)= 0
m(4,2)= 0
m(5,2)= 0
m(6,2)= -1.562500000000000

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 9.000000000000000 12.000000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 25.9687500000
ETAPA 3=
13
72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 9.000000000000000 12.000000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 25.9687500000
m(4,3)= 0.500000000000000
m(5,3)= 0
m(6,3)= 0

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 0 7.500000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 25.9687500000
ETAPA 4=

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 0 7.500000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 25.9687500000
m(5,4)= 0
m(6,4)= 0

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 0 7.500000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 25.9687500000
ETAPA 5=

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
14
0 0 0 7.500000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 25.9687500000
m(6,5)= 0

72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 0 7.500000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 25.9687500000
SOLUCI?N:
La matriz Ab final:
72.000000000000000 0 0 9.000000000000000 0
0 2.880000000000000 0 0 0 -4.5000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 0 7.500000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 25.9687500000
X1= 0.036111111111111
X2= 0.521460088247092
X3= 0.088888888888889
X4= -0.066666666666667
X5= 0.036363636363636
X6= 0.222623345367028

A=[72 0 0 9 0 0;0 2.88 0 0 0 -4.5;...
0 0 18 9 0 0; 0 0 9 12 0 0;0 0 0 0 33 0;0 -4.5 0 0 0 33]
A = 6×6
72.000000000000000 0 0 9.000000000000000
0 2.880000000000000 0 0
0 0 18.000000000000000 9.000000000000000
15
0 0 9.000000000000000 12.000000000000000
0 0 0 0
0 -4.500000000000000 0 0
b=[2;0.5;1;0;1.2;5]
b = 6×1
2.000000000000000
0.500000000000000
1.000000000000000
0
1.200000000000000
5.000000000000000
Ab = 6×7
72.000000000000000 0 0 9.000000000000000
0 2.880000000000000 0 0
0 0 18.000000000000000 9.000000000000000
0 0 9.000000000000000 12.000000000000000
0 0 0 0
0 -4.500000000000000 0 0
72.000000000000000 0 0 9.000000000000000 0
0 -4.500000000000000 0 0 0 33.0000000000
0 0 18.000000000000000 9.000000000000000 0
0 0 0 7.500000000000000 0
0 0 0 0 33.000000000000000
0 0 0 0 0 16.6200000000
sol = 6×6
0.036111111111111 0 0 0
0 0.521460088247092 0 0
0 0 0.088888888888889 0
0 0 0 -0.066666666666667
0 0 0 0
0 0 0 0
y=diag(sol)
y = 6×1
0.036111111111111
0.521460088247092
0.088888888888889
-0.066666666666667
0.036363636363636
0.222623345367028
n=6;
for i=1:n
Xi=X(1,i);
fprintf('\nX%g=',i)
disp(Xi);
end
X1=
0.036111111111111
X2=
0.521460088247092
X3=
16
0.088888888888889
X4=
-0.066666666666667
X5=
0.036363636363636
X6=
0.222623345367028

A=[72 0 0 9 0 0;0 2.88 0 0 0 -4.5;...
0 0 18 9 0 0; 0 0 9 12 0 0;0 0 0 0 33 0;0 -4.5 0 0 0 33]
A = 6×6
72.000000000000000 0 0 9.000000000000000
0 2.880000000000000 0 0
0 0 18.000000000000000 9.000000000000000
0 0 9.000000000000000 12.000000000000000
0 0 0 0
0 -4.500000000000000 0 0
b=[2;0.5;1;0;1.2;5]
b = 6×1
2.000000000000000
0.500000000000000
1.000000000000000
0
1.200000000000000
5.000000000000000

GaussSeidel

1 0.0278 0.1736 0.0556 -0.0417 0.0364 0.1752 0.2603
2 0.0330 0.4473 0.0764 -0.0573 0.0364 0.2125 0.2775
3 0.0349 0.5057 0.0842 -0.0632 0.0364 0.2205 0.0597
4 0.0357 0.5181 0.0871 -0.0653 0.0364 0.2222 0.0131
5 0.0359 0.5207 0.0882 -0.0662 0.0364 0.2225 0.0030
6 0.0360 0.5213 0.0886 -0.0665 0.0364 0.2226 0.0008
7 0.0361 0.5214 0.0888 -0.0666 0.0364 0.2226 0.0002
8 0.0361 0.5215 0.0889 -0.0666 0.0364 0.2226 0.0001
17

Taller

Cargado por

Copyright:

Formatos disponibles

Taller

Cargado por

Información del documento

Descripción original:

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Taller

Cargado por

Copyright:

Formatos disponibles

Taller 1: Sistemas de Ecuaciones

Resolver los ejercicios planteados:

Método de Gauss sin pivoteo

La Matriz C, que corresponte a la matriz aumentada [Ab] es =

0.000006000000000 -0.000002857142857 -0.000003076923077 0

La matriz correspondiente a esta etapa antes del proceso:

0.000006000000000 -0.000002857142857 -0.000003076923077 0

Los Multiplicadores correspondientes a esta etapa son:

La matriz correspondiente a esta etapa despues del proceso:

0.000006000000000 -0.000002857142857 -0.000003076923077 0

La matriz correspondiente a esta etapa antes del proceso:

0.000006000000000 -0.000002857142857 -0.000003076923077 0

La matriz correspondiente a esta etapa despues del proceso:

0.000006000000000 -0.000002857142857 -0.000003076923077 0

0.000006000000000 -0.000002857142857 -0.000003076923077 0

La solucion de X1 hasta Xn es:

Método de Gauss con pivoteo

sol = gausspivoteo(A,b) %script con funcion

0.000008000000000 0.000008571428571 0.000009230769231 1.274000000000000

Método de Gauss Seidel

n=3; %número de ecuaciones

K X(1) X(2) X(3) Norma

Método de Gauss sin pivoteo

La Matriz C, que corresponte a la matriz aumentada [Ab] es =

La matriz correspondiente a esta etapa antes del proceso:

La matriz correspondiente a esta etapa despues del proceso:

La matriz correspondiente a esta etapa antes del proceso:

Los Multiplicadores correspondientes a esta etapa son:

La matriz correspondiente a esta etapa despues del proceso:

La solucion de X1 hasta Xn es:

Método de Gauss con pivoteo

sol = gausspivoteo(A,b) %script con funcion

Método de Gauss Seidel

K X(1) X(2) X(3) Norma

Método de Gauss sin pivoteo

0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000

La matriz correspondiente a esta etapa antes del proceso:

0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000

Los Multiplicadores correspondientes a esta etapa son:

La matriz correspondiente a esta etapa despues del proceso:

0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000

La matriz correspondiente a esta etapa antes del proceso:

0.033300000000000 0.048100000000000 -0.022200000000000 0.050000000000000

Los Multiplicadores correspondientes a esta etapa son:

La matriz correspondiente a esta etapa despues del proceso:

0.003330000000000 0.004810000000000 -0.002220000000000 0.005000000000000

0.003330000000000 0.004810000000000 -0.002220000000000 0.005000000000000

La solucion de X1 hasta Xn es:

Método de Gauss con pivoteo

sol = gausspivoteo(A,b) %script con funcion

0.048100000000000 0.100000000000000 7.450000000000000 0

Método de Gauss Seidel

n=3; %número de ecuaciones

K X(1) X(2) X(3) Norma

% A=[72 0 0 9 0 0;0 2.88 0 0 0 -4.5;...

La Matriz C, que corresponte a la matriz aumentada [Ab] es =

La matriz correspondiente a esta etapa antes del proceso:

La matriz correspondiente a esta etapa despues del proceso: