Объясните построение тетраэдра по двум зад точкам

Program ger02;
Uses CRT, Graphabc;
  Type TR= array[1..3000] of real; 
    Procedure PR(X, Y, Z, Zpr, Zf: real; Var Xpi, Ypi: real);
       Var t: real; 
         Begin
           t:=(Zpr-Z)/(Z-Zf); 
           Xpi:=X+X*t;       
           Ypi:=Y+Y*t;        
         end;   
  procedure CIRC(Var a,b:TR; N:integer; fi:real);
      Var r,c,s:real; i:integer;
       Begin
         c:=cos(fi);
         s:=sin(fi);
      for i:=1 to n do
       begin
         r:=c*a[i]-s*b[i];
         b[i]:=a[i]*s+b[i]*c;
         a[i]:=r;
       end;
       end;
 Var t,x, y, z, Xs, Ys: TR;
     fi, Zpr, Zf,  xc, yc, zc, Rk, k, Pk,
     Hx,Hy,Hz,sx,sy,sz,px,py,ax,ay,az,rx,ry,rz,pz,qx,qy,qz,nx,ny,nz, fx,fy,fz,
     nk,mx,my,mz,s,l,r: real;
     i,  N, L0,tx,ty: integer; kod:char;
      BEGIN    [COLOR="#8b0000"]{программное присвоение исходных данных}[/COLOR]
 
 
    N:=8;
      x[1]:=0; x[2]:=0;
        y[1]:=0; y[2]:=0;
         z[1]:=2; z[2]:=6;
  [COLOR="DarkOrange"]  l:=sqrt((x[2]-x[1])*(x[2]-x[1])+(y[2]-y[1])*
             (y[2]-y[1])+(z[2]-z[1])*(z[2]-z[1]));
           ax:=(x[2]-x[1]);
           ay:=(y[2]-y[1]);
           az:=(z[2]-z[1]);
 px:=(x[1]+x[2])/2; py:=(y[1]+y[2])/2; pz:=(z[1]+z[2])/2;
       if az=0 then  begin rx:=0; ry:=0; rz:=1  end
       else   begin ry:=-az;  rx:=-az; rz:=ax+ay end;
 
        s:=sqrt(rx*rx+ry*ry+rz*rz);
 
          mx:=rx/s*sqrt(3)*l/2;
          my:=ry/s*sqrt(3)*l/2;
          mz:=rz/s*sqrt(3)*l/2;
         x[3]:=px+mx;
         y[3]:=py+my;
         z[3]:=pz+mz;
           sx:=(x[1]+x[2]+x[3])/3;
           sy:=(y[1]+y[2]+y[3])/3;
           sz:=(z[1]+z[2]+z[3])/3;
      fx:=(x[1]-sx); fy:=(y[1]-sy); fz:=(z[1]-sz);
      qx:=(x[3]-sx); qy:=(y[3]-sy);qz:=(z[3]-sz);
         nx:=fy*qz-fz*qy;
         ny:=-fx*qz+fz*qx;
         nz:=fx*qy-fy*qx;
      nk:=sqrt(nx*nx+ny*ny+nz*nz);
          Hx:=nx/nk*l*sqrt(2/3);
          Hy:=ny/nk*l*sqrt(2/3);
          Hz:=nz/nk*l*sqrt(2/3);
      x[4]:=sx+Hx;
       y[4]:=sy+Hy;
         z[4]:=sz+Hz;[/COLOR]
    x[5]:=x[1];   y[5]:=y[1];  z[5]:=z[1];
    x[6]:=x[3];   y[6]:=y[3];  z[6]:=z[3];
    x[7]:=x[4];   y[7]:=y[4];  z[7]:=z[4];
    x[8]:=x[2];   y[8]:=y[2];  z[8]:=z[2];
 
 L0:=500;
 
SetWindowSize(L0,L0);
 
             xc:=0; yc:=0; zc:=0;
            for i:=1 to N do
             begin
               xc:=xc+x[i]/N;
                yc:=yc+y[i]/N;
                 zc:=zc+z[i]/N;
             end;
            for i:=1 to N do
             begin
               x[i]:=x[i]-xc;
               y[i]:=y[i]-yc;
               z[i]:=z[i]-zc;
             end;
              Rk:=0;
          for i:=1 to N do
            if Rk<sqrt(x[i]*x[i]+y[i]*y[i]+z[i]*z[i]) then
               Rk:=sqrt(x[i]*x[i]+y[i]*y[i]+z[i]*z[i]);
              Zpr:=rk*5; Zf:=rk*4;
            Pk:=Rk*(Zf-Zpr)/Sqrt(Zf*Zf-Rk*Rk);
            k:=0.5*L0/Pk;
 
while true do
begin   изображения}[/COLOR]
            
 for i:=1 to N do PR(x[i], y[i], z[i], Zpr, Zf, Xs[i], Ys[i]);
 
 Setpencolor(RGB(0,0,200));  
          for i:=1 to N-1 do    
                Line(Round((Xs[i]+Pk)*k),    Round((Pk-Ys[i])*k),
                Round((Xs[i+1]+Pk)*k),   Round((Pk-Ys[i+1])*k));
   kod:=ReadKey;
     fi:=0.11;
     case  kod  of
     #27:exit;
      'a':CIRC(x,y,N,fi);
      'z':CIRC(x,y,N,-fi);
      's':CIRC(y,z,N,fi);
      'x':CIRC(z,y,N,fi);    
      'd':CIRC(x,z,N,fi);
      'c':CIRC(x,z,N,-fi);
      end;
       clearwindow;
end;
END.