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