disp(' KONTROLA CELNIHO SOUKOLI S PRIMYMI A SIKMYMI ZUBY')
m=input('Zadej normalovy modul          m [ mm ] = ');
beta=input('Zadej uhel sklonu zubu       beta [ ° ] = ');
psi=input('Zadej soucinitel sirky ozubeni      psi = ');
mi=input('Zadej soucinitel presnosti ozubeni   mi = ');
Lh=input('Zadej dobu behu              Lh [ hod ] = ');
P=input('Zadej prenaseny vykon          P [ kW ] = ');
n1=input('Zadej otacky pastorku      n1 [ 1/min ] = ');
z1=input('Zadej pocet zubu pastorku            z1 = ');
z2=input('Zadej pocet zubu kola                z2 = ');
b=psi*m
D1=z1*m/cos(beta*pi/180)
v=0.001*pi*D1*n1/60
n2=n1/(z2/z1)
Mk1=1000*P/(pi*n1/30)
zn1=z1/(cos(beta*pi/180))^3
u=(m/10)^0.2
sigmaDo1=input('Zadej dovolene napeti pro ohyb pro pastorek     sigmaDo1 [ MPa ] = ');
sigmaDd1=input('Zadej dovolene napeti pro otlaceni pro pastorek sigmaDd1 [ MPa ] = ');
sigmaDo2=input('Zadej dovolene napeti pro ohyb pro kolo         sigmaDo2 [ MPa ] = ');
sigmaDd2=input('Zadej dovolene napeti pro otlaceni pro kolo     sigmaDd2 [ MPa ] = ');
ro1=input('Zadej rychlostni soucinitel pro ohyb pro pastorek     ro1 = ');
ro2=input('Zadej rychlostni soucinitel pro ohyb pro kolo         ro2 = ');
yo1=input('Zadej tvarovy soucinitel pro ohyb pro pastorek        yo1 = ');
yo2=input('Zadej tvarovy soucinitel pro ohyb pro kolo            yo2 = ');
rd1=input('Zadej rychlostni soucinitel pro otlaceni pro pastorek rd1 = ');
rd2=input('Zadej rychlostni soucinitel pro otlaceni pro kolo     rd2 = ');
yd=input('Zadej tvarovy soucinitel pro otlaceni                  yd = ');
co1=sigmaDo1*ro1/yo1
co2=sigmaDo2*ro2/yo2
cd1=sigmaDd1*rd1/(yd*u)
cd2=sigmaDd2*rd2/(yd*u)
cmin=min([co1 co2 cd1 cd2])
Fd=pi*cmin*b*m*mi
Fo=2000*Mk1/D1