2 Симпсон әдісінің программа листингі
program SIMPSON;
uses crt;
label 1,2,3;
var
n,i,j:integer;
a,b,eps,y,x,h,s,k,f1,t:real;
function f(x:real):real;
begin
f:=(sin(x*x)/cos(x*x))/(x*x+1);
end;
begin
window(1,1,80,25);
textbackground(1);
clrscr;
writeln('****************************************************');
writeln('* *');
writeln('* Данной программой вычисляется *');
writeln('* определенный интеграл *');
writeln('* методом Стмпсона *');
writeln('* *');
writeln('****************************************************');
writeln;
window(15,5,65,15);
textbackground(0);
clrscr;
textcolor(13);
writeln('Введите a,b,n,eps');
readln(a,b,n,eps);
h:=(b-a)/n;
s:=0;
k:=n/2;
j:=trunc(k);
1: f1:=f(a);
y:=f1;
f1:=f(b);
y:=(y-f1)/2;
x1:=a+h;
for i:= 1 to j do
begin
f1:=f(x1);
y:=y+2*f1;
f1:=f(x1+h);
y:=y+f1;
x1:=x1+2*h
end;
y:=y*2*h/3;
t:=abs(y-s);
s:=y;
h:=h/2;
j:=j*2;
if t>eps then goto 1;
window(15,20,65,23);
textbackground(6);
clrscr;
textcolor(13);
writeln('Значение интеграла I=',y:7:4);
readln
end.
4 қосымша
1 Қарапайым дифференциалдық теңдеулерді
сандық әдістермен шешу
Оқу құралында негізінен бірінші ретті кәдімгі дифференциал-дық теңдеулер үшін Коши есебі қарастырылады:
(9)
Теңдеудің оң жағы, екі айнымалыдан тәуелді функция – қандайда бір облыста анықталған деп есептейміз.
Достарыңызбен бөлісу: |