#include #include main(){ int N, i, j, k; double h,t,c; double yps,minyps,ypsr,minminyps; double x,r,hl; double xi,xil,yj; /* printf("N: "); scanf("%d",&N);*/ N=500; r = 0.5; yps = 1.0; ypsr=1.0e-14; minminyps=1.0; for(k = 2 ;k<=N; k++){ k = k*2; h = 1.0/k; minyps = 1.0; for (i = 1 ; i < k/2 ; i++){ xi = r-i*h; xil = r-(i-1)*h; for (j = 1 ; jr){ yps=(sqrt(r*r-yj*yj)-xi)/h; if( yps < minyps ){ minyps = yps; if(minyps < minminyps ){ minminyps = minyps; } } } } } if((k/2)%5==0){ /*printf("N=%3d εの最小値=%g ", k/2 , minyps);*/ printf("%d %g\n",k/2,minyps); } k=k/2; } /*printf ("\n 最小最小ε=%g \n",minminyps);*/ }