DRAW_PRO.m

function[Xr,Yr,Zr,XC1,YC1,ZC1,XC2,YC2,ZC2,XC12,YC12,ZC12,XC22,YC22,ZC22,XYZT,XYZC1,XYZC2,XA,YA,ZA,XC1i,YC1i,ZC1i,XC2i,YC2i,ZC2i,n,col]=DRAW_PRO(XYZt,XYZc1,XYZc2,MTC1,MTC2,D,r,F)
fac=1; 
 ZZ0=r(1:110)*fac;

  XYZ=[XYZt(:,1) XYZt(:,2) XYZt(:,3)+ZZ0-40];
%GRAF beam
for i=2:length(XYZ(:,1))
    p1=(XYZ(i,:)-XYZ(i-1,:));
    angle(i-1)=atan(p1(1,3)/p1(1,2));
end
Angle(1)=angle(1);
for i=2:length(angle)
    Angle(i)=(angle(i-1)+angle(i))/2;
end  
Angle(i+1)=angle(end);
Xt=[225 225 -225 -225 -175 -175 175 175 225];
Yt=zeros(1,9);
Zt=[-35 25 25 -35 -35 0 0 -35 -35];
for i=1:length(Angle)
    th=Angle(i);
    MT=[1 0 0;0 cos(th) -sin(th);0 sin(th) cos(th)];
    for j=1:length(Xt)
        VR=MT*[Xt(j);Yt(j);Zt(j)];
        Xr(i,j)=VR(1)+XYZ(i,1);Yr(i,j)=VR(2)+XYZ(i,2);Zr(i,j)=VR(3)+XYZ(i,3);
    end
end
%GRAF COLUMN1
Xc1=[55 55 -55 -55 55];
Yc1=[-55 55 55 -55 -55];
Zc=[0 0 0 0 0];
 YY=[0;r(111:143)]*fac;
 ZZ=[0;r(177:209)]*fac;
  XYZ2=[XYZc1(:,1) XYZc1(:,2)+YY XYZc1(:,3)+ZZ];
for i=2:(length(XYZ2(:,1)))
    p1=(XYZ2(i,:)-XYZ2(i-1,:));
    angle2(i-1)=-atan(p1(1,2)/p1(1,3));
end
Angle2(1)=0;
for i=2:length(angle2)
    Angle2(i)=(angle2(i-1)+angle2(i))/2;
end  
Angle2(i+1)=angle2(end);
facx=[ones(1,8),linspace(1,0.8182,17),ones(1,9)*0.8182];
ampx=[linspace(180,290,6),290,linspace(290,45,18),ones(1,9)*45];
for i=1:length(Angle2)
    th=Angle2(i);
    MT2=[1 0 0;0 cos(th) -sin(th);0 sin(th) cos(th)];
    for j=1:5
        VR2=MT2*[Xc1(j)*facx(i);Yc1(j);Zc(j)];
       XC1(i,j)=VR2(1)+XYZ2(i,1)+ampx(i);
       YC1(i,j)=VR2(2)+XYZ2(i,2);
       ZC1(i,j)=VR2(3)+XYZ2(i,3);
       XC2(i,j)=VR2(1)+XYZ2(i,1)-ampx(i);
       YC2(i,j)=VR2(2)+XYZ2(i,2);
       ZC2(i,j)=VR2(3)+XYZ2(i,3);
   end
end
XC1i=[[XC1(6,1:2) XC1(8,2) XC1(8,1) XC1(6,1)]-110;[XC2(6,4) XC2(6,3) XC2(8,3) XC2(8,4) XC2(6,4)]+110];
YC1i=[YC1(6,1:2) YC1(8,2) YC1(8,1) YC1(6,1);YC2(6,4) YC2(6,3) YC2(8,3) YC2(8,4) YC2(6,4)];
ZC1i=[ZC1(6,1:2) ZC1(8,2) ZC1(8,1) ZC1(6,1);ZC2(6,4) ZC2(6,3) ZC2(8,3) ZC2(8,4) ZC2(6,4)];

%GRAF COLUMN2
 YY1=[0;r(144:176)]*fac;
 ZZ1=[0;r(210:242)]*fac;
XYZ3=[XYZc2(:,1) XYZc2(:,2)+YY1 XYZc2(:,3)+ZZ1];
for i=2:(length(XYZ3(:,1)))
    p1=(XYZ3(i,:)-XYZ3(i-1,:));
    angle3(i-1)=-atan(p1(1,2)/p1(1,3));
end
Angle3(1)=0;
for i=2:length(angle3)
    Angle3(i)=(angle3(i-1)+angle3(i))/2;
end  
Angle3(i+1)=angle3(end);
facx=[ones(1,8),linspace(1,0.8182,17),ones(1,9)*0.8182];
ampx=[linspace(180,290,6),290,linspace(290,45,18),ones(1,9)*45];

for i=1:length(Angle3)
    th=Angle3(i);
    MT3=[1 0 0;0 cos(th) -sin(th);0 sin(th) cos(th)];
    for j=1:5
        VR3=MT3*[Xc1(j)*facx(i);Yc1(j);Zc(j)];
       XC12(i,j)=VR3(1)+XYZ3(i,1)+ampx(i);
       YC12(i,j)=VR3(2)+XYZ3(i,2);
       ZC12(i,j)=VR3(3)+XYZ3(i,3);
       XC22(i,j)=VR3(1)+XYZ3(i,1)-ampx(i);
       YC22(i,j)=VR3(2)+XYZ3(i,2);
       ZC22(i,j)=VR3(3)+XYZ3(i,3);
   end
end
XC2i=[[XC12(6,1:2) XC12(8,2) XC12(8,1) XC12(6,1)]-110;[XC22(6,4) XC22(6,3) XC22(8,3) XC22(8,4) XC22(6,4)]+110];
YC2i=[YC12(6,1:2) YC12(8,2) YC12(8,1) YC12(6,1);YC22(6,4) YC22(6,3) YC22(8,3) YC22(8,4) YC22(6,4)];
ZC2i=[ZC12(6,1:2) ZC12(8,2) ZC12(8,1) ZC12(6,1);ZC22(6,4) ZC22(6,3) ZC22(8,3) ZC22(8,4) ZC22(6,4)];

XYZT=[XYZt(:,1) XYZt(:,2) XYZt(:,3)+ZZ0];
XYZC1=[XYZc1(:,1) XYZc1(:,2)+YY XYZc1(:,3)+ZZ];
XYZC2=[XYZc2(:,1) XYZc2(:,2)+YY1 XYZc2(:,3)+ZZ1];
col=[zeros(50,1),linspace(0,1,50)',ones(50,1);
     zeros(50,1),ones(50,1),linspace(1,0,50)';
     linspace(0,1,50)',ones(50,1),zeros(50,1);
     ones(51,1),linspace(1,0,51)',zeros(51,1)];
 n=length(col(:,1));

XA1=[40 40 -40 -40 40];
YA1=[-40 40 40 -40 -40];
ZA1=[1 1 1 1 1];
XA2=[1 1 1 1 1];
YA2=[-40 40 40 -40 -40];
ZA2=[-40 -40 40 40 -40];
Xa1=[[XA1;XA1;XA1;XA1;XA1]-190];
Ya1=[YA1;YA1;YA1;YA1;YA1];
Za1=[ZA1*100;ZA1*200;ZA1*300;ZA1*400;ZA1*500];
Xa2=[[XA1;XA1;XA1;XA1;XA1]+190];
Ya2=[YA1;YA1;YA1;YA1;YA1];
Za2=[ZA1*500;ZA1*400;ZA1*300;ZA1*200;ZA1*100];
Xa3=[XA2*-230;XA2*-150;XA2*-50;XA2*50;XA2*150;XA2*230];
Ya3=[YA2;YA2;YA2;YA2;YA2;YA2];
Za3=[ZA2;ZA2;ZA2;ZA2;ZA2;ZA2]+540;
XA=[Xa1;Xa3;Xa2];
YA=[Ya1;Ya3;Ya2];
ZA=[Za1;Za3;Za2];