(* ::Package:: *)

(*525.m*)
(*Eq. 5.25 is evaluated in this subroutine*)

tBy={};

For[i=0,i<=NStepx,i++,
	x = xini + i*Deltax;
	By=0;Bydum=0;
    	For[k=kini,k<=klast,k++,
		Bydum=(mu0I/(4*Pi))*(x*Deltaz)/( (k*Deltaz)^2 + x^2)^(3/2);
		By=By+Bydum;
		];(*end for k*)

tBy=Append[tBy,{x,By}];
];(*end for i*)
(*Print[tBy];*)
(*Print[ListPlot[tBy,Joined->True,AxesLabel->{"x","By"},PlotMarkers->Automatic]];*)
