/*f - jakaś tam funkcja, mniejsza
 a - początek przedziału
 b - koniec przedziału
 s - środek przedziału
 h - krok
 krok - to jest takie gówienko, które się przyda przy "przesuwaniu" przedziałów
 przedzialow - jak sama nazwa wskazuje :P
 C - obliczona całka
 C_stare - całka obliczona w poprzednim przebiegu pętli
 N - dokładność
*/

double a, b, h, s, tmp1, tmp2;
double C = 0.0, C_stare = 0.0;
int przedzialow = 0;
//tutaj użytkownik podaje przedział

do
{
    C_stare = 0.0;
    C = 0.0;
    przedzialow += 1;
    for(int i = 1; i <= przedzialow; i++)
{
krok = (a + b) / przedzialow;
tmp1 = a;
tmp2 = a + krok;
h = (tmp2 - tmp1) / 2;
s = (tmp1 + tmp2) / 2;
//tutaj obliczam wartości funkcji w tych punktach
C = C + //ten wzór simpsona na całkę
}
}while(C_stare - C > N)