Bonjour,

j'ai essayé d'effectuer le code suivant pour représenter un ellipsoïde en utilisant une équation implicite, mais sans succès du coup ça serait pour avoir un avis

Code:
%d'après wikipedia
%(x^2/ a^2)+(y^2/ b^2)+(z^2/ c^2) = 1

N = length(x) 
N1= length(x) 

x = zeros(N1,N)
y = zeros(N1,N)
z = zeros(N1,N) 
for i = 1:N
for j = 1:N1
    x = 0 : 0.01 : 0.05
    y = -0.05 : 0.01 : 0.05
    z = -0.05 : 0.01 : 0.05
    a = 0.05
    b = 0.05
    c = 0.05
    [X,Y]  = meshgrid(x,y) 
    x(j,i) = sqrt(1-(y(j,i).^2/b^2)-(Z(j,i).^2/c.^2))*a.^2
    y(j,i) = sqrt(1-(X(j,i).^2/b^2)-(Z(j,i).^2/c.^2))*b.^2
    z(j,i) = sqrt(1-(X(j,i).^2/a^2)-(Y(j,i).^2/b.^2))*c.^2
 
end 
end

Code:
je précise que j'ai déjà fais la représentation paramétrique pour une demi ellipsoïde
N1 = 181
N = 180
x1 = zeros(N1,N);
y1 = zeros(N1,N);
z1 = zeros(N1,N);
for i = 1 : N
for j = 1 :N1
    u1 = 0 : 1: 180; %phi
    v1 = 0 : 1 : 180; %teta
    u = pi*u1/180;
    v = pi*u1/180;
    [U V] = meshgrid(u,v);  
    x1(j,i) = 10*cos(V(j,i)).*cos(U(j,i));
    y1(j,i)= 20*cos(V(j,i)).*sin(U(j,i)); 
    z1(j,i) = 20*sin(V(j,i)); 
end 
end
Bien cordialement