求两球交圆

function [yuanxin,r,xiangliang] = liangqiuxiangjiao(base1,base2,r1,r2)
%两球相交，求交圆 圆点yuanxin及半径r,以及法向量xiangliang
%   r1为圆点为base1的球的半径，r2为圆点为base2的球的半径
a=r1;b=r2;c=norm(base2-base1)
xiangliang=(base2-base1)/norm(base2-base1);
costheta=(a^2+c^2-b^2)/(2*a*c);%大臂与GHhand夹角
yuanxin=a*costheta*xiangliang+base1;%圆心
r=r1*(1-costheta^2)^0.5;%半径

end