disp(' CELNI SOUKOLI SE SIKMYMI ZUBY ( uhly zadavany ve stupnich )')
alfa=input('Zadej uhel zaberu      alfa = ');
m=input('Zadej normalovy modul        m = ');
beta=input('Zadej uhel sklonu zubu beta = ');
z1=input('Zadej pocet zubu pastorku  z1 = ');
z2=input('Zadej pocet zubu kola      z2 = ');
i=z2/z1
D1=z1*m/cos(beta*pi/180)
D2=z2*m/cos(beta*pi/180)
at=(D1+D2)/2
aw=2*(fix((D1+D2)/4))+2
alfat=atan(tan(alfa*pi/180)/cos(beta*pi/180))
alfawt=acos(at*cos(alfat)/aw)
invalfat=tan(alfat)-alfat
invalfawt=tan(alfawt)-alfawt
SIGMAx=(z1+z2)*(invalfawt-invalfat)/(2*tan(alfa*pi/180))
x1=SIGMAx/(1+(1/i))
x2=x1/i
Da1=D1+2*m+2*x1*m
Da2=D2+2*m+2*x2*m
Df1=D1-2.5*m+2*x1*m
Df2=D2-2.5*m+2*x2*m