点の3次元プロットを作成する

Question

図に示すように、scilabのグラフィカル関数のどれをグラフで得ることができるかを指定します。

plots:-pointplots3d from Maple

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変数です。そして、私は上記のグラフを得たいと思います。

Answer 2

clf; 
drawlater() 
param3d(x,y,z); 
set(gce(),"mark_mode","on","mark_style",2,"mark_size_unit","point","mark_size",4,"line_mode","off") 
drawnow()