from numpy import *
%pylab inline
U = 3.6e5
eps = 8.854e-12
v = 10
d = 1e-6
def capacitor_area(u,v,d):
return 2*u*d/v/v/eps
print "You need %.2e m^2 plate area"%capacitor_area(U,v,d)
print "That is %.2e meters high!"%(d*capacitor_area(U,v,d))
mu = 1.26e-6
l = 2.
r = 1.
B = 10
print "Outward force: %.2e N"%(pi*l*r*B*B/mu)