画长方体函数

function boxplot3(x0,y0,z0,Lx,Ly,Lz) 
    %(x0,y0,z0)是第一个顶点的位置; (Lx,Ly,Lz)是长方体的长宽高.
    x=[x0-Lx/2 x0-Lx/2 x0-Lx/2 x0-Lx/2 x0+Lx/2 x0+Lx/2 x0+Lx/2 x0+Lx/2];
    y=[y0-Ly/2 y0-Ly/2 y0+Ly/2 y0+Ly/2 y0-Ly/2 y0-Ly/2 y0+Ly/2 y0+Ly/2];
    z=[z0-Lz/2 z0+Lz/2 z0+Lz/2 z0-Lz/2 z0+Lz/2 z0-Lz/2 z0-Lz/2 z0+Lz/2];
    index=zeros(6,5);
    index(1,:)=[1 2 3 4 1];  %按一定顺序得到长方体角点的位置
    index(2,:)=[5 6 7 8 5];
    index(3,:)=[2 1 6 5 2];
    index(4,:)=[4 3 8 7 4];
    index(5,:)=[1 6 7 4 1];
    index(6,:)=[8 5 2 3 8];
    for k=1:6
        plot3(x(index(k,:)),y(index(k,:)),z(index(k,:)),'r');
        fill3(x(index(6,:)),y(index(6,:)),z(index(6,:)),'b');  %填充多边形函数
    hold on
    end
end