Название: Вычислительная математика
Автор: Валерий Федорович Альмухаметов
Издательство: ЛитРес: Самиздат
Жанр: Руководства
isbn:
isbn:
WRITELN (‘ИНТЕГРАЛ РАВЕН ’, S);
END.
Подпрограмма на языке Basic
10 DEF FNA(X)=EXP(-Х*Х)
20 PRINT "N,B0,B9,H1"; : INPUT N,B0,B9,H1
30 C=2/SQR(PI) : A=0 : S1=0
40 FOR B=B0 TO B9 SТЕР H1
50 GOSUB 100
60 S1=S1+S : A=B
70 PRINT B,C*S1 : NEXT В
90 GOTO 10
100 H=(B-A)/N : S1=(FNA(A)+FNA(B))/2
110 FOR I = 1 TO N-1 : S=S+FNA(A+ I *H) : NEXT I
120 S=S*H
190 RETURN
110 REM МЕТОД ТРАПЕЦИИ
120 INPUT . “Входные переменные A,B,N =”;A,B,N
130 H = (B-A)/N
140 S = O
150 X = A
160 FOR I = 1 TO N-1
170 X = X+H
180 S = S + FNY (X)
190 NEXT I
200 S = H*(FNY(A) + FNY(B) + 2*S) / 2
210 RETURN
Подпрограмма на языке Pascal
VAR N,I,K:INTEGER; A,B,B0,B9,H,C,S,S1:REAL;
FUNCTION F(X:REAL):REAL;BEGIN
F:=EXP(-X*X);END;
PROCEDURE TRAP(VAR A,B:REAL;
N:INTEGER; FUNCTION F:REAL;S:REAL);
VAR I:INTEGER;H:REAL;
BEGIN H:=(B-A)/N;S:=(F(A)+F(B))/2;
FOR I:=1 TO N-1 DO S:=S+F(A+I*H);
S:=S*H;END;
BEGIN C:=2/SQRT(3.14159265);
REPEAT WRITE('N,B0,B9,H?'); READLN(N,B0,B9,H);
K:=ROUND((B9-B0)/H+1.0); B:=B0; A:=0.0; S1:=0.0;
FOR I:= 1 TO K DO BEGIN TRAP(A,B,N,F,S); S1:=S1+S; A:=1.;
WRITELN(B,' ',C*S1); B:=B*H; END; UNTIL FALSE;END.
Формула СИМПСОНА
S= (b-a)/(6N)(f (x0) + f (x2N) + ∑i=12N-1 (3 + (-1) i-1) f (xi)) шаг:h=(b-a)/2N
Геометрическая интерпретация метода:
Пример программы на языке C#
Конец ознакомительного фрагмента.
Текст предоставлен ООО «ЛитРес».
Прочитайте эту книгу целиком, купив полную легальную версию на ЛитРес.
Безопасно оплатить книгу можно банковской картой Visa, MasterCard, Maestro, со счета мобильного телефона, с платежного терминала, в салоне МТС или Связной, через PayPal, WebMoney, Яндекс.Деньги, QIWI Кошелек, бонусными картами или другим удобным Вам способом.