ShaftStresses:

% Computes shaft stresses.

 Type ShaftStresses for a demo.

--------------------------------------------------------------------------
   Form:
  [tau, alpha, j] = ShaftStresses( t, rI, rO, l, g )
--------------------------------------------------------------------------

   ------
   Inputs
   ------
   t         (1,1)   Torque (Nm)
   rO        (1,1)   Outer radius (m)
   rI        (1,1)   Inner radius (m)
   l        	(1,:)   Length (m)
   g         (1,1)   Modulus of rigidity (Pa)

   -------
   Outputs
   ------- 
   tau       (1,1)   Torsional stress (Pa)
   alpha     (1,1)   Twist angle
   j         (1,1)   Moment of inertia (m^4)

--------------------------------------------------------------------------
   Reference:  J. Dugundji, "16.20 Structural Mechanics Course Notes,"
               MIT, 1976.
--------------------------------------------------------------------------