function y=funcao1 (~,x)

load parameters  m b L g 
y=[0 1;(-g/L) -(b/(m*L^2))]*x;

end