(* ::Package:: *)

(*initialisation.m*)

(*Fix the distance x0 from the wire.*)
(*Choose the size of the interval*)
Deltaz=0.1;(* in dimensionless unit*)


(*The number of intervals NStep is determined by *)
NStep = IntegerPart[(2*L)/Deltaz];

z[0] = -L; z[NStep] = L; 
mu0=1.0;I0=1.0;mu0I=mu0*I0;
kini=0;klast=NStep;

xini=0.001;xlast=1.0;NStepx=20;
Deltax=(xlast-xini)/NStepx;
