`ϵ` := 0.6104000000e23/(Pi*R^2); 1; R := 8000000; 1 

8000000 

K := 0.6672000000e-10; 1 

0.6672000000e-10 

phi := 2*Pi*K*`ϵ`*(sqrt(R^2+y^2)-y); 1 

.1272684000*(64000000000000+y^2)^(1/2)-.1272684000*y 

plot(phi, y = 0 .. 100000000); 1 

 

Plot 

 

 

plot(phi, y = 0 .. 1000000); 1 

 

Plot