#include double em(double x) { return(exp(-0.1*x)*sin((1.0 - 0.2*x)*x + 2.0)); } double emp(double x) { double w; w = 1.0 - 0.2*x; return(exp(-0.1*x)*(-0.1*sin(w*x + 2.0) + cos(w*x + 2.0)*(-0.2*x + w))); }