2017-09-15 13 views
0

Avec laquelle des fonctions graphiques de scilab pourrait être obtenu un graphique comme indiqué sur les figures.créer un graphique tridimensionnel des points

plots:-pointplots3d from Maple

x, y, z sont le résultat de la fonction ode. Merci.

J'essaie de reproduire cet article dans Scilab Solving Algebraic Equations by the Dragilev Method

+0

Salut, bienvenue à [stackoverflow] (https://stackoverflow.com/tour). S'il vous plaît améliorer vos questions en ajoutant plus de détails: quelles sont la taille de x, y et z? Quel est l'appel à ode? Veuillez fournir un [MWE] (https://stackoverflow.com/help/mcve). – PTRK

+0

Veuillez prendre le temps d'écrire une question bien formulée. Nous ne lirons pas votre article pour savoir ce que vous faites. Ce n'est pas le but de ce site Web. Jetez un oeil à [comment poser une question.] (Https://stackoverflow.com/help/how-to-ask) – PTRK

+0

J'essaie de représenter graphiquement, les résultats obtenus selon ce script: – HerClau

Répondre

0

J'essaie de représenter graphiquement les résultats obtenus selon ce script:

//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) 

Les résultats à projeter sont dans la variable LL. Et je voudrais obtenir le graphique ci-dessus

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