Array1D< real > z(l.dim());
int i;
int N=l.dim();
z[0]=t[0]/l[0];
for(i=1 ; i<N ; i++){
z[i]=t[i]-(m[i-1]*z[i-1]);
z[i]=z[i]/l[i];
}
return z;
}
/** HACER ALUMNO */
int i;
int N=z.dim();
Array1D< real >v(N);
v[N-1]=z[N-1];
for(i=N-2;i>=0;i--){
v[i]=z[i]-u[i]*v[i+1];
}
return v;
}
/** HACER ALUMNO */
Array1D< real > l(a.dim());
Array1D< real > m(c.dim());
Array1D< real > u(b.dim());
int aux = crout_descomposicion(a,b,c,l,m,u);
if (aux==-1)return Array1D< real > ();
Array1D< real > z = crout_descenso(l,m,t);
Array1D< real > v = crout_remonte(u,z);
retu