выч методы не разобраны! / Лаба№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.