disp(' NAVRH CELNIHO SOUKOLI S PRIMYMI A SIKMYMI ZUBY')
beta=input('Zadej uhel sklonu zubu      beta [ ° ] = ');
psi=input('Zadej soucinitel sirky ozubeni     psi = ');
c=input('Zadej materialovou konstantu c [ MPa ] = ');
i=input('Zadej prevodovy pomer                i = ');
P=input('Zadej prenaseny vykon         P [ kW ] = ');
n1=input('Zadej otacky pastorku     n1 [ 1/min ] = ');
z1=input('Zadej pocet zubu pastorku           z1 = ');
z2=round(i*z1)
i_skutecny=z2/z1
rozdil_v_procentech=100*abs(i-i_skutecny)/i
Mk1=1000*P/(pi*n1/30)
if beta == 0
m=8.6025*(Mk1/(c*psi*z1))^(1/3)
else
m=7.515*(Mk1*cos(beta*pi/180)/(c*psi*z1))^(1/3)
end
m_norm=input('Zadej normalizovany modul m_norm = ');
D1=z1*m_norm/cos(beta*pi/180)
v=0.001*pi*D1*n1/60
b=psi*m_norm