function dx=masares(x,t)
  m=1;
  k=1;
  b=100;
  F=1;
  dx=[x(2);1/m*(-k*x(1)-b*x(2)+F)];
end