	  for C1 = 0:7
              for C0 = 0:7 
                if mod(C0*S1L+S0L,7)==0 && mod(C0*S2L+S1L,7)==0 && mod(C0*S3L+S2L,7)==0;
                    for C1F= 1:7
                        if mod((C0*C1F),7)==1
                            if S0L==0 && S1L==0 && S2L==0 && S3L==0
                                C1F=0;
                            end
                            cx=(1+C1F*x);
                            F0L=mod(-(C1F*S3L),7);
                            F1L=mod(-(C1F*F0L),7);
                        end
                    end
                elseif mod(C0*S2L+C1*S1L+S0L,7)==0 && mod(C0*S3L+C1*S2L+S1L,7)==0;
                    for C2= 1:7
                        if mod((C0*C2),7)==1
                            C1=mod(C1*C2,7)
                            cx=(1+C1*x+C2*x^2);
                            F0L=mod(-(C1*S3L+C2*S2L),7);
                            F1L=mod(-(C1*F0L+C2*S3L),7);
                        end
                    end
                end 
             end   
            end									