#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)-2+sqrt(1+h))/(h*h)+1.0/4); } int main(){ double h=1; for(int i=0;i<13;i++){ printf("%g %g\n",h,Eh(h)); h/=2; } return 0; }