%resonance02.m
%Plot Amplitude

w=[0:0.01:2];

fo=1;
k=1;
m=1;

a1=fo./(5*k/-10/9*m*w.^2);

figure(1)
plot(w,a1,'r-')
axis([0 2 -20 20])
xlabel('\omega')
ylabel('a')
title('Amplitude vs. Frequency')
hold on
plot([1 1],[-20,20],'k:')
plot([0 2],[0,0],'k-')
hold off
