function dx = f(x,t)
	dx(1) = -3.0*x
end

t=linspace(0,10,101);
x=lsode("f",1,t);
plot (t,x,";x(t);");
pause
