1.程序代码 function x=majacobi(A,b,x0,ep,N) n=length(b); if nargin<5 N=500; end if nargin<4 ep=1e-6; end if nargin<3