我编了下面一个程序，是二叉树结构划分的，为什么每次调用就出现下面的错误 One or more output arguments not assigned during call to 'parsubsample'？请求高手指点！谢谢 

 function [Y1,f,Y2,g,Y3,h,Y4,l,Y5,o,Y6,p,Y7,q,Y8,r,P,c]=parsubsample(N,z1,z2,z3,m) 
    %the order of the splitting variable both in the left and right branch is z1,z2,z3
   %to find the corresponding partion of the tree and the corresponding
   %subsamples when svar(1,8)=z1,svar(2,4)=svar(2,,12)=z2,svar(3,2)=svar(3,6)=svar(3,10)=svar(3,14)=z3
  
 
 
   m=length(z1);
    %m is the number of elements in z1, z2 and z3
 
   c(1)=0;
    %to find the corresponding partion of the tree and the corresponding
   %subsamples when svar(1,8)=z1,svar(2,4)=svar(2,,12)=z2,svar(3,2)=svar(3,6)=svar(3,10)=svar(3,14)=z3
   if N(1,8)=='terminal' 
       %if the root node is the terminal node
       P=sym('[ ]');
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           f=f+1;
           Y1(f,:)=[z1(t),z2(t),z3(t)];
           c(t)=1;
       end
       Y2=zeros(0,3);Y3=zeros(0,3);Y4=zeros(0,3);Y5=zeros(0,3);
       Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
   end
  
      
   %there is no terminal node in the first row of the tree
   if N(2,4)=='terminal'&N(2,12)=='terminal'
       %if the left and right nodes in the second row are both terminal nodes
       P1=sym('[z1<=0]');
       P2=sym('[z1>0]');
       P=[P1,P2];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;      
       for t=2:m
           if z1(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           else
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           end
        end
       Y3=zeros(0,3);Y4=zeros(0,3);Y5=zeros(0,3);
       Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
   end
  
  if N(2,4)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
      %if the left in the second row and the third and fourth node in the third  row are all terminal nodes
       P1=sym('[z1<=0]');
       P2=sym('[z1>0z2<=0]');
       P3=sym('[z1>0z2>0]');
       P=[P1,P2,P3];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;     
       for t=2:m
           if z1(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)>0&z2(t-1)<=0
                g=g+1;
                Y2(g,:)=[z1(t),z2(t),z3(t)];
                c(t)=2;
           elseif z1(t-1)>0&z2(t-1)>0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
           end
       end
       Y4=zeros(0,3);Y5=zeros(0,3);Y6=zeros(0,3);
       Y7=zeros(0,3);Y8=zeros(0,3);
  end 
  
   if N(2,4)=='terminal'&N(4,9)=='terminal'&N(3,14)=='terminal'
      %if the left in the second row ,the fourth node in the third row and the fifth and sixth node in the fourth row are all terminal nodes
       P1=sym('[z1<=0]');
       P2=sym('[z1>0z2<=0z3<=0]');
       P3=sym('[z1>0z2<=0z3>0]');
       P4=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;     
       for t=2:m
           if z1(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
            elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
           elseif z1(t-1)>0&z2(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
           end
       end
       Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
   end 
 
    if N(2,4)=='terminal'&N(3,10)=='terminal'&N(4,13)=='terminal'
      %if the left in the second row ,the third node in the third row and the seventh and eighth node in the fourth row are all terminal nodes
       P1=sym('[z1<=0]');
       P2=sym('[z1>0z2<=0]');
       P3=sym('[z1>0z2>0z3<=0]');
       P4=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;     
       for t=2:m
           if z1(t-1)<=0
                f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
              c(t)=1;
           elseif z1(t-1)>0&z2(t-1)<=0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
            elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
           elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
           end
       end
       Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
    end 
  
    if N(2,4)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
      %if the left in the second row ,the fifth,sixth,seventh and eighth node in the fourth row are all terminal nodes
       P1=sym('[z1<=0]');
       P2=sym('[z1>0z2<=0z3<=0]');
       P3=sym('[z1>0z2<=0z3>0]');
       P4=sym('[z1>0z2>0z3<=0]');
       P5=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;      
       for t=2:m
           if z1(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
            elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
           elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
           elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
                o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
              c(t)=5;
           end
       end
       Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
    end  
   
   
    if N(2,12)=='terminal'&N(3,2)=='terminal'&N(3,6)=='terminal'
      %if the right in the second row and the first and second node in the third  row are all terminal nodes
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0]');
       P3=sym('[z1>0]');
       P=[P1,P2,P3];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;     
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)>0
                h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
           end
       end
       Y4=zeros(0,3);Y5=zeros(0,3);Y6=zeros(0,3);
       Y7=zeros(0,3);Y8=zeros(0,3);
    end 
 
    if N(2,12)=='terminal'&N(4,1)=='terminal'&N(3,6)=='terminal'
      %if the right in the second row ,the second node in the third row and the first and second node in the fourth row are all terminal nodes
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0]');
       P4=sym('[z1>0]');
       P=[P1,P2,P3,P4];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;      
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
            elseif z1(t-1)<=0&z2(t-1)>0
                h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
           elseif z1(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
           end
       end
       Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
    end 
  
   if N(2,12)=='terminal'&N(4,5)=='terminal'&N(3,2)=='terminal'
      %if the right in the second row ,the first node in the third row and the third and fourth node in the fourth row are all terminal nodes
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0z3<=0]');
       P3=sym('[z1<=0z2>0z3>0]');
       P4=sym('[z1>0]');
       P=[P1,P2,P3,P4];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;      
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
            elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
           elseif z1(t-1)>0
                l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
              c(t)=4;
           end
       end
       Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
   end

   if N(2,12)=='terminal'&N(4,1)=='terminal'&N(4,5)=='terminal'
      %if the right in the second row ,the first,second,third and fourth node in the fourth row are all terminal nodes
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0z3<=0]');
       P4=sym('[z1<=0z2>0z3>0]');
       P5=sym('[z1>0]');
       P=[P1,P2,P3,P4,P5];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;     
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
            elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
                h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
                l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
              c(t)=4;
           elseif z1(t-1)>0
                o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
              c(t)=5;
           end
       end
       Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
   end

  
   % there is no terminal node in the first and second row of the tree
   if N(3,2)=='terminal'&N(3,6)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
       %if the first,second,third and fourth node in the third row are
       %terminal node
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0]');
       P3=sym('[z1>0z2<=0]');
       P4=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0
                g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
              c(t)=2;
          elseif z1(t-1)>0&z2(t-1)<=0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)>0&z2(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
           end
       end
       Y5=zeros(0,3);Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
   end
      
   if N(4,1)=='terminal'&N(3,6)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
       %if the second,third and fourth node in the third row and the first and second node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0]');
       P4=sym('[z1>0z2<=0]');
       P5=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4,P5];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
              g=g+1;
                Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
          elseif z1(t-1)<=0&z2(t-1)>0
                h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)>0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
           end
       end
       Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
   end
  
    if N(4,5)=='terminal'&N(3,2)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
       %if the first,third and fourth node in the third row and the third and fourth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0z3<=0]');
       P3=sym('[z1<=0z2>0z3>0]');
       P4=sym('[z1>0z2<=0]');
       P5=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4,P5];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
              f=f+1;
                Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
          elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
                h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)>0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
           end
       end
       Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
    end
  
    if N(4,9)=='terminal'&N(3,2)=='terminal'&N(3,6)=='terminal'&N(3,14)=='terminal'
       %if the first,second and fourth node in the third row and the fifth and sixth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0]');
       P3=sym('[z1>0z2<=0z3<=0]');
       P4=sym('[z1>0z2<=0z3>0]');
       P5=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4,P5];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
              f=f+1;
                Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)>0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
           end
       end
       Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
    end 
   
     if N(4,13)=='terminal'&N(3,2)=='terminal'&N(3,6)=='terminal'&N(3,10)=='terminal'
       %if the first,second and third node in the third row and the seventh and eighth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0]');
       P3=sym('[z1>0z2<=0]');
       P4=sym('[z1>0z2>0z3<=0]');
       P5=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
          elseif z1(t-1)>0&z2(t-1)<=0
                h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
           end
       end
       Y6=zeros(0,3);Y7=zeros(0,3);Y8=zeros(0,3);
     end

     if N(4,1)=='terminal'&N(4,5)=='terminal'&N(3,10)=='terminal'&N(3,14)=='terminal'
       %if the third and fourth node in the third row and the first,second,third and fourth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0z3<=0]');
       P4=sym('[z1<=0z2>0z3>0]');
       P5=sym('[z1>0z2<=0]');
       P6=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4,P5,P6];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)<=0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
           end
       end
       Y7=zeros(0,3);Y8=zeros(0,3);
     end 
    
      if N(4,1)=='terminal'&N(4,9)=='terminal'&N(3,6)=='terminal'&N(3,14)=='terminal'
       %if the second and fourth node in the third row and the first,second,fifth and sixth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0]');
       P4=sym('[z1>0z2<=0z3<=0]');
       P5=sym('[z1>0z2<=0z3>0]');
       P6=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4,P5,P6];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
                g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
              c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0
                h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0
                p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
              c(t)=6;
           end
       end
       Y7=zeros(0,3);Y8=zeros(0,3);
      end 
    
      if N(4,1)=='terminal'&N(4,13)=='terminal'&N(3,6)=='terminal'&N(3,10)=='terminal'
       %if the second and third node in the third row and the first,second,seventh and eighth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0]');
       P4=sym('[z1>0z2<=0]');
       P5=sym('[z1>0z2>0z3<=0]');
       P6=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5,P6];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
           end
       end
       Y7=zeros(0,3);Y8=zeros(0,3);
      end 
     
      if N(3,2)=='terminal'&N(4,5)=='terminal'&N(4,9)=='terminal'&N(3,14)=='terminal'
       %if the first and fourth node in the third row and the third,fourth,fifth and sixth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0z3<=0]');
       P3=sym('[z1<=0z2>0z3>0]');
       P4=sym('[z1>0z2<=0z3<=0]');
       P5=sym('[z1>0z2<=0z3>0]');
       P6=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4,P5,P6];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
                o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
              c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
           end
       end
       Y7=zeros(0,3);Y8=zeros(0,3);
      end 
     
      if N(3,2)=='terminal'&N(4,5)=='terminal'&N(3,10)=='terminal'&N(4,13)=='terminal'
       %if the first and third node in the third row and the third,fourth,seventh and eighth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0z3<=0]');
       P3=sym('[z1<=0z2>0z3>0]');
       P4=sym('[z1>0z2<=0]');
       P5=sym('[z1>0z2>0z3<=0]');
       P6=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5,P6];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
                g=g+1;
                Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
           end
       end
       Y7=zeros(0,3);Y8=zeros(0,3);
      end 
     
      if N(3,2)=='terminal'&N(3,6)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
       %if the first and second node in the third row and the fifth,sixth,seventh and eighth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0]');
       P3=sym('[z1>0z2<=0z3<=0]');
       P4=sym('[z1>0z2<=0z3>0]');
       P5=sym('[z1>0z2>0z3<=0]');
       P6=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5,P6];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               h=h+1;
                Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
           end
       end
       Y7=zeros(0,3);Y8=zeros(0,3);
      end
     
     
      if N(3,2)=='terminal'&N(4,5)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
       %if the first node in the third row and the third,fourth,fifth,sixth,seventh and eighth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0]');
       P2=sym('[z1<=0z2>0z3<=0]');
       P3=sym('[z1<=0z2>0z3>0]');
       P4=sym('[z1>0z2<=0z3<=0]');
       P5=sym('[z1>0z2<=0z3>0]');
       P6=sym('[z1>0z2>0z3<=0]');
       P7=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5,P6,P7];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0
              f=f+1;
                Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               q=q+1;  
               Y7(q,:)=[z1(t),z2(t),z3(t)];
               c(t)=7;
           end
       end
       Y8=zeros(0,3);
      end
     
      if N(4,1)=='terminal'&N(3,6)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
       %if the second node in the third row and the first,second,fifth,sixth,seventh and eighth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0]');
       P4=sym('[z1>0z2<=0z3<=0]');
       P5=sym('[z1>0z2<=0z3>0]');
       P6=sym('[z1>0z2>0z3<=0]');
       P7=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5,P6,P7];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0
                h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
              c(t)=3;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               q=q+1;   
               Y7(q,:)=[z1(t),z2(t),z3(t)];
               c(t)=7;
           end
       end
       Y8=zeros(0,3);
      end
     
      if N(4,1)=='terminal'&N(4,5)=='terminal'&N(3,10)=='terminal'&N(4,13)=='terminal'
       %if the third node in the third row and the first,second,third,fourth,seventh and eighth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0z3<=0]');
       P4=sym('[z1<=0z2>0z3>0]');
       P5=sym('[z1>0z2<=0]');
       P6=sym('[z1>0z2>0z3<=0]');
       P7=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5,P6,P7];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
                f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
              c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)<=0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
                p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
              c(t)=6;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               q=q+1;  
               Y7(q,:)=[z1(t),z2(t),z3(t)];
               c(t)=7;
           end
       end
       Y8=zeros(0,3);
      end
     
      if N(4,1)=='terminal'&N(4,5)=='terminal'&N(4,9)=='terminal'&N(3,14)=='terminal'
       %if the fourth node in the third row and the first,second,third,fourth,fifth and sixth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0z3<=0]');
       P4=sym('[z1<=0z2>0z3>0]');
       P5=sym('[z1>0z2<=0z3<=0]');
       P6=sym('[z1>0z2<=0z3>0]');
       P7=sym('[z1>0z2>0]');
       P=[P1,P2,P3,P4,P5,P6,P7];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
               l=l+1;
              Y4(l,:)=[z1(t),z2(t),z3(t)];
                c(t)=4;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
          elseif z1(t-1)>0&z2(t-1)>0
               q=q+1; 
               Y7(q,:)=[z1(t),z2(t),z3(t)];
               c(t)=7; 
           end
       end
       Y8=zeros(0,3);
      end
     
     
      %there is no terminal node in the first,second and third row in the
      %tree
       if N(4,1)=='terminal'&N(4,5)=='terminal'&N(4,9)=='terminal'&N(4,13)=='terminal'
       %if the first,second,third,fourth,fifth,sixth,seventh and eighth node in the fourth row are
       %terminal node
       P1=sym('[z1<=0z2<=0z3<=0]');
       P2=sym('[z1<=0z2<=0z3>0]');
       P3=sym('[z1<=0z2>0z3<=0]');
       P4=sym('[z1<=0z2>0z3>0]');
       P5=sym('[z1>0z2<=0z3<=0]');
       P6=sym('[z1>0z2<=0z3>0]');
       P7=sym('[z1>0z2>0z3<=0]');
       P8=sym('[z1>0z2>0z3>0]');
       P=[P1,P2,P3,P4,P5,P6,P7,P8];
       f=0;g=0;h=0;l=0;o=0;p=0;q=0;r=0;
       for t=2:m
           if z1(t-1)<=0&z2(t-1)<=0&z3(t-1)<=0
               f=f+1;
               Y1(f,:)=[z1(t),z2(t),z3(t)];
               c(t)=1;
           elseif z1(t-1)<=0&z2(t-1)<=0&z3(t-1)>0
               g=g+1;
               Y2(g,:)=[z1(t),z2(t),z3(t)];
               c(t)=2;
           elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)<=0
               h=h+1;
               Y3(h,:)=[z1(t),z2(t),z3(t)];
               c(t)=3;
          elseif z1(t-1)<=0&z2(t-1)>0&z3(t-1)>0
               l=l+1;
               Y4(l,:)=[z1(t),z2(t),z3(t)];
               c(t)=4;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)<=0
               o=o+1;
               Y5(o,:)=[z1(t),z2(t),z3(t)];
               c(t)=5;
          elseif z1(t-1)>0&z2(t-1)<=0&z3(t-1)>0
               p=p+1;
               Y6(p,:)=[z1(t),z2(t),z3(t)];
               c(t)=6;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)<=0
               q=q+1;   
               Y7(q,:)=[z1(t),z2(t),z3(t)];
               c(t)=7;
          elseif z1(t-1)>0&z2(t-1)>0&z3(t-1)>0
               r=r+1;
              Y8(r,:)=[z1(t),z2(t),z3(t)];
                c(t)=8;
           end
       end
   end