Respons på Heavisides enhetssprang
> u:=t->Heaviside(t):
> masse_fjaer(1,5,1,0,0,10,10*u);
Hastighet og posisjon i faseplanet