from em_examples import sphereElectrostatic_example as electrostatic_sphere

sig0 = 10.**-3.
sig1 = 10.**-2.
sig2 = 1.310344828 * 10**-3.
R0    = 20.
R1 = 40.
E0   = 1.
n = 50
xr = np.linspace(-100, 100, n)
yr = xr.copy()
zr = np.r_[0]
XYZ = electrostatic_sphere.ndgrid(xr,yr,zr)
xstart = -100.
ystart = 50.
xend = 100.
yend = 50.
nb_dipole = 11
electrode_spacing = 20.
PlotOpt = 'Total'

#Plot Configuration, Potential and Data
electrostatic_sphere.two_configurations_comparison(XYZ,sig0,sig1,sig2,R0,R1,E0,xstart,ystart,xend,yend,nb_dipole,electrode_spacing,PlotOpt)

plt.tight_layout()
plt.show()