Solution of Possion2D equation with method ofJacobi iteration.
%setparameterN1=31;N2=41;L1=1000;L2=1500;rp=1.0;%relaxationparameterdeltax=L1/(N1-1)deltay=L2/(N2-1)figure(1)subplot(1,2,1)output_interval=10%initializegrids[x,y]=meshgrid(0:L2/(N2-1):L2,0:L1/(N1-1):L