octave lsode function plotting

function xdot = f (x, t)
xdot(1) = 77.27 * (x(2) - x(1)*x(2) + x(1) - 8.375e-06*x(1)^2);
xdot(2) = (x(3) - x(1)*x(2) - x(2)) / 77.27;
xdot(3) = 0.161*(x(1) - x(3));
endfunction


x0 = [ 4; 1.1; 4 ];
t = linspace (0, 500, 1000);
y = lsode ("f", x0, t);

这个命令可以用octave的gnuplot模块看到y的输出函数数值
plot (y);

输出一个文件，这个文件不是xyz坐标的，则是矩阵排列的。

save ("ode_1.dat", "y");

使用openSCAD可以观看三维的效果，也可以用octave的plot3d （y)来实现。

//openSCAD to have a look
color("red")
surface(file = "ode_1.dat");