формула трапеций для вычисления интегралов

Ответ от Максим Петренко[новичек]
Даже не пытайся...

Ответ от Александр Сидорченко[гуру]
для первого примера//метод прямоугольников#include "stdio.h"#include "math.h"double f(double a) { return (a*a); } // вычисление значения функции в точкеvoid main(){double a; // нижний предел интегрированияdouble b; // верхний предел интегрированияdouble d; // количество отрезков разбиенияdouble h; /// длина отрезка разбиенияdouble a1; // начало отрезка разбиенияdouble b1; // конец отрезка разбиенияdouble s; // интегральная суммаdouble fs; // определенный интегралint j;a = 1.0;b = 4.0;d = 15.0;fs = 0.0;h = (b-a)/d;printf("%f
", h);for(j = 0; j < d; j++){a1 = a + (double)j * h;b1 = a1 + h;s = f ( (a1 + b1) / 2.0 ) * ( b1 - a1);printf("%d a+j*h:%f a1 +h:%f %f
", j, a1, b1, s);fs +=s;}printf("%f
", fs);} //метод трапеций#include "stdio.h"#include "math.h"double f(double a) { return (a*a); } // вычисление значения функции в точкеvoid main(){double a; // нижний предел интегрированияdouble b; // верхний предел интегрированияdouble d; // количество отрезков разбиенияdouble h; /// длина отрезка разбиенияdouble a1; // начало отрезка разбиенияdouble b1; // конец отрезка разбиенияdouble s; // интегральная суммаdouble fs; // определенный интегралint j;a = 1.0;b = 4.0;d = 15.0;fs = 0.0;h = (b-a)/d;printf("%f
", h);for(j = 0; j < d; j++){a1 = a + (double)j * h;b1 = a1 + h;s = ((f (a1) + f (b1)) / 2.0 ) * (b1 - a1);printf("%d a+j*h:%f a1 +h:%f %f
", j, a1, b1, s);fs +=s;}printf("%f
", fs);} // Метод Симпсона#include "stdio.h"#include "math.h"double f(double a) { return (a*a); } // вычисление значения функции в точкеvoid main(){double a; // нижний предел интегрированияdouble b; // верхний предел интегрированияdouble d; // количество отрезков разбиенияdouble h; /// длина отрезка разбиенияdouble a1; // начало отрезка разбиенияdouble b1; // конец отрезка разбиенияdouble s; // интегральная суммаdouble fs; // определенный интегралint j;a = 1.0;b = 4.0;d = 15.0;fs = 0.0;h = (b-a)/d;printf("%f
", h);for(j = 0; j < d; j++){a1 = a + (double)j * h;b1 = a1 + h;s = ( (b1-a1)/6.0 ) * ( f (a1) + 4.0 * f ( (a1 + b1) / 2.0 ) + f(b1) );printf("%d a+j*h:%f a1 +h:%f %f
", j, a1, b1, s);fs +=s;}printf("%f
", fs);} //для второго по Симпсонуdouble f(double a) { return (a); } // вычисление значения функции в точкеvoid main(){double a; // нижний предел интегрированияdouble b; // верхний предел интегрированияdouble d; // количество отрезков разбиенияdouble h; /// длина отрезка разбиенияdouble a1; // начало отрезка разбиенияdouble b1; // конец отрезка разбиенияdouble s; // интегральная суммаdouble fs; // определенный интегралint j;a = 4.0;b = 6.0;d = 10.0;fs = 0.0;h = (b-a)/d;printf("%f
", h);for(j = 0; j < d; j++){a1 = a + (double)j * h;b1 = a1 + h;s = ( (b1-a1)/6.0 ) * ( f (a1) + 4.0 * f ( (a1 + b1) / 2.0 ) + f(b1) );printf("%d a+j*h:%f a1 +h:%f %f
", j, a1, b1, s);fs +=s;}printf("%f
", fs);} //для третьего по Симпсонуdouble f(double a) { return (a*a*a); } // вычисление значения функции в точкеvoid main(){double a; // нижний предел интегрированияdouble b; // верхний предел интегрированияdouble d; // количество отрезков разбиенияdouble h; /// длина отрезка разбиенияdouble a1; // начало отрезка разбиенияdouble b1; // конец отрезка разбиенияdouble s; // интегральная суммаdouble fs; // определенный интегралint j;a = 0.0;b = 2.0;d = 10.0;fs = 0.0;h = (b-a)/d;printf("%f
", h);for(j = 0; j < d; j++){a1 = a + (double)j * h;b1 = a1 + h;s = ( (b1-a1)/6.0 ) * ( f (a1) + 4.0 * f ( (a1 + b1) / 2.0 ) + f(b1) );printf("%d a+j*h:%f a1 +h:%f %f
", j, a1, b1, s);fs +=s;}printf("%f
", fs);} //для четвертого по Симпсонуdouble f(double a) { return (exp(a)); } // вычисление значения

Ответ от Петя Росин[активный]
Сложновато

Ответ от Azamat kurman[новичек]
+

Ответ от Andrey Wayne[активный]
тут думать надо

Ответ от 3 ответа[гуру]