11 int main(int argc, char **argv) {
17 double delta=2*EPSILON;
22 printf("usage: %s file\n"argv[0]);
27 // for(l=0:l<MAX_L;l++)
30 for(r=0;r<MAX_R;r+=R_STEP) {
33 // fd=gp_init(argv[1]);
41 while(delta>EPSILON) {
46 delta=absolute_value(t/Jl);
49 printf("%f %f %f (i=%d)\n",r,Jl,Jl_1,i);