3_теориЛеготкина / выч методы / метод сеток
.docxProgram Metod_Setok;
Const h=0.1; a=0; b=1; n=round((b-a)/h);
Var result :text;
arr :array [0..n, 0..n] of real;
i,k :integer;
x,t :real;
Function F(x:real):real;
begin
f:=0.5*x*(x+1);
end;
Function Ft(x:real):real;
begin
ft:=x*cos(x);
end;
Function Q(t:real):real;
begin
q:=2*t*t;
end;
Function G(t:real):real;
begin
g:=1;
end;
Begin
for i:=0 to n do begin
x:=a+i*h;
arr[i,0]:=F(x);
end;
for k:=0 to n do begin
t:=a+k*h;
arr[0,k]:=Q(t);
end;
for k:=0 to n do begin
t:=a+k*h;
arr[10,k]:=G(t);
end;
for i:=1 to n-1 do begin
x:=a+i*h;
arr[i,1]:=arr[i,0]+h*Ft(x);
end;
for k:=2 to n do begin
for i:=1 to 9 do begin
arr[i,k]:=arr[i+1,k-1]+arr[i-1,k-1]-arr[i,k-2];
end;
end;
assign(result, 'Table_of_Result.txt'); (сюда запишется результат)
rewrite(result);
write (result, ' i ');
for i:=0 to n do begin
write(result, i:8);
end;
writeln(result);
writeln(result,'k');
for k:=0 to n do begin
write(result,' ');
write(result, k:2);
for i:=0 to 10 do begin
write(result, arr[i,k]:8:3);
end;
writeln(result);
end;
close(result);
readln;
End.