Форум программистов, компьютерный форум, киберфорум
Наши страницы
Pascal ABC
Войти
Регистрация
Восстановить пароль
 
Рейтинг 4.60/5: Рейтинг темы: голосов - 5, средняя оценка - 4.60
DummyGirl
3 / 3 / 1
Регистрация: 12.10.2015
Сообщений: 202
1

Вычислить интеграл методом трапеции

03.03.2016, 07:42. Просмотров 954. Ответов 3
Метки нет (Все метки)

Найти
методом трапеций, разбивая

отрезок интегрирования на n = 22.

http://www.cyberforum.ru/cgi-bin/latex.cgi?<br />
\int_{0}^{\frac n2} \sin ^2x<br />

Так же найти с заданной точностью (выбрать e=0,001¸0,00001).
Есть наброски, но могу кинуть вечером.
Должны быть процедурки и менюшка. У меня там не получается, меньше и меньше числа.
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
Similar
Эксперт
41792 / 34177 / 6122
Регистрация: 12.04.2006
Сообщений: 57,940
03.03.2016, 07:42
Ответы с готовыми решениями:

Вычислить интеграл методом трапеций
помогите пожалуйста, буду очень благодарен, паскаль Вычислить интеграл методом...

Вычислить интеграл методом прямоугольников
Помогите пожалуйста, уже не первый день сижу: С погрешностью 0.001 вычислить...

Поиск интеграла методом трапеции
Найти интеграл cos(x) на отрезке (-1;2) методом трапеции. Заранее спасибо!

Вычисление интеграла методом трапеции
Помогите пожалуйста написать программу, решение функции методом трапеции=) ...

Реализовать программу вычисления методом трапеции
Реализовать программу вычисления методом трапеции интеграла интеграл от 0 до пи...

3
Puporev
Модератор
54582 / 42088 / 29061
Регистрация: 18.05.2008
Сообщений: 99,260
03.03.2016, 10:11 2
Лучший ответ Сообщение было отмечено DummyGirl как решение

Решение

В принципе решается так.
Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
uses crt;
function f(x:real): real;
begin
f:=sqr(sin(x));
end;
var a,b,e,h,s,x,y: real;
    i,n: integer;
begin
a:=0;
b:=pi/2;
n:=22;
s:=0;
h:=(b-a)/n;
for i:=1 to n-1 do
 begin
  x:=a+i*h;
  s:=s+f(x);
 end;
y:=h*((f(a)+f(b))/2+s);
writeln('При n=22 y=',y:0:5);
e:=0.001;
n:=round((b-a)/e);
s:=0;
for i:=1 to n-1 do
 begin
  x:=a+i*e;
  s:=s+f(x);
 end;
y:=e*((f(a)+f(b))/2+s);
writeln('При e=0.001 y=',y:0:5);
e:=0.00001;
n:=round((b-a)/e);
s:=0;
for i:=1 to n-1 do
 begin
  x:=a+i*e;
  s:=s+f(x);
 end;
y:=e*((f(a)+f(b))/2+s);
writeln('При e=0.00001 y=',y:0:5);
end.
0
DummyGirl
3 / 3 / 1
Регистрация: 12.10.2015
Сообщений: 202
03.03.2016, 16:32  [ТС] 3
Цитата Сообщение от Puporev Посмотреть сообщение
В принципе решается так.
Пыталась в процедуры засунуть, но чет не получается, где ошибка?
Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
//uses crt;
function f(x:real): real;
begin
f:=sqr(sin(x));
end;
var a,b,e,h,s,x,y: real;
    i,n: integer;
    procedure qwe(h:real);
    var i,j:integer;
begin
a:=0;
b:=pi/2;
n:=22;
s:=0;
h:=(b-a)/n;
for i:=1 to n-1 do
 begin
  x:=a+i*h;
  s:=s+f(x);
 end;
 end;
 
 procedure asd(Y:real);
 var i,j:integer;
 begin
y:=h*((f(a)+f(b))/2+s);
writeln('При n=22 y=',y:0:5);
e:=0.001;
n:=round((b-a)/e);
s:=0;
for i:=1 to n-1 do
 begin
  x:=a+i*e;
  s:=s+f(x);
 end;
y:=e*((f(a)+f(b))/2+s);
writeln('При e=0.001 y=',y:0:5);
e:=0.00001;
n:=round((b-a)/e);
s:=0;
for i:=1 to n-1 do
 begin
  x:=a+i*e;
  s:=s+f(x);
 end;
y:=e*((f(a)+f(b))/2+s);
writeln('При e=0.00001 y=',y:0:5);
end;
begin
qwe(H);
asd(y);
end.
Добавлено через 3 минуты
Вот так правильно будет?
Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
//uses crt;
function f(x:real): real;
begin
f:=sqr(sin(x));
end;
var a,b,e,h,s,x,y: real;
    i,n: integer;
 
procedure ndd(y:real);
var i,j:integer;
begin
a:=0;
b:=pi/2;
n:=22;
s:=0;
h:=(b-a)/n;
for i:=1 to n-1 do
 begin
  x:=a+i*h;
  s:=s+f(x);
 end;
y:=h*((f(a)+f(b))/2+s);
writeln('При n=22 y=',y:0:5);
end;
procedure eoo(y:real);
var i,j:integer;
begin
e:=0.001;
n:=round((b-a)/e);
s:=0;
for i:=1 to n-1 do
 begin
  x:=a+i*e;
  s:=s+f(x);
 end;
y:=e*((f(a)+f(b))/2+s);
writeln('При e=0.001 y=',y:0:5);
end;
 
procedure eooo(y:real);
var i,j:integer;
begin
e:=0.00001;
n:=round((b-a)/e);
s:=0;
for i:=1 to n-1 do
 begin
  x:=a+i*e;
  s:=s+f(x);
 end;
y:=e*((f(a)+f(b))/2+s);
writeln('При e=0.00001 y=',y:0:5);
end;
begin
ndd(y);
eoo(y);
eooo(y);
end.
0
Puporev
Модератор
54582 / 42088 / 29061
Регистрация: 18.05.2008
Сообщений: 99,260
03.03.2016, 16:58 4
Лучший ответ Сообщение было отмечено Памирыч как решение

Решение

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
const a=0;
      b=pi/2;
      n=22;
function f(x:real): real;
begin
f:=sqr(sin(x));
end;
procedure ndd;
var h,x,s,y:real;
    i:integer;
begin
s:=0;
h:=(b-a)/n;
for i:=1 to n-1 do
 begin
  x:=a+i*h;
  s:=s+f(x);
 end;
y:=h*((f(a)+f(b))/2+s);
writeln('При n=',n,' y=',y:0:5);
end;
procedure eoo(e:real);
var m,i:integer;
    x,s,y:real;
begin
m:=round((b-a)/e);
s:=0;
for i:=1 to m-1 do
 begin
  x:=a+i*e;
  s:=s+f(x);
 end;
y:=e*((f(a)+f(b))/2+s);
writeln('При e=',e:0:5,' y=',y:0:5);
end;
 
begin
ndd;
eoo(0.001);
eoo(0.00001);
end.
0
MoreAnswers
Эксперт
37091 / 29110 / 5898
Регистрация: 17.06.2006
Сообщений: 43,301
03.03.2016, 16:58

Вычислить площадь равнобедренной трапеции
Всем добрый вечер, срочно нужна помощь. Помогите кто может пожалуйста составить...

Написать функцию, которая считала бы интеграл на промежутке [a, b] методом трапеции
Добрый день. Знаю что раньше эта тема обсуждалась но все же обращусь ко всем за...

Методом трапеции вычислить площадь фигуры
Написать программу, которая методом трапеции вычислить площадь фигуры: тут...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
4
Ответ Создать тему
Опции темы

КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Рейтинг@Mail.ru