def mk_line(n1,theta): x = np.linspace(0,n1-1,n1) y = np.copy(x)*0 count = 0 for t in x: y[count] = t*np.tan(theta*np.pi/180.) count +=1 return x,y-y[-1]/2