i++)
b=f(t+h,k); trapez = (a+b)*h/2; pole = pole + trapez; t=t+h; a=b;
cout << setprecision(7) << fixed; cout << " k="<<k<<" alfa="<<alfa<<"'"<<" całka E(k)=" << pole << endl;
return 0; }
Obliczenia numeryczne całek eliptycznych niezupełnych
(1) W celu obliczenia całki niezupełnej 2. rodzaju wystarczy w powyższym kodzie zmienić parametr t2 z wartości M_PI/2 (oznaczającego liczbę /2) na odpowiednią wartość parametru (linia 24 kodu).