Matlab - SUSTITUCIONES REGRESIVAS MATLAB

function [x]=regresivas(A,b)

[f,c]=size(A);

if f~=c

    error('la matriz de coeficientes no es cuadrada')

end

x=zeros(f,1);

for i=f-1:-1:1

    x(i)=b(i);

    for j=i+1:f

        x(i)=x(i)+A(i,j)*x(j);

    end

    x(i)=x(i)/A(i,i);

end

function [x]=regresivas(A,b)

[f,c]=size(A);

if f~=c

    error('la matriz de coeficientes no es cuadrada')

end

x=zeros(f,1);

for i=f-1:-1:1

    x(i)=b(i);

    for j=i+1:f

        x(i)=x(i)-A(i,j)*x(j); %es aquí donde me refiero que deberías tener un menos, y no un más.

    end

    x(i)=x(i)/A(i,i);

end

function x=regresivas_Ax_b(A,b)

if nargin==0

    A=[1 4 5; 0 3 2; 0 0 3]

    b=[1; 2; 6]

end

[f,c]=size(A);

if f~=c

    error('la matriz de coeficientes no es cuadrada')

end

x=zeros(f,1);

for i=f:-1:1

    x(i)=b(i);

    for j=i+1:f

        x(i)=x(i)-A(i,j)*x(j);

    end

    x(i)=x(i)/A(i,i);

end

>> A=[1 4 5; 0 3 2; 0 0 3],b=[1; 2; 6],x=regresivas_Ax_b(A,b)

A =

     1     4     5

     0     3     2

     0     0     3

b =

     1

     2

     6

x =

   -6.3333

   -0.6667

    2.0000