float eq_sphere(float q , float x , float y , float a);
float eq_dotXY(float q , float x , float y);
float get_disk(float x , float y);
float eq_sphere(float q , float x , float y , float a){
float r;
r = get_disk(x , y);
if (r < a){
r = a;
}
return I4PE*q/fabs(r);
}
float ep_dotXY(float q , float x , float y){
float r;
r=get_disk(x , y);
return I4PE*q/fabs(r);
}
float get_disk(float x , float y){
return sqrt(x*x+y*y);
}
