Underdempning

> f:=t->Heaviside(t):

> masse_fjaer(2,5,2,0,0,15,15*f);