Вычислять определенные интегралы с заданной точностью эпсилон

[sefs]

Помогите с кодом, буду благодарен.
1. Методом Симпсона вычислить значение интеграла, разбивая отрезок интегрирования на 100 частей.
2. Вычислить значение интеграла с точностью 10 в -4 степени.
3. Сделать вывод, с какой точностью получен ответ в первом задании.

Снимок331.PNG

[kta87]

sefs, какую часть кода тебе помочь реализовать? Или может у тебя с математикой проблемы и ты не знаешь как решать определенные интегралы? а может проблемы и с кодом и с математикой? Суть в том, что нужно выкладывать попытки что то сделать и уже тогда будем разбираться что не получается!

[challengerr]

#include <stdio.h>
#include <math.h>

#define PI 3.14159265359

double f(double a) { return (sin(a) - cos(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 = PI;
d = 100.0;
fs = 0.0;
h = (b-a)/d;

printf("%f\n", 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 %f\n", j, s);

fs +=s;

}

printf("%f\n", fs);




}