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

x0= [904.8395;
0.5641;
166.9104;
0.1756;
199.8356;
0.0200;
0.1545;
99.9363;
0.0837;
0.1545;
95.2973;
9.0081;
38.4593;
0.0351;
39.9671;
0.0040;
0.0309;
19.9873;
0.0167;
0.0309;
0.1047;
0.1990;
2.0000;
0.1531;
1.9984;
0.0027;
2.0000;
1.9984;
2.0000;
0.0027;
1.9987;
1.9987;
0.0027;
55705.0;
0.0617;
0.0008;
0.0478;
0.0478;
0;
1003.0;
0;
1.8953;
1.8010;
1.8469;
1.9973;
1.9973;
1.9973;
0.0016;
0.0016;
0.0000;
0.0000;
0.0013;
0.0013;
0;
0];

x = lsode (dxdt,x0,t);

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

endfunction