Ayuda de Matlab, dibuja cuatro esferas en el sistema de coordenadas y encuentra las coordenadas centrales del área de intersección

clc

r=100;% radio

[x,y,z]=esfera(50);%coordenada de esfera unitaria

%bola central

x=r*x;y=r*y;z=r*z;

%bola uno

x1= x+50;y1=y;z1=z;

%Bola 2

x2=x+25;y2=y+25*sqrt(3);z2=z;

%Bola 3

x3=x+25;y3=y+30*sqrt(3);z3=z+50;

figura

figura

p>

espera

I1=(x-50).^2+y.^2+z.^2

I2=(x -25).^2+(y-25*sqrt(3)).^2+z.^2

I3=(x-25).^2+(y -25*sqrt(3)).^2+(z-50).^2

I=buscar(I1 & I2 & I3);

dispersión3(x(I),y(I),z(I),'.b')

I1=x1.^2+y1 .^2+z1.^2

I2=(x1-25) .^2+(y1-25*sqrt(3)).^2+z1.^2< r*r;

I3= (x1-25).^2+(y1-25*sqrt(3)).^2+(z1-50).^2

I=buscar(I1 & I2 & I3);

dispersión3(x1(I),y1(I),z1(I),'.c')

I1=x2.^2+y2.^ 2+z2.^2

I2=(x2- 50).^2+y2.^2+z2.^ 2

I3=(x2-25).^2+(y2-25*sqrt(3)).^2+(z2-50).

I=buscar(I1 & I2 & I3);

tamaño(I)

dispersión3(x2(I),y2(I),z2(I),'.m ')

I1=x3. ^2+y3.^2+z3.^2

I2=(x3-50).^2+y3. ^2+z3.^2

I3=(x3-25).^2+(y3-25*sqrt(3)).^2+z3.^2

I=buscar (I1 & I2 & I3);

dispersar3(x3(I),y3(I),z3(I),'.g')

%grid(x, y,z)

%grid(x1,y1,z1)

%grid(x2,y2,z2)

%grid(x3 ,y3,z3)

Los ejes son iguales

view(3)