disp（‘*********直齒輪的設計******‘）;
disp（‘********初步計算********‘）;
disp（‘輸入已知條件‘）;
PI=pi/180;%角度變換成弧度的系數；
alpha=20*PI;%齒輪壓力角；
P1=input（‘齒輪所傳遞的功率P1=‘）;
n1=input（‘小齒輪的轉速n1=‘）;
n2=input（‘大齒輪的轉速n2=‘）;
t=input（‘輸入壽命t=‘）;
CLCL=input（‘選擇大小齒輪的材料‘）;
disp（‘查圖2.1-16至2.1-27取得相關的參數‘）;
sigma_Hlim1=input（‘小齒輪的試驗齒輪的接觸疲勞極限sigma_Hlim1=‘）;
sigma_Hlim2=input（‘大齒輪的試驗齒輪的接觸疲勞極限sigma_Hlim2=‘）;
sigma_Flim1=input（‘小齒輪的試驗齒輪的彎曲疲勞極限sigma_Flim1=‘）;
sigma_Flim2=input（‘大齒輪的試驗齒輪的彎曲疲勞極限sigma_Flim2=‘）;
T1=9550*P1/n1;%計算小齒輪的轉矩；
K=input（‘初選載荷系數K=‘）;
u=input（‘輸入齒數比u=I;u=‘）;
if?sigma_Hlim1>=sigma_Hlim2
????????sigma_Hlim=sigma_Hlim2;
????else
????????sigma_Hlim=sigma_Hlim1;
end
????SHmin=input（‘查表達式2.1-45初選接觸應力安全系數SHmin=‘）;
????sigma_HP=sigma_Hlim/SHmin;%sigma_HP為許用接觸應力；
?????if?sigma_Flim1>=sigma_Flim2
????????sigma_Flim=sigma_Flim2;o=2;
????else
????????sigma_Flim=sigma_Flim1;o=1;
?????end
?????SFmin=input（‘查表達式2.1-45初選彎曲應力安全系數SFmin=‘）;
????sigma_FP=sigma_Flim/SFmin;%sigma_FP為彎曲應力；
????SJGS=input（‘輸入所使用的設計公式：‘JCYL’（接觸應力），‘WQYL’（彎曲應力）‘）;
if?SJGS==‘JCYL‘
????Ja=input（‘齒面接觸強度計算系數Ja=‘）;
????phi_a=input（‘初選齒寬系數（參考148頁）；phi_a=‘）;
????if?sigma_Hlim1>=sigma_Hlim2
?????o=2;
????else
?????o=1;
????end
????sigma_HP=sigma_Hlim/SHmin;%sigma_HP為許用接觸應力；
????disp（‘***按齒面接觸強度計算****‘）;
????LH=input（‘輸入嚙合情況：（外嚙合）‘+’；（內嚙合）‘-’；LH=‘）;
????if?LH==‘+‘
????????amin=Ja*（u+1）*（K*T1/（phi_a*u*sigma_HP^2））^（1/3）;%amin為最小中心距；
????????a=input（‘選取中心距a>=amin;a=‘）;
????elseif?LH==‘-‘
????????amin=Ja*（u-1）*（K*T1/（phi_a*u*sigma_HP^2））^（1/3）;%amin為最小中心距；
????????a=input（‘選取中心距a>=amin;a=‘）;????????
????end
????disp（‘根據經驗公式：平穩載荷m=（.007~.01）a;中等沖擊m=（.01~.015）a;‘）;
????????disp（‘較大沖擊m=（.015~.02）a;‘）;
????????m=input（‘確定齒輪模數（查表2.1-487頁）m=‘）;
elseif?SJGS==‘WQYL‘
????phi_m=input（‘初選齒寬系數（參數141頁）；phi_m=‘）;
????Z1=input（‘初選小齒輪齒數（查表2.1-39147頁）；Z1=‘）;
????YFS=input（‘初選大小齒輪復合齒形系數（查2.1-342.1-35）;YFS=‘）;
????if?sigma_Flim1>=sigma_Flim2
?????o=2;
????else
???????o=1;
????end
????SFmin=input（‘查表達式2.1-45初選接觸應力安全系數SFmin=‘）;
????sigma_FP=sigma_lim/SFmin;%sigma_FP為彎曲應力；
????disp（‘****按齒根彎曲強度計算****‘）;
????mmin=12.5*（K*T1*YFS/（phi_m*Z1^2*sigma_FP））^（1/3）;%mmin為最小模數；
????m=input（‘選取模數m（參考2.1-487頁）；m=‘）;
????disp（‘根據經驗公式：平穩載荷m=（.007~.01）a;中等沖擊m=（.01~.015）a;‘）;
????disp（‘較大沖擊m=（.015~.02）a;‘）;
????a=input（‘確定中心距a=‘）;
end
Z1=2*a/（m*（u+1））;
fprintf（1‘小齒輪齒數Z1約等于%3.3f\n‘Z1）;
Z1=input（‘選取小齒輪齒數Z1=‘）;
Z2=round（u*Z1）;%Z2為大齒輪的齒數；
d1=m*Z1;%d1為小齒輪分度圓直徑；
if?SJGS==‘JCYL‘
????b=phi_a*a;%b?為齒輪寬度；
elseif?SJGS==‘WQYL‘
????b=phi_m*m;
end
disp（‘****校核齒面接觸疲勞強度****‘）;
disp（‘***修正載荷系數計算***‘）;
KA=input（‘使用系數（查表2.1-46161頁）;KA=‘）;
K1=input（‘計算Kv中的K1（查表2.1-48161頁）；K1=‘）;
K2=input（‘計算Kv中的K2（查表2.1-48161頁）；K2=‘）;
Ft=2000*9550*P1/（d1*n1）;%Ft為分度圓切向力；
Kv=1+（K1*b/（KA*Ft）+K2）*Z1*pi*n1*d1/1.2e7*（u^2/（u^2+1））^.5;%Kv為動載系數；
s=input（‘s為齒輪與軸中心（齒輪在支點中間）或軸支點（齒輪在支點外）的距離；s=‘）;
l=input（‘l為軸支點間的距離；l=‘）;
if?s/l<.1
????disp（‘齒輪的結構布局為對稱支承‘）;
elseif?s/l>.1&&s/l<.3
????disp（‘齒輪的結構布局為非對稱支承‘）;
elseif?s