%MATLAB-skript: mfd_ode45.m
clear
t0=0; %Starttid
t1=50; %Sluttid
x0=[0;0]; %Initialtilstand
global m Kf D;
m=20;Kf=2;D=4;
[t,x]=ode45(@mfd,[t0,t1],x0);
figure(1)
plot(t,x(:,1)) %Plotter x1=y som funksjon av t.
grid
title('Simulering av masse-fjær-demper-system')
xlabel('t [sek]')
ylabel('y = x1(t) (posisjon)')