191
}
void
trap(){
double
integral=0, t;
for
(t=a+h; t<=b-h; t=t+h){
integral=integral+f(t);
}
integral=h*((f(a)+f(b))/2.0+integral);
cout<<
"integral trap="
<<integral<<endl;
}
void
simps(){
double
s1=0,s2=0,integral,t;
for
(t=a+h; t<=b-h; t=t+2.0*h)
s1=s1+f(t);
for
(t=a+2.0*h; t<=b-2.0*h; t=t+2.0*h)
s2=s2+f(t);
integral=h/3.0*(f(a)+f(b)+4.0*s1+2.0*s2);
cout<<
"integral simps="
<<integral<<endl;
}
void
centr(){
double
integral=0,tc;
for
(tc=a+h/2.0; tc<b; tc=tc+h){
integral=integral+f(tc);
}
integral=integral*h;
cout<<
"integral centr="
<<integral<<endl;
}
void
main(){
prav();
lev();
trap();
simps();
centr();
}
Результат выполнения программы представлен на рис. 6.11.
Рис. 6.11. Результат выполнения программы
Пример 6.8
Производство оборудования некоторого вида характеризуется
темпом роста его выпуска
yt
y
K
1
, где
y
- прирост выпуска этого
оборудования за промежуток времени
t
, а
y
– уровень его производства
за единицу времени на момент времени
t
. Найти общее количество
оборудования, произведенного за 10 лет, полагая, что
К
– известная
постоянная величина, единицей времени является год, а начальный
момент времени
t
=0.
