Matlab - Intersección de superficies

clear all

figure(gcf)

theta=0:pi/80:2*pi;

shi=0:pi/80:pi;

r=1;

[THETA, SHI]=meshgrid(theta,shi);

x=r*sin(SHI).*cos(THETA);

y=r*sin(SHI).*sin(THETA);

z=r*cos(SHI);

surf(x,y,z); shading flat

axis('equal')

hold on

x=-r:0.01:r;

y=x;

[X,Y]=meshgrid(x,y);

Z=0.5*ones(size(X));

surf(X,Y,Z,'edgecolor','r');

 clear all

 syms x y z

ec1=  x^2+y^2+z^2-1

ec2=z-0.5

 R=solve(ec1,ec2,x,y,z);

 x=double(R.x);

 y=double(R.y);

 z=double(R.z);

 plot3(x,y,z,'*')

 hold off

>> interseccion_esfera_con_plano