untitled.m

这是用来解线性方程的

A=[1 2 1;2 2 3;-1 -3 0];b=[1 2 3];
[L,U]=lu(A);
[x1,L1,U1,P1]=col_gauss(A,b,0);
B=[1 2 6;2 5 15;6 15 46];
[L1,U1]=lu(B);
C=[1 0 0 ;2 1 0 ;-1 -1/2 1];
D=[1 2 1;0 2 -1;0 0 1];
C*D;
[L,U]=luNopiv(A,50)
[L,U,pv]=luPiv(A,0.0001);