%   FY1001/TFY4145 Mekanisk fysikk november 2013
%   Oving 13
%   Faseplott for udempet harmonisk oscillator

% Tabell med tidspunkter fra 0 til 100, i alt 1001 verdier
% Sluttid 100 kan reduseres til under 2*pi 
% slik at mindre enn en periode vises
tid = linspace(0,6,1001);

figure(1);

subplot(1,1,1);
% Udempet oscillator, omega = 1, m = 1, slik at
% x(t)=sin(t) resulterer i p(t)=mv(t)=cos(t).
% A=1:
x11 = sin(tid);
p11 = cos(tid);
% A=2:
x12 = 0.5.*x11;
p12 = 0.5.*p11;

plot(x11,p11,x12,p12);
xlabel('x');
ylabel('p');
axis([-1 1 -1 1]);
title('Udempet');
