import os
os.sys.path.append("../pyspire")
from SpiresCoaxiales import *
import numpy
from matplotlib.pyplot import *
            

max=5.0
spire = SystemeSpires(-max,max,-max,max)
spire.ajouter(Spire(1.0,0.0,1.0))

            

points=[]
for k in range(7):
    points.append([k*0.15,0.0])
figure(figsize=(7,7))
spire.plot_lignes(points,'b')
axis([-max,max,-max,max])
xlabel('z')
ylabel('r')
grid()
            

max=20.0
spire.bornes(-max,max,-max,max)
points=[]
for k in range(10):
    points.append([k*0.1,0.0])
figure(figsize=(7,7))
spire.plot_lignes(points,'b')
axis([-max,max,-max,max])
xlabel('z')
ylabel('r')
grid()
            

z = numpy.arange(-5.0,5.0,0.01)
bz = spire.Bz_z(0.0,z)
figure()
plot(z,bz)
xlabel("z")
ylabel("Bz")
title("r=0")
grid()
            

x = numpy.arange(-3.0,3.0,0.01)
bz = spire.Bz_x(x,0.0)
figure()
plot(x,bz)
xlabel("r")
ylabel("Bz")
title("z=0")
axis([-3,3,-50,50])
grid()