Hej :) mam problem z całka trygonometryczna.
Zadanie: Stosując wzór trapezów, Simpsona obliczyć całkę z (cosx)^2+1 na przedziale <0,2Pi>
Nie wiem jak zapisać tą całkę w programie.
//METODA TRAPEZOW
#include <iomanip>
#include <iostream>
#define M_PI 3.141 /pi/
using namespace std;
// Tutaj definiuje funkcję
double fun(double x);
{
return((cosx*cosx)+1);// nie wiem czy ten zapis jest poprawny
}
//Program główny
int main()
{
const int N = 2; //liczba punktów/trapezów podziałowych
double xp,xk,suma,h;
int i;
cout.precision(3); // 3 cyfry po przecinku
cout <<
"Podaj poczatek przedzialu calkowania\n\n"
"xp = ";
cin >> xp;
cout << "\nPodaj koniec przedzialu calkowania\n\n"
"xk = ";
cin >> xk;
cout << endl;
suma = 0;
h = (xk - xp) / N;
for(i = 1; i < N; i++)
{suma += fun(xp + i * h);}
suma = (suma + (fun(xp) + fun(xk)) / 2) * h;
cout << "Wartosc calki wynosi : " << setw(8) << suma
<< "\n\n";
system("PAUSE"); return 0;
}