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

t=[0:0.1:10];      # 横軸 t の作成．0 から 10 を0.1 おきに． t は横ベクトル． 
x=lsode("f",1,t);  # 1 は初期値 x(0)． こう書けば解いてくれる．
z = [t', x];       # t' は t の転置
save("-ascii","ode001.dat","z" );
plot (t,x,";x(t);");
system("gnuplot gp001 > kadai001.eps");
pause
