Ciao a tutti. Devo implementare una function che mi risolve questo tipo di matrici, ecco come l'ho fatta, non capisco dove sbaglio, ma non dà il giusto risultato, neanche il primo termine del vettore x che dovrebbe essere giusto sempre per come l'ho scritto. Qualcuno mi può dare una mano? Grazie mille
function [x] = triangolareinf( A,b )
%Metodo delle sostituzioni successive in avanti per la risoluzione di un
%sistema triangolare inferiore.
%IMPUT
% A--> Matrice dei coefficienti
% b--> vettore termini noti
%OUTPUT
% x--> soluzione del sistema
[n,m]=size(A);
if n~=m
error ('la matrice deve essere quadrata')
end
for i= 1:n
if A(i,i)==0
error ('la matrice è singolare')
end
end
x(1) = b(1) / A(1,1);
c = 1;
t=0;
for j= 2:n
for l= 1:c
s = t + (A(j,l).*x(l));
t = (A(j,l).*x(l));
end
t=0;
x(j) = ( b(j) - s ) / A(j,j);
c = c+1;
end
end