%
%  file for plotting output and getting average amplitudes
%
index = floor(length(time)/2);
plot(time(index:end),theta1(index:end)); grid;
disp('Average Amplitude is ')
A = 0.5*(max(theta1(index:end))-min(theta1(index:end)))