Нахождение интегралов в среде Pascal
Нахождение интегралов в среде Pascal
Нахождение интегралов в среде Pascal Текст задания: найти интеграл функции sin (x+10)+x4=0 используя два метода: метод прямоугольников и метод трапеций; пользователь вводит значения и выбирает метод решения. Текст программы: program LW_5; uses Crt; var dx, x1, x2, e, i, c, s:real; function Fx (x:real):real; begin Fx:=sin (x+10)+x*x*x*x; {function} end; procedure CountViaBar; var xx1, xx2:real; c:longint; begin writeln ('-'); writeln ('->MEthod PRYAMOUGOL''NIKOVVVV'); writeln ('vsego interaciy:', round (abs(x2-x1)/e)); i:=0; for c:=1 to round (abs(x2-x1)/e) do begin write ('integration ', c, chr(13)); xx1:=Fx (x1+c*e); xx2:=Fx (x1+c*e+e); i:=i+abs (xx1+xx2)/2*e; end; writeln ('-'); writeln ('INTEGRAL=', i); end; procedure CountViaTrap; var xx1, xx2, xx3:real; c:longint; begin writeln ('-'); writeln ('->MEthod TRAPECIY'); writeln ('vsego interaciy:', round (abs(x2-x1)/e)); i:=0; for c:=1 to round (abs(x2-x1)/e) do begin write ('Integration ', c, chr(13)); xx1:=Fx (x1+c*e); xx2:=Fx (x1+c*e+e); if xx2>xx1 then xx3:=xx1 else xx3:=xx2; i:=i+abs (xx2-xx1)*e+abs(xx3)*e; end; writeln ('-'); writeln ('INTEGRAL=', i); end; begin clrscr; writeln ('vvedite ishodniye zna4eniya:'); write ('NA4AL''NOYE ZNE4ENIYE x (x1)='); Readln(x1); write ('KONE4NOYE ZNA4ENIYE x (x2)='); Readln(x2); write ('TO4NOST'' VbI4ISLENIY (e)='); Readln(e); writeln ('-======vvedite metod====== - :'); writeln ('DLYA METODA PRYAMOUGOL''NIKOV VVEDITE 1, DLYA METODA TRAPECIY VVEDITE 2'); //writeln ('DLYA METODA TRAPECIY VVEDITE 2'); readln (s); if s=1 then begin CountViaBar; end; if s=2 then begin CountViaTrap; end; readln; end.
|