МММ / ммм / упр4

tau=0.2;

N=5/tau;

for n=1:(N-1)

tn=n*tau;

U(n)=(1+tn)*exp(-tn);

U(n+1)=(1-tau*tn/(1+tn))*U(n)

plot(tn,U(n+1),'r*')

hold on

end

for k=1:N

tk=k*tau;

U(k)=(1+tk)*exp(-tk);

U(k+1)=((1+tk)/(1+(tau+1)*tk))*U(k);

plot(tk,U(k+1),'bo')

hold on

end

for r=2:(N-1)

tr=r*(1-tau);

U(r)=(1+tr)*exp(-tr);

U(r+1)=U(r-1)-2*U(r)*tr*tau/(1+tr);

plot(tr,U(r+1),'g*')

hold on

end