2017-09-15 144 views
0

scilab的哪些圖形函數可以通過圖中所示的圖獲得。創建一個點的三維圖

plots:-pointplots3d from Maple

在x,y,z座標是賦函數的結果。 謝謝。

我嘗試重現本文中的Scilab Solving Algebraic Equations by the Dragilev Method

+0

嗨,歡迎來到[stackoverflow](https://stackoverflow.com/tour)。請通過添加更多細節來增強您的問題:x,y和z的大小是多少?什麼是頌讚?請提供[MWE](https://stackoverflow.com/help/mcve)。 – PTRK

+0

請花時間寫出一個很好的問題。我們不會閱讀你的文章,以瞭解你在做什麼。這不是本網站的目的。看看[如何提出問題。](https://stackoverflow.com/help/how-to-ask) – PTRK

+0

我嘗試繪製,根據此腳本獲得的結果: – HerClau

回答

0

我嘗試圖表,根據該腳本獲得的結果:

//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變量。我想獲得上面顯示的圖表

0
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()