em hw 3

import numpy
import pylab

I = 2.0
A = 1.0

x_point = 1.0
y_point = 3.0
z_point = 1.0

R = 1.0

mu_not = 4*numpy.pi*10**-7

length = numpy.sqrt(x_point**2 + y_point**2 + z_point**2)

xhat = (x_point)/length
yhat = (y_point)/length
zhat = (z_point)/length

N = 100.

width = (2*numpy.pi)/N

dphi = numpy.arange(0,2*numpy.pi+0.001,width)
phi = 2*numpy.pi

intx = numpy.zeros(dphi.size)
inty = numpy.zeros(dphi.size)
intz = numpy.zeros(dphi.size)

for i in range(dphi.size):
    scriptr = ((R*xhat - A*numpy.cos(dphi[i])/length)**2 + (R*yhat - A*numpy.sin(dphi[i])/length)**2 + (R*zhat)**2)**(3./2.)

    intx[i] = width*((mu_not*I)/(4*numpy.pi))*(R*xhat - A*numpy.cos(dphi[i])/length)/scriptr

    inty[i] = width*((mu_not*I)/(4*numpy.pi))*(R*yhat - A*numpy.sin(dphi[i])/length)/scriptr

    intz[i] = width*((mu_not*I)/(4*numpy.pi))*((R*zhat)**2)/scriptr

x_sum = numpy.sum(intx)
y_sum = numpy.sum(inty)
z_sum = numpy.sum(intz)

print ('Coordinates of B field are '),x_sum,('T x^ ,'),y_sum,('T y^ ,'),z_sum,('T z^')