выч методы не разобраны! / Лаба№3 / Мишины / PROGON

Program m_progonki;
Uses Crt;
const h=0.05; n=6;
a=0.2; b=0.5;
a0=1; a1=0;
b0=0.5; b1=-1;
a3=2; b3=1;
type
massiv=array[0..n] of real;
Var
x,y,c,d:massiv;
i:integer;
out:text;
function p(x :real):real;
begin
p:=2*x
end;

function q(x :real):real;
begin
q:=-1/(x*x)
end;

function f(x :real):real;
begin
f:=3+x
end;

function m(x:real):real;
begin
m:=(2*p(x)*h)/(2*sqr(h)*q(x)-4)
end;

function k(x:real):real;
begin
k:=(2-h*p(x))/(2*sqr(h)*q(x)-4)
end;

function f1(x:real):real;
begin
f1:=2*h*h*f(x)/(2*sqr(h)*q(x)-4)
end;

{-----------------------MAIN PROGRAM--------------------------}

Begin
clrscr;
assign(out,'progon.dat');
rewrite(out);
writeln(out,'ЊҐв®¤ Їа®Ј®­ЄЁ : ');
writeln('ЊҐв®¤ Їа®Ј®­ЄЁ : ');
writeln;
x[0]:=a;x[n]:=b;
c[0]:=(m(x[0])*(a0*h-a1))/(a0*h-a1*(1+k(x[0])));
d[0]:=(f1(x[0])*(a0*h-a1)-k(x[0])*a3*h)/(m(x[0])*(a0*h-a1));
writeln(out,'‚室­лҐ ¤ ­­лҐ : ');
writeln( '‚室­лҐ ¤ ­­лҐ : ');
writeln;
writeln(out,' ®в१®Є [',a:1:1,',',b:1:1,'] и Ј h=',h:1:2);
writeln( ' ®в१®Є [',a:1:1,',',b:1:1,'] и Ј h=',h:1:2);
writeln(out,' a0=',a0,' a1=',a1,' a3=',a3);
writeln( ' a0=',a0,' a1=',a1,' a3=',a3);
writeln(out,' b0=',b0:1:1,' b1=',b1,' b3=',b3);
writeln( ' b0=',b0:1:1,' b1=',b1,' b3=',b3);
writeln(out,' p(x)=x q(x)=-1/(x*x) f(x)=3+x ');
writeln( ' p(x)=x q(x)=-1/(x*x) f(x)=3+x ');

for i:=1 to n-1 do
begin
x[i]:=x[i-1]+h;
c[i]:=m(x[i])/(1-k(x[i])*c[i-1]);
d[i]:=(f1(x[i])-k(x[i])*c[i-1]*d[i-1])/m(x[i])
end;

y[n]:=(b3*h+b1*c[n-1]*d[n-1])/(b0*h+b1*(c[n-1]*d[n-1]+1));
for i:=n-1 downto 0 do
y[i]:=c[i]*(d[i]-y[i+1]);

y[0]:=(a3*h-a1*y[1])/(a0*h-a1);
writeln;
writeln(out,'‚л室­лҐ ¤ ­­лҐ :');
writeln( '‚л室­лҐ ¤ ­­лҐ :');
writeln;
for i:=0 to n do
begin
writeln(out, 'x[',i,']=',x[i]:1:4,' y[',i,']=',y[i]:1:4);
writeln( 'x[',i,']=',x[i]:1:4,' y[',i,']=',y[i]:1:4);
end;
close(out);
readln
end.