model Pelotita
  parameter Real m=1;
  parameter Real k=1e6, b=30, g=9.81;
  Real y(start=1), v;
equation
  der(y) = v;
  der(v) = if y>0 then -g else -g-(k*y+b*v)/m;
end Pelotita;