#include #include /* set logscale x set logscale y plot "deriv.dat" pt 5 ps 2,x*x lc 0 */ double Eh(double h){ return fabs((sqrt(1+h)-sqrt(1-h))/(2*h)-1.0/2); } int main(){ double h=1; for(int i=0;i<20;i++){ printf("%g %g\n",h,Eh(h)); h/=2; } return 0; }