数值线性代数课程设计 1、 >> p=eig(redgongwithmass)； >> p=max(abs(p)) p = 2.0414 >> q=max(abs(eig(redgongwithstiff))) q = 9.8915e+009 2、 function [l,m]=mineig(A,epsilon) %按模最小 u0=diag(eye(size(A))); a=norm(u0); y=u0/a; A=inv(A); u=A*y;b0=y'*u;b=3*b0; while abs((1/b-1/b0)/(1/b))>epsilon u0=u; b0=b; a=nor