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

a=0.5

b=1.5

k1=1

k2=0.3

k3=0

k4=1

c=0.6

d=1.8

t=abs(b-a)/50;

A(1)=0;

B(1)=k1-k2/t;

c(1)=k2/t;

t(1)=a;

for i=2:51

T(i)=a+t*i

A(i)=1/t^2-1/(2*t*T(i))

B(i)=-0.4-2/t^2

C(i)=1/t^2+1/(2*t*T(i))

R(i)=2*T(i)

end;

A(51)=-k4/t;

B(51)=k3+k4/t;

C(51)=0;

al(1)=0;

be(1)=0;

for j=2:51

al(j)=-C(j-1)/(B(j-1)+A(j-1)*al(j-1))

be(j)=(R(j-1)-A(j-1)*be(j-1))/(A(j-1)*al(j-1)+B(j-1))

end;

y(51)=(R(51)-A(51)*be(51))/(B(51)+A(51)*al(51));

j=51;

while j>1

y(j-1)=al(j)*y(j)+be(j)

j=j-1

end;

plot(T,y)