figure (1)
clf;
clear;

robot = tf(1, [1 8 25]);		%robot with m=1,b=8,k=25

Kp = 25;			% defines the proportional controller gain

t=0:0.01:2; %sample at 100Hz for 2 seconds

perturbation = 0;       %no perturbation

open_loop = Kp*robot+perturbation;

figure(1)
step(open_loop,t)		% give a step input with magnitude of 1 meter and plot the output
title('open loop')

