Скачиваний:
329
Добавлен:
04.03.2014
Размер:
1 Кб
Скачать
program kray;
uses crt;
 const n=10;
  var   a,b,c,f:array[1..10] of real;
         alfa, beta ,y,yan:array[0..10] of real;
          i,j,kapa1,kapa2,p,q,a1,b1:integer;
           h,aa,bb,mu1,mu2,x:real;
   begin
    clrscr;
    a1:=0;b1:=1;
     h:=(b1+a1)/n;
      aa:=1;bb:=1/4*exp(4)+1/3*exp(3)+5/12;
       p:=-7;q:=12;
   for i:=1 to n-1 do begin
      a[i]:=1/(h*h)-p/(2*h);
       b[i]:=-2/(h*h)+q;
        c[i]:=1/(h*h)+p/(2*h);
         f[i]:=5;    end;
    kapa1:=0 ;

    kapa2:=0;
      mu1:=aa;              mu2:=bb;
    alfa[1]:=kapa1;
     beta[1]:=mu1;
      for i:=1 to n-1 do
       begin
        alfa[i+1]:=-c[i]/(a[i]*alfa[i]+b[i]);
         beta[i+1]:=(f[i]-a[i]*beta[i])/(a[i]*alfa[i]+b[i]);end;
          y[n]:=(kapa2*beta[n]+mu2)/(1-kapa2*alfa[n]);
           for i:=n-1 downto 0 do
             y[i]:=alfa[i+1]*y[i+1]+beta[i+1];

     for i:=0 to n do begin
      x:=a1+i*h;
       yan[i]:=1/4*exp(4*x)+1/3*exp(3*x)+5/12;end;

     writeln('  ','X','         ','YAN','           ','y','  ');
    for i:=0 to n do
     writeln((a1+i*h):3:1,'      ',yan[i]:7:2,'     ',y[i]:7:2);
        readln;
    end.







Соседние файлы в папке Alex