//Example 3.16 //Jacobi Method //Page no. 72 clc;clear;close; x0=0;y0=0;z0=0; deff('x=f1(y,z)','x=(y-z+10)/5') deff('y=f2(x,z)','y=(-2*x+z+11)/8') deff('z=f3(x,y)','z=(x-y+3)/4') for i=1:13 x1=f1(y0,z0); y1=f2(x0,z0); z1=f3(x0,y0); printf('\tx(%i) = %g\n\n\ty(%i) = %g\n\n\tz(%i) = %g\n\n\n\n',i,x1,i,y1,i,z1) x0=x1;y0=y1;z0=z1; end printf('Thus we find that solution converges to %g,%g and %g',x0,y0,z0)