untitled4.m

这是用来解线性方程的

A=[5 10 15 20 25;7 21 42 63 84;6 16 32 56 66;6 17 37 76 96;7 19 39 73 93];
b=[15 49 24 27 33];
[x,L,U,P]=col_gauss(A,b,0.0001)
[L1,U1]=lu(A);
[L2,U2]=luNopiv(A,0.1);