#include #include double f(double x){ return exp(-x*x); } double df(double x){ return -2*x*exp(-x*x); } double adf(double x,double h){ return (f(x+h)-f(x-h))/(2*h); } double a4df(double x,double h){ return (4*adf(x,h)-adf(x,2*h))/3; } int main(){ double x=1,h=1; for(int n=0;n<40;n++){ //double y=f(x); double dy=df(x); double ady=a4df(x,h); printf("%20.10e %20.10e %20.10e\n", h,ady,ady-dy); h/=2; } return 0; }