%FY1001/TFY4145 Mekanisk fysikk høsten 2011.
%Oving 4, Oppgave 3c.
clear all;
%N = antall punkter langs omega-aksen
N = 200;
%g = tyngdens akselerasjon
g = 9.81;
%R = ringens radius
R = 0.10;
%tabell
omega_max = 40;
theta_max = pi/2;
omega = linspace(0.0,omega_max,N);

for i = 1:N
    nevner(i) = R*omega(i).^2;
    if nevner(i) < g
        theta(i) = 0.0;
    else
        theta(i) = acos(g/nevner(i));
    end;
end;
% Plotter omega(theta)
plot(omega,theta*180/pi,'linewidth',2);
axis([0.0 omega_max 0.0 theta_max*180/pi]);
% Aksetekster:
xlabel('Omega (1/s)');
ylabel('Theta (grader)');