function io = integrate_octave(start,step,steps,x0,dxdt)
t = linspace(start,start+steps*step,steps+1);

x = stoch_sim (dxdt,x0,t);

io.t = t;
io.x = x;

endfunction
