r = 1;      % [m]
l = 3;      % [m]
m = 100;    % [kg]
theta = 20; % [deg ]nutation angle
% How long does it take the cylinder to precess through 180 ° if the
% spin rate is 2*pi/radians per minute
ws = 2*pi;                      % [rad/min]
C = m*r^2                       % [kg*m^2]
A = 1/2*m*r^2 + 1/12*m*l^2

C =

   100

A =

  125.0000

A > C, Direct (prograde) precession

wp = C/(A - C)*ws/cosd(theta)   % [rad/min]

wp =

   26.7457

The time for the spin axis to precess through an angle of 180 ° is

t = pi/wp*60                    % [sec]

t =

    7.0477

Published with MATLAB® 7.10