Matlab - eliminación gaussiana

Hola, tengo un problema para dejar una matriz diagonal. Soy capaz de sacar la matriz triangular superior pero de ahí en adelante tengo algo mal en el código que no sé qué es. Por favor, alguien que me diga que tengo que corregir.

Este es el código que estoy utilizando








row = input('Ingrese número de filas ');
column = input ('Ingrese número de columnas ');



for a = 1 : row
for b = 1 : column

valor = input('Ingrese valor de constante ');
matriz(a , b) = valor;
end
end
disp('La matriz es')
matriz


for c = 1 : row - 1
for f = 1 : column
if(c == f)
if matriz (c , f) == 0 && matriz (c + 1 , f) ~= 0
for a = c : row - 1
for b = 1 : column

arreglo2(b) = matriz(a , b);
arreglo3(b) = matriz(a + 1 , b);
end
end

for a = c : row - 1
for b = 1 : column

matriz(c , b) = arreglo3(b);
matriz(a + 1 , b) = arreglo2(b);
end
end

sprintf('Hubo un cambio entre la fila %d y %d', c, c + 1)
end
for a = c : row - 1
for b = 1 : column

arreglo(b) = matriz(c , b);
arreglo1(b) = matriz(a + 1 , b);
end

d = arreglo(c);
e = arreglo1(c);

for b = 1 : column

arreglo(b) = (arreglo(b)/d)*e;
arreglo1(b) = arreglo1(b) - arreglo(b);
end

for b = 1 : column
matriz(a + 1 , b) = arreglo1(b);
end
end
end
end
end

disp('La matriz con pivote parcial es')
matriz


for c = row : -1 : 1
for f = column : -1 : 1
if(c == f)
if( matriz(c , f) == 0)
continue;


for a = c: -1 : 1
for b = 1 : column

arreglo(b) = matriz(c , b);
arreglo1(b) = matriz (a - 1 , b);
end

d = arreglo(c);
e = arreglo1(c);

for b = 1 : column

arreglo(b) = (arreglo(b)/d)*e;
arreglo1(b) = arreglo1(b) - arreglo(b);
end

for b = 1 : column
matriz (a - 1 , b) = arreglo1(b);

end
end
end
end
end
end



disp('La matriz con pivote total es')
matriz



clear

1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

17

18

19

20

21

22

23

24

25

26

27

28

29

>> eliminacion_gaussiana1_foro

La matriz es

matriz =

    0.3400    0.8467    0.2461    0.5815    0.9377

    0.0478    0.0540    0.0206    0.6815    0.5986

    0.1140    0.7962    0.6179    0.0702    0.0693

    0.1360    0.7889    0.0924    0.2379    0.2436

    0.1048    0.8584    0.6982    0.7337    0.6505

La matriz con pivote parcial es

matriz =

    0.3400    0.8467    0.2461    0.5815    0.9377

         0   -0.0650   -0.0140    0.5997    0.4668

         0         0    0.4253    4.5993    3.4320

         0         0         0    5.2679    3.9291

         0         0         0         0    0.1253

matriz =

    0.3400         0         0         0         0

         0   -0.0650         0         0         0

         0         0    0.4253         0         0

         0         0         0    5.2679         0

         0         0         0         0    0.1253