youqi.m

这是用来解线性方程的

A=[80 0 30 10;0 80 10 10;16 20 60 72;4 0 0 8];
b=[40 27 31 2];
[x,L,U,P]=col_gauss(A,b,0.1)
[L1,U1]=lu(A)
[L2,U2]=luNopiv(A,0.1)