Добавил:
bagiwow
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз:
Предмет:
Файл:mettod cetok
.docProgram 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:=(x+0.5)*(x+0.5); end; Function Ft(x:real):real; begin ft:=(x+1)*sin(x); end; Function Q(t:real):real; begin q:=0.5*(0.5+t); end; Function G(t:real):real; begin g:=2.25; 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-2 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.
Соседние файлы в предмете Вычислительные методы