Question

可以获得scilab的哪些图形函数，如图所示。

x，y，z坐标是ode函数的结果。谢谢。

我尝试在Scilab中重现这篇文章 Solving Algebraic Equations by the Dragilev Method

Answer 1

我尝试绘制图表，根据此脚本获得的结果：

//Draghilev's Method
clear;
function Sys=Q(x,z)
Sys(1)=(x(1)-2)^2+(x(2)+2)^2+z^2-9;
Sys(2)=x(1)^6+x(2)^6+z^6-12;
endfunction
z=0;
[j,v,info]=fsolve([2;-0.1],list(Q,z))
disp(j,v,info)

//
N=100;
smin=0.0;
smax=0.046;
h=0.001;
x01=3.9691163496*10^-12;
x02=0.353580783;
x03=-1.530442808;
ics=[x01; x02; x03];

disp(["x(3)^2+(x(2)-2)^2+(x(1)-2)^2-9";"x(3)^6+x(2)^6+x(1)^6-12"],"systema")

disp(["2*(x(1)-2),2*(x(2)-2),2*x(3)";"6*x(1)^5,6*x(2)^5,6*x(3)^5"],"jacobian(Vm,[x(1),x(2),x(3)])"); //

//D1=12*x(2)^5*x(3)-12*(x(2)-2)*x(3)^5;
//D2=12*(x(1)-2)*x(3)^5-12*x(1)^5*x(3);
//DD=12*x(1)^5*(x(2)-2)-12*(x(1)-2)*x(2)^5;
function dydt=odes(t,x)
//dydt=zeros(x);
dydt(1)=-12*x(2)^5*x(3)-12*(x(2)-2)*x(3)^5;
dydt(2)=12*(x(1)-2)*x(3)^5-12*x(1)^5*x(3);
dydt(3)=12*x(1)^5*(x(2)-2)-12*(x(1)-2)*x(2)^5;
endfunction
step=0.0000005;
t=[smin:step:smax];
t0=0
atol=h/100000;
LL= ode( ics,t0, t,atol,odes)

要绘制的结果在LL变量中。我想得到上面显示的图表

创建一个三维的积分图

2 个答案: