double integrand(double x, double param){ return x * sin(x * param); } double primal(double x, double param){ return sin(x * param) / pow(param, 2) - x * cos(x * param) / param; }