%   FY1001/TFY4145 Mekanisk fysikk november 2012
%   Oving 13
%   Faseplott for dempet harmonisk oscillator

% Tabell med tidspunkter fra 0 til 100, i alt 1001 verdier
tid = linspace(0,100,1001);

figure(1);

subplot(1,1,1);
% Udempet oscillator, A = 1 og 0.5, omega = 1, m = 1
% slik at x(t)=sin(t) resulterer i p(t)=mv(t)=cos(t):
x11 = sin(tid);
p11 = cos(tid);
x12 = 0.5.*x11;
p12 = 0.5.*p11;
plot(x11,p11,x12,p12);
xlabel('x');
ylabel('p');
axis([-1 1 -1 1]);
title('Udempet');
