function section(t)
overlap = 0;
n = 1000; %生成椭球个数
ellipsoida = ones(1,n);
ellipsoidb = ones(1,n);
ellipsoidc = ones(1,n); %各个椭球的三半轴长矩阵
ellipsoidx0 = ones(1,n);
ellipsoidy0 = ones(1,n);
ellipsoidz0 = ones(1,n); %各个椭球的中心坐标矩阵
ellipsoidthetax = ones(1,n);
ellipsoidthetay = ones(1,n);
ellipsoidthetaz = ones(1,n); %各个椭球的旋转矩阵
r = ones(1,n); %各个椭球的有效半径
[x,y,z]=sphere(20); %先生成第一个椭球
abc = (rand(1,3) + 1) * 0.5; %三半轴长为0.5~1
ellipsoida(1) = abc(1);
ellipsoidb(1) = abc(2);
ellipsoidc(1) = abc(2);
[x,y,z] = f_scal(x,y,z,ellipsoida(1),ellipsoidb(1),ellipsoidc(1)); %拉伸
angle = pi * (2 * rand([1,3]) - 1);
ellipsoidthetax(1)