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
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
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
J'essaie de représenter graphiquement, les résultats obtenus selon ce script: – HerClau