t=0:0.5:20;

function dhdt=dhdt(h)
    dhdt=3-0.5*h;
endfunction

h=lsode("h",1.5,t);

plot(t,h)