Метод трапеций
Наименьший допускаемый диаметр базового отверстия дискового резца из условия виброустойчивости
Окончание табл. 15
Контрольная работа по предмету «Вычислительная математика»
Выполнил: Студент группы 2-70а Кишалов. Д.В Шифр 911105
Иваново - 2013
Задание №1 Численное интегрирование
#include <stdio.h> #include <math.h> const double Pi = 3.141593;
double func (double x) { double w; w =pow(2.0,3.0*x); // Указать подинтегральную функцию return w; }
int main () { FILE * f; fopen_s (&f, "result_v1.dat", "w"); double a, b, h = 0.01; a =0; // Указать значение a b =1.0; // Указать значение b int n = (int) ((b-a)/h); // Метод прямоугольников в среднем double x = a + h, integr = 0.0; for (int i=0; i<n; i++) { integr += h * func (x-0.5*h); x += h; } printf ("Metod prjam - integr = %2.8f\n",integr); fprintf (f, "Метод прямоугольников - интеграл = %2.8f\n",integr); // Метод трапеций integr = 0.0; x = a + h; for (int i=0; i<n-1; i++){ integr += h * func(x); x += h; } integr += 0.5 * h * (func(a) + func(b)); printf ("Metod trap - integr = %2.8f\n",integr); fprintf (f, "Метод трапеций - интеграл = %2.8f\n",integr); // Метод Симпсона integr = 0.0; x = a + h; if (n & 1) n++; double sum0 = 0.0, sum1 = 0.0; for (int i=1; i<n; i++) { if (i & 1) sum1 += func(x); else sum0 += func(x); x += h; } integr = h*(func(a)+4.0*sum1+2.0*sum0+func(b))/3.0; printf ("Metod Simpson - integr = %2.8f\n",integr); fprintf (f, "Метод Симпсона - интеграл = %2.8f\n",integr); fclose (f); return 0; }
Результаты вычислений:
Метод прямоугольников - интеграл = 3.36622778 Метод трапеций - интеграл = 3.36640973 Метод Симпсона - интеграл = 3.36628843
Блок схема Метод прямоугольников всреднем
Метод трапеций
|