ЧМ (ЭКТ-3) / Лабы / Всякие сделанные варианты / other / lab6

M=1

T=1

a=1

t=0

x=0

step_by_x=6

step_by_t=6

h=a/step_by_x

tau=T/step_by_t

U1=zeros(step_by_t+1,step_by_x+1)

%Вычислим x

for i=1:step_by_x+1

x(1,i)=((i-1)*a)/step_by_x

end

%Вычислим t

for i=1:step_by_t+1

t(i,1)=((i-1)*T)/step_by_t

end

%Вычислим U(x,0)

U1(1,:)=x(1,:).^2+1

%Вычислим U(0,t)

U1(:,1)=t(:,1).^3+1

U2=U1

U3=U1

for j=2:step_by_t+1

for i=2:step_by_x+1

tt=t(j-1,1)+h/2

xx=x(1,i)-tau/2

fi_i = (3*(tt^2))+( M*2*xx)

U1(j,i) = (fi_i-(M/h)*(U1(j-1,i)-U1(j-1,i-1))+(U1(j-1,i)/tau))*tau

U2(j,i) = (fi_i - ((1/tau)-(M/h))*U2(j,i-1)+(1/tau)*U2(j-1,i-1))*(h/M)

U3(j,i) = (fi_i+(M/h)*(U3(j,i-1))+(1/tau)*(U3(j-1,i)))/((M/h)+(1/tau))

end

end

Y=0

for j=1:step_by_t+1

for i=1:step_by_x+1

Y(j,i) = t(j,1)^3+x(1,i)^2+1

end

end

mesh(x,t,U1)

pause

mesh(x,t,U2)

pause

mesh(x,t,U3)

pause

mesh(x,t,Y)

pause

mesh(x,t,U1-Y)

pause

mesh(x,t,U2-Y)

pause

mesh(x,t,U3-Y)

pause