= x + (b-a)/n
begin
y: = 0; x: = a;
y = y * (b-a)/n
for i: = 0 to n do
begin
y: = y + f (x);
x: = x + (b-a)/n;
end;
y: = y * (b-a)/n;
writeln ('y =', y: 8:5);
readln;
end.
ВІДПОВІДЬ:
y = 0.28099
В
у = 0; х = а
PROGRAM TRAPEZIA;
h = (b-a)/n
x = a + h
CONST a = 0.4; b = 0.8; n = 100;
var x, y, h, s: real;
s: = (f (a) + f (b))/2
n: integer;
function f (x: real): real;
begin
s = s + f (x)
f: = (sin (x * x +0.5)/cos (x * x +0.5))/(1 +2 * x * x);
end;
x: = x + h
begin
h: = (b-a)/n;
y: = s * h
x: = a + h;
s: = (f (a) + f (b))/2;
while x <= b-h do
begin
s: = s + f (x);
x: = x + h;
end;
y: = s * h;
writeln ('y =', y: 8:5);
readln;
end.
ВІДПОВІДЬ:
y = 0.28173
В В
h = (b-a)/2n
В
PROGRAM SIMPSON;
s: = f (a) + f (b)
LABEL 10;
CONST a = 0.4; b = 0.8; n = 100;
c: = 1
x: = x + h
var x, y, h, s: real;
c, n: integer;
s: = s + (3 + c) * f (x)
function f (X: real): real;
begin
x: = x + h
f: = (Sin (x * x +0.5)/cos (x * x +0.5))/(1 +2 * x * x);
c =-c
end;
begin
h: = (b-a)/(2 * n);
s: = f (a) + f (b);
y: = s * h/3
c: = 1;
x: = a + h;
10: s: = s + (3 + c) * f (x);
x: = x + h;
c: =-c;
if x <= b-h then goto 10;
y: = s * h/3;
writeln ('y =', y: 8:5);
readln;
end.
ВІДПОВІДЬ:
y = 0.27919
Рішення інтеграла в ППП "Eureka"
y = integ ((sin (x ^ 2 +0.5)/cos (x ^ 2 +0.5))/(1 +2 * x ^ 2), x, 0.4,0.8)
y = 0.2823564
2.4 Методи рішення диференціальних рівнянь
При використанні різних протікають у часі процесах першим етапом є складання диференціального рівняння, описує процес, а другим - пошук рішення цього рівняння. Диференціальним рівнянням називається рівність, що пов'язує значення аргументу, невідомої функції деяких її похідних. Найвищий порядок входять в рівняння похідних називається порядком диференціального рівняння.
