Source

wd_pys60 / ntbc.py

Full commit

a=1.7
b=0.18655
c=3.2
w1=4.087
from math import *

def ntbc(a,b,c,w1):
	print "  w         vd         v         ad        a\n"
	for w in range(0,360,10):
		p=w*pi/180
		q=b*sin(p)/(b*cos(p)+a)
		x=atan(q)
		n=cos(x)
		m=b*b+2*b*a*cos(p)+a*a
		w2=w1*(b*b+a*b*cos(p))/m
		y1=w2*c
		y2=n*c*w2
		t=a*b*w1*w1*sin(p)*(b*b-a*a)
		z0=t/(m*m)
		z1=z0*c
		z2=z0*c*n
		print "%4.4f    %4.4f     %4.4f    %4.4f    %4.4f"%(x*180/pi,y1,y2,z1,z2)

ntbc(a,b,c,w1)