0
scilab的哪些圖形函數可以通過圖中所示的圖獲得。創建一個點的三維圖
plots:-pointplots3d from Maple
在x,y,z座標是賦函數的結果。 謝謝。
我嘗試重現本文中的Scilab Solving Algebraic Equations by the Dragilev Method
scilab的哪些圖形函數可以通過圖中所示的圖獲得。創建一個點的三維圖
plots:-pointplots3d from Maple
在x,y,z座標是賦函數的結果。 謝謝。
我嘗試重現本文中的Scilab Solving Algebraic Equations by the Dragilev Method
我嘗試圖表,根據該腳本獲得的結果:
//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變量。我想獲得上面顯示的圖表
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()
嗨,歡迎來到[stackoverflow](https://stackoverflow.com/tour)。請通過添加更多細節來增強您的問題:x,y和z的大小是多少?什麼是頌讚?請提供[MWE](https://stackoverflow.com/help/mcve)。 – PTRK
請花時間寫出一個很好的問題。我們不會閱讀你的文章,以瞭解你在做什麼。這不是本網站的目的。看看[如何提出問題。](https://stackoverflow.com/help/how-to-ask) – PTRK
我嘗試繪製,根據此腳本獲得的結果: – HerClau