Bonjour,

j'ai un compte rendu de TP en analyse numérique matricielle à faire sur les méthodes de gradient je dois faire un algorithme à pas fixe.
J'ai réalisé cet algorithme qui fonctionnait très bien en labo ce matin mais quand l'ai testé chez moi (version 64bits) ça ne fonctionne plus alors que je n'ai fait qu'un copier-coller...
Quelqu'un peut m'aider et me dire ce qu'il ne va pas?


function [x,n,resi]=GPF(A,b,x0,eta,Imax,rho)
    resi=zeros(1,1);
    x=x0;
    resi(1,1)=norm(A*x0-b)
    nr =norm(A*x-b)
    n=1
    while (n<=Imax) & (nr>eta)
        n=n+1;
        x=x -rho*(A*x-b);
        resi(n,1)=norm(A*x-b);
        nr=norm(A*x-b)
    end  
endfunction

E=diag(ones(19,1),1)
A=21*21*(2*eye(20,20)-E-E')
b=ones(20,1)
x0=zeros(20,1)
n=0
//rho=2/max(spec(A))

Merci d'avance à mon(ma) sauveur(se)