partialX(C,a,b,c,d,e,f) =
{return(subst(subst(subst(subst(subst(subst(deriv(C,x),x,a),y,b),z,c),u,d),v,e),w,f))}
partialY(C,a,b,c,d,e,f) =
{return(subst(subst(subst(subst(subst(subst(deriv(C,y),x,a),y,b),z,c),u,d),v,e),w,f))}
partialZ(C,a,b,c,d,e,f) =
{return(subst(subst(subst(subst(subst(subst(deriv(C,z),x,a),y,b),z,c),u,d),v,e),w,f))}
partialU(C,a,b,c,d,e,f) =
{return(subst(subst(subst(subst(subst(subst(deriv(C,u),x,a),y,b),z,c),u,d),v,e),w,f))}
partialV(C,a,b,c,d,e,f) =
{return(subst(subst(subst(subst(subst(subst(deriv(C,v),x,a),y,b),z,c),u,d),v,e),w,f))}
partialW(C,a,b,c,d,e,f) =
{return(subst(subst(subst(subst(subst(subst(deriv(C,w),x,a),y,b),z,c),u,d),v,e),w,f))}

find_singular(C,p) = {

a=0;
b=0;
c=0;
d=0;
e=0;
f=1;

while(a<=p-1,
  while(b<=p-1,
    while(c<=p-1,
   while(d<=p-1,
     while(e<=p-1,
       while(f<=p-1,
          if(partialX(C,a,b,c) %p == 0 &&
                partialY(C,a,b,c,d,e,f) %p == 0 &&
                partialZ(C,a,b,c,d,e,f) %p == 0 &&
                partialU(C,a,b,c,d,e,f) %p == 0 &&
                partialV(C,a,b,c,d,e,f) %p == 0 &&
                partialW(C,a,b,c,d,e,f) %p == 0,
                print("(",a,",",b,",",c,",",d,",",e,",",f,")")
               ); 
            f++);
         e++;f=0);
      d++;e=0;f=0);
   c++;d=0;e=0;f=0);
  b++;c=0;d=0;e=0;f=0);
a++;b=0;c=0;d=0;e=0;f=0;
)
}