clear
t0=0; %Starttid
t1=20; %Sluttid
x0=[0.5;0]; %Initialtilstand
global m D K R;
m=1;
D=0.2;
K=5;
R=0.5;
[t,x]=ode45(@mfd,[t0,t1],x0);
figure(1)
plot(t,x(:,1)) %Plotter x1=y som funksjon av t.
xlabel('t [sek]')
ylabel('x1(t) (posisjon)')