T=[0:0.1:40];
[X, istate, MSG]=lsode("foscilw",[0 1],T)