Скачиваний:
289
Добавлен:
04.03.2014
Размер:
1.18 Кб
Скачать
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