ВТМО - Введение в Теорию Моделирования и Оптимизации

a=1.1; b=1.8; d=2.5; l=2; h=0.05;

s=0.16666;

n=1/h+1;

x=0:h:1;

c1=(4/l^2)*(-d+(a+b)/2);

c2=(2/l)*(2*d-(3*a+b)/2);

u=c1*x.^2+c2*x+a*ones(1,n);

v=u;

u0=u;

m=25;

for j=1:m,

for i=2:(n-1),

v(i)=u(i)+s*(u(i-1)-2*u(i)+u(i+1));

end;

u=v;

end;

u1=u;

for j=1:125,

for i=2:(n-1),

v(i)=u(i)+s*(u(i-1)-2*u(i)+u(i+1));

end;

u=v;

end;

u2=u;

for j=1:100,

for i=2:(n-1),

v(i)=u(i)+s*(u(i-1)-2*u(i)+u(i+1));

end;

u=v;

end;

u3=u;

for j=1:250,

for i=2:(n-1),

v(i)=u(i)+s*(u(i-1)-2*u(i)+u(i+1));

end;

u=v;

end;

u4=u;

plot(x,u0,x,u1,x,u3,x,u4)

grid

prtsc