ЧМ (ЭКТ-3) / Лабы / Лаба 7 (ДЗ) / Laba#7 ВАР№8

Программа:

n=50;

h=1/n;

x=0:h:1; y=0:h:1;

u=zeros(n);

u(n/2:n,1)=0; u(n/2:n,n)=0; u(n,:)=0;

for k=1:n/2,

u(k,1)=10*exp((k+1)*0.01), u(k,n)=10*exp((k+n)*0.01);

end;

for m=2:n ,

u(1,m)=10*exp((1+m)*0.01);

end;

f=zeros(n);

for i=1:1:n,

for j=1:1:n,

f(i,j)=0;

end;

end;

mesh(f);

%pause;

eps=0.001

pogr=1000;

u0=zeros(n);

while (pogr>eps),

pogr=0;

for i=2:1:n-1,

for j=2:1:n-1,

u(i,j)=(u(i+1,j)+u(i-1,j)+u(i,j+1)+u(i,j-1))/4;

pogr0=abs(u(i,j)-u0(i,j));

if(pogr0>pogr),

pogr=pogr0;

end;

u0(i,j)=u(i,j);

end;

end;

end;

mesh(u);