disp("Gear Speeds");


%Initial Paramaters (Gear radii in inches)
Nb=5.3256;
Na=0.66703;
Nc=0.55586;
Nd=1.9685;
in_speed= 4500; %[rpm]

%calculated prameters
g_ratio2= Nb/Na %ratio between gear A and B
g_ratio1= Nd/Nc %ratio betwen gear B and C

%RPM calculations for gear ratios
g_speed_A= in_speed
g_speed_B= g_speed_A/g_ratio2
g_speed_C= g_speed_B
g_speed_D= g_speed_C/g_ratio1

%Calculations to convert from RPM speed to Feet per minute
g_speed_A_fm= (g_speed_A*0.001885)*(2*Na*2.54) *54.681
g_speed_B_fm= (g_speed_B*0.001885)*(2*Nb*2.54) *54.681
g_speed_C_fm= (g_speed_C*0.001885)*(2*Nc*2.54) *54.681
g_speed_D_fm= (g_speed_D*0.001885)*(2*Nd*2.54) *54.681

%Gear quality equation specified in report.
Qv = "quality number"
%SOLVED WITH WOLFRAM
A.Qv = 1.92996;
C.Qv = -3.4;
%thus, gear quality number must be greater than max(A.Qv, B.Qv);

%Qv chosen value as 7
Qv=7
BC= 0.25*(12-Qv)^(2/3)
AC= 50 + 56*(1-BC)

%Dynamic Factor calculations given individual gear speeds and quality factor
Kv_a=((AC+sqrt(g_speed_A_fm))/AC)^BC
Kv_b=((AC+sqrt(g_speed_B_fm))/AC)^BC
Kv_c=((AC+sqrt(g_speed_C_fm))/AC)^BC
Kv_d=((AC+sqrt(g_speed_D_fm))/AC)^BC