5 réponses
Il te faut trois fichiers:
main.m qui contient:
%programme principal pour résoudre l'equation Ax=b
n=input ('Donner la taille de la matrice, n = ');
n
b=input('Donner le second membre, b = ');
b
A=input('Entrer la matrice, A = ');
A
[A,b]=triangle(A,b,n);
disp('La solution du système est')
b = solve_triangle(A,b,n)
triangle.m qui contient:
function [A,b]=triangle(A,b,n)
for k=1:n
for i=k+1:n
A(i,k)=A(i,k)/A(k,k);
b(i)=b(i)-A(i,k)*b(k);
for j=k+1:n
A(i,j)=A(i,j)-A(i,k)*A(k,j);
end
end
end
solve_triangle.m qui contient:
function b=solve_triangle(A,b,n)
for k=n:-1:1
for j=k+1:n
b(k)=b(k)-A(k,j)*b(j);
end
b(k)=b(k)/A(k,k);
end
main.m qui contient:
%programme principal pour résoudre l'equation Ax=b
n=input ('Donner la taille de la matrice, n = ');
n
b=input('Donner le second membre, b = ');
b
A=input('Entrer la matrice, A = ');
A
[A,b]=triangle(A,b,n);
disp('La solution du système est')
b = solve_triangle(A,b,n)
triangle.m qui contient:
function [A,b]=triangle(A,b,n)
for k=1:n
for i=k+1:n
A(i,k)=A(i,k)/A(k,k);
b(i)=b(i)-A(i,k)*b(k);
for j=k+1:n
A(i,j)=A(i,j)-A(i,k)*A(k,j);
end
end
end
solve_triangle.m qui contient:
function b=solve_triangle(A,b,n)
for k=n:-1:1
for j=k+1:n
b(k)=b(k)-A(k,j)*b(j);
end
b(k)=b(k)/A(k,k);
end
