%MATLAB-skript: mfd_tf_step.m
U=4; %Spranghoyden i u.
m=20; Kf=2; D=4;
%Definerer en LTI-modell med navn mfd vha. tf-funksjonen:
tellerpolynom=[1]; %Koeffisientene i tellerpolynomet:
nevnerpolynom=[m,D,Kf];  %Koeffisientene i tellerpolynomet:
mfd=tf(tellerpolynom,nevnerpolynom); %tf står for transfer function
tstopp=50;
[y1,t]=step(mfd,tstopp);%y1 er utgangsvektor.%t er tidsvektor.
%step-funksjonen antar at inngangen er enhetssprang.
y=y1*U; %Gir respons som svarer til inngangsspranghøyde U
%Plotting:
figure(1)
plot(t,y),ylabel('[m]'),xlabel('t [sek]'),title('Posisjon, y')
grid