sust progresiva

1137 days ago by etlopez

%auto def sustitucion_progresiva(L_aumentada): n=L_aumentada.nrows() x=vector(RDF,n) for i in [0,..,(n-1)]: x[i]=L_aumentada[i,n]; for j in [0..i-1]: x[i]=x[i]-L_aumentada[i,j]*x[j] x[i]=x[i]/L_aumentada[i,i] return x