APPUNTI MATLAB CALCOLO NUMERICO
%metodo sostituzione all'indietro
x=length(b)
x(n)=b(n)/A(n,n)
s=0
for i=n-1:-1:1
s=A(i,i+1:n)*x(i+1:n)
end
x(i)=(b(i)-s)/A(i,i)
end
%GAUSS
function[x]=gauss_noscambi(A,b)
n=length(b)
for k=1:n-1
for i=k+1:n
A(i,k)=A(i,k)-A(i,k)*A(k,k)
end
for j=k+1:n
b(i)=b(i)-A(i,k)*b(k)
end
end
x=zeros(n,1)
x(n)=b(n)/A(n,n)
for i=n-1:-1:1
x(i)=b(i)-A(i,i+1:n)*x(i+1:n)/A(i,j)
end
%FATTORIZZAZIONE PA=LU (risoluzione SISTEMI LINEARI)
[L U P]=lu(A)
y=L\(P*b)
x=U\y
%calcolo dell'inversa
Ainv=Linv*Uinv*P
%FATTORIZZAZIONE DI CHOLESKI A=R'R (risoluzione SISTEMI
LINEARI)
