Script y Ejemplo Para Matlab Para Resolver Una Matriz Por El Metodo De Gauss

function [b]= gauss(A, b)
[m n]= size (A);

for p=1:n
d = A(p,p);
for c=1:n
A(p,c)= (A(p,c))/d;
end
b(p)= b(p)/d;
for r=1:n
k= A(r,p);

if r ~= p
for c=1:n
A(r,c)=A(r,c)-(k* (A(p,c)));
end
b(r)=b(r)-(k* (b(p)));

end

endend

end

>> A= [415583.301 179878.041 0 0 0 0; 149878.041 1055343.579 -8081.22 168358.757 0 0; 0 -8081.22 209.730 -5675.944 2405.277 0; 0 168358.757 -5675.944 810938.797 69382.142 0; 0 0 2405.277 69352.142 212214.93 310.060; 0 0 0 0 310.060 41039.411]

A =

1.0e+006 *

0.4156 0.1799 0 0 0 0
0.1499 1.0553 -0.0081 0.1684 00
0 -0.0081 0.0002 -0.0057 0.0024 0
0 0.1684 -0.0057 0.8109 0.0694 0
0 0 0.0024 0.0694 0.2122 0.0003
0 0 0 0 0.0003 0.0410

>> b= [-110; 75; -1.5; 2.5; 144; 96]

b =

-110.0000
75.0000
-1.5000
2.5000
144.0000
96.0000

>> resultado =gauss(A,b)

resultado =

-0.0002
-0.0001
-0.0317
-0.0003
0.0011
0.0023
>> A= [70 0 0 0 -35 0 0 0 0 0 0 0 0;0 105 0 -105 0 0 0 0 0 0 0 0 0;0 0 201.61 80.64 -40.32 26.88 -80.64 -53.76 0 0 0 0 0; 0 -105 80.64 194.61 26.88 -17.92 -53.76 -35.84 0 0 0 0 0; -35 0 -40.32 26.88 150.64 0 0 0 -35 0 -40.32 -26.88 0; 0 0 26.88 -17.92 0 88.34 0 -52.5 0 0 -26.88 -17.92 0; 0 0 -80.64-53.76 0 0 161.29 0 0 0 -80.64 53.76 0; 0 0 -53.76 -35.84 0 -52.5 0 124.18 0 0 53.76 -35.84 0; 0 0 0 0 -35 0 0 0 70 0 0 0 -35; 0 0 0 0 0 0 0 0 0 105 0 -105 0; 0 0 0 0 -40.32 -26.88 -80.64 53.76 0 0 201.61 -80.64 -80.64; 0 0 0 0 -26.88 -17.92 53.76 -35.84 0 -105 -80.64 194.61 53.76; 0 0 0 0 0 0 0 0 -35 0 -80.64 53.76 115.64]

A =

Columns 1 through 3

70.0000 0 00 105.0000 0
0 0 201.6100
0 -105.0000 80.6400
-35.0000 0 -40.3200
0 0 26.8800
0 0 -80.6400
0 0 -53.7600
0 0 0
0 0 0
0 0 0
0 0 0
0 0 0

Columns 4 through 60 -35.0000 0
-105.0000 0 0
80.6400 -40.3200 26.8800
194.6100 26.8800 -17.9200
26.8800 150.6400 0
-17.9200 0 88.3400
-53.7600 0 0
-35.8400 0 -52.5000
0 -35.0000 0
0 0 0
0 -40.3200 -26.8800
0 -26.8800 -17.9200
00 0

Columns 7 through 9

0 0 0
0 0 0
-80.6400 -53.7600 0
-53.7600 -35.8400 0
0 0 -35.0000
0 -52.5000 0
161.2900 0 0
0 124.1800 0
0 0 70.0000
0 0 0
-80.6400 53.7600 053.7600 -35.8400 0
0 0 -35.0000

Columns 10 through 12

0 0 0
0 0 0
0 0 0
0 0 0
0 -40.3200 -26.8800
0 -26.8800 -17.9200
0 -80.6400 53.7600
0 53.7600 -35.8400
0 0 0
105.0000 0-105.0000
0 201.6100 -80.6400
-105.0000 -80.6400 194.6100
0 -80.6400 53.7600

Column 13

0
0
0
0
0
0
0
0
-35.0000
0
-80.6400
53.7600
115.6400

>> b= [0;0;0;-4;0;0;0;-4;0;0;0;-4;0]

b =

0
0
0
-4
0
0
0
-4
0...