Метод Симпсона (численное интегрирование)
Primary tabs
#include stdio.h #include conio.h #include iostream #include cmath using namespace std; double f_1(double x){ double S = log(sin(x) + x); return S; } double Simp(double a, double b, double h){ double sum = f_1(a) + f_1(b); for (double x = a + h; x < b; x += 2 * h) { if (x + h >= b) sum += 2 * f_1(x); else sum += 2 * f_1(x) + 4 * f_1(x + h); } return sum * h / 3.0; } int main(int argc, char* argv[]){ cout << Simp(2, 4, 0.001) << endl; getch (); return 0; }
- Log in to post comments
- 4990 reads
vedro-compota
Tue, 04/22/2014 - 11:55
Permalink
no comments in code? =))
no comments in code? =))
_____________
матфак вгу и остальная классика =)