Examen practicas

158 days ago by Baston

 
       

b=matrix(QQ,3,1,[1,1,1])
A=matrix(QQ,3,3,[1,4,2,5,2,1,1,1,3])
P,L,U= mLU(A)
show(html("$"+latex(P)+latex(A)+"="+latex(L)+latex(U)+"$"))
#caculamos Ly=Pb
y=bajada(L,P*b)
show(html("$\\text{y = }" +latex(y.n()) + "$"))
#calulamos Ux=y
x=subida(U,y)
show(html("$\\text{x = }" +latex(x.n()) + "$"))
#solución SAGE
sol=A\b
show(sol.n())
show("error relativo: " + ((sol-x).norm(oo)/sol.norm(oo)).str())