Hyperbolic section 

from P. Olver, Pick invariant, in preparation. 

fr > Section:=[x=0,y=0,u[0,0]=0,u[1,0]=0,u[0,1]=0,u[1,1]=0,u[2,0]=1,u[0,2]=-1,u[0,3]=0,u[2,1]=0, u[1,2]=u[3,0]]:stair(Section);
 

Plot 

fr > C, Invariantizations :=dmf(LieAlg,Section, fr, [op(1..10,Section), u[3,0]=Pk,u[4,0]=e-3*c,u[3,1]=d,u[2,2]=c,u[1,3]=-12*Pk*S1-3*d,u[0,4]=12*Pk*S2-3*c],'COM'):
 

"Section :" 

[u[0, 3] = 0, u[1, 2] = u[3, 0], u[2, 1] = 0, u[0, 2] = -1, u[1, 1] = 0, u[2, 0] = 1, u[0, 1] = 0, u[1, 0] = 0, u[0, 0] = 0, y = 0, x = 0]
[u[0, 3] = 0, u[1, 2] = u[3, 0], u[2, 1] = 0, u[0, 2] = -1, u[1, 1] = 0, u[2, 0] = 1, u[0, 1] = 0, u[1, 0] = 0, u[0, 0] = 0, y = 0, x = 0]
 

"Transversality condition :"   -u[3   0]
"Invariantizations:"
 

x = 0, y = 0, u[0, 0] = 0, u[1, 0] = 0, u[0, 1] = 0, u[1, 1] = 0, u[2, 0] = 1, u[0, 2] = -1, u[0, 3] = 0, u[2, 1] = 0, u[3, 0] = Pk, u[4, 0] = e-3*c, u[3, 1] = d, u[2, 2] = c, u[1, 3] = -12*Pk*S1-3*d,...
x = 0, y = 0, u[0, 0] = 0, u[1, 0] = 0, u[0, 1] = 0, u[1, 1] = 0, u[2, 0] = 1, u[0, 2] = -1, u[0, 3] = 0, u[2, 1] = 0, u[3, 0] = Pk, u[4, 0] = e-3*c, u[3, 1] = d, u[2, 2] = c, u[1, 3] = -12*Pk*S1-3*d,...
x = 0, y = 0, u[0, 0] = 0, u[1, 0] = 0, u[0, 1] = 0, u[1, 1] = 0, u[2, 0] = 1, u[0, 2] = -1, u[0, 3] = 0, u[2, 1] = 0, u[3, 0] = Pk, u[4, 0] = e-3*c, u[3, 1] = d, u[2, 2] = c, u[1, 3] = -12*Pk*S1-3*d,...
 

"Commutation rules" 

[D1, D2] = -1/12*(3*u[3, 1]+u[1, 3])*D1/u[3, 0]-1/12*(-u[0, 4]-3*u[2, 2])*D2/u[3, 0] 

"Syzygies" 

[-36*Pk[0, 0]^2*S1[0, 0]-9*Pk[0, 0]*d[0, 0]-4*d[0, 0]*S2[0, 0]-e[0, 0]*S1[0, 0]-e[0, 1]+3*c[0, 1]+d[1, 0], -S2[0, 0]*e[0, 0]-6*Pk[0, 0]^2*S2[0, 0]+4*d[0, 0]*S1[0, 0]+24*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*...
[-36*Pk[0, 0]^2*S1[0, 0]-9*Pk[0, 0]*d[0, 0]-4*d[0, 0]*S2[0, 0]-e[0, 0]*S1[0, 0]-e[0, 1]+3*c[0, 1]+d[1, 0], -S2[0, 0]*e[0, 0]-6*Pk[0, 0]^2*S2[0, 0]+4*d[0, 0]*S1[0, 0]+24*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*...
[-36*Pk[0, 0]^2*S1[0, 0]-9*Pk[0, 0]*d[0, 0]-4*d[0, 0]*S2[0, 0]-e[0, 0]*S1[0, 0]-e[0, 1]+3*c[0, 1]+d[1, 0], -S2[0, 0]*e[0, 0]-6*Pk[0, 0]^2*S2[0, 0]+4*d[0, 0]*S1[0, 0]+24*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*...
[-36*Pk[0, 0]^2*S1[0, 0]-9*Pk[0, 0]*d[0, 0]-4*d[0, 0]*S2[0, 0]-e[0, 0]*S1[0, 0]-e[0, 1]+3*c[0, 1]+d[1, 0], -S2[0, 0]*e[0, 0]-6*Pk[0, 0]^2*S2[0, 0]+4*d[0, 0]*S1[0, 0]+24*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*...
[-36*Pk[0, 0]^2*S1[0, 0]-9*Pk[0, 0]*d[0, 0]-4*d[0, 0]*S2[0, 0]-e[0, 0]*S1[0, 0]-e[0, 1]+3*c[0, 1]+d[1, 0], -S2[0, 0]*e[0, 0]-6*Pk[0, 0]^2*S2[0, 0]+4*d[0, 0]*S1[0, 0]+24*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*...
[-36*Pk[0, 0]^2*S1[0, 0]-9*Pk[0, 0]*d[0, 0]-4*d[0, 0]*S2[0, 0]-e[0, 0]*S1[0, 0]-e[0, 1]+3*c[0, 1]+d[1, 0], -S2[0, 0]*e[0, 0]-6*Pk[0, 0]^2*S2[0, 0]+4*d[0, 0]*S1[0, 0]+24*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*...
[-36*Pk[0, 0]^2*S1[0, 0]-9*Pk[0, 0]*d[0, 0]-4*d[0, 0]*S2[0, 0]-e[0, 0]*S1[0, 0]-e[0, 1]+3*c[0, 1]+d[1, 0], -S2[0, 0]*e[0, 0]-6*Pk[0, 0]^2*S2[0, 0]+4*d[0, 0]*S1[0, 0]+24*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*...
[-36*Pk[0, 0]^2*S1[0, 0]-9*Pk[0, 0]*d[0, 0]-4*d[0, 0]*S2[0, 0]-e[0, 0]*S1[0, 0]-e[0, 1]+3*c[0, 1]+d[1, 0], -S2[0, 0]*e[0, 0]-6*Pk[0, 0]^2*S2[0, 0]+4*d[0, 0]*S1[0, 0]+24*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*...
 

"Syzygies with attempt to eliminate the y's" 

y0[0, 0] = -1/Pk[0, 0] 

y1[0, 0] = Pk[0, 0] 

d[1, 0] = 36*Pk[0, 0]^2*S1[0, 0]+9*Pk[0, 0]*d[0, 0]+4*d[0, 0]*S2[0, 0]+e[0, 0]*S1[0, 0]+e[0, 1]-3*c[0, 1]
d[1, 0] = 36*Pk[0, 0]^2*S1[0, 0]+9*Pk[0, 0]*d[0, 0]+4*d[0, 0]*S2[0, 0]+e[0, 0]*S1[0, 0]+e[0, 1]-3*c[0, 1]
 

c[1, 0] = S2[0, 0]*e[0, 0]+6*Pk[0, 0]^2*S2[0, 0]-4*d[0, 0]*S1[0, 0]-24*Pk[0, 0]*S1[0, 0]^2-3*Pk[0, 0]*c[0, 0]+d[0, 1]
c[1, 0] = S2[0, 0]*e[0, 0]+6*Pk[0, 0]^2*S2[0, 0]-4*d[0, 0]*S1[0, 0]-24*Pk[0, 0]*S1[0, 0]^2-3*Pk[0, 0]*c[0, 0]+d[0, 1]
 

Pk[1, 0] = 1/4*e[0, 0] 

S1[1, 0] = 1/12*(-8*d[0, 0]*S2[0, 0]-96*Pk[0, 0]^2*S1[0, 0]-24*Pk[0, 0]*d[0, 0]+12*Pk[0, 0]*S1[0, 0]*S2[0, 0]+8*c[0, 1]-6*e[0, 0]*S1[0, 0]-3*e[0, 1])/Pk[0, 0]
S1[1, 0] = 1/12*(-8*d[0, 0]*S2[0, 0]-96*Pk[0, 0]^2*S1[0, 0]-24*Pk[0, 0]*d[0, 0]+12*Pk[0, 0]*S1[0, 0]*S2[0, 0]+8*c[0, 1]-6*e[0, 0]*S1[0, 0]-3*e[0, 1])/Pk[0, 0]
 

S2[1, 0] = -c[0, 0]+S2[0, 0]^2+3/2*Pk[0, 0]*S2[0, 0]-S1[0, 0]^2+1/8*e[0, 0]-S1[0, 1] 

Pk[0, 1] = -2*d[0, 0]-9*Pk[0, 0]*S1[0, 0] 

fr > COM;
 

[[y, x] = [-S1[0, 0], -S2[0, 0]]] 

We can express [c,d,e], the invariantization of [u[2, 2], u[3, 1], u[4, 0]]in terms of [S1, S2, Pk] the invariantizations of  

[1/12*Typesetting:-mo( 

> R := differential_ring(ranking=[[y0,y1],[c,d,e],[S2,S1,Pk]], derivations=[x,y],
          commutations=COM,notation=vjet);
CC := Rosenfeld_Groebner(equations(C[1]), [Pk[0,0]], R);
for g in rewrite_rules(CC[1]) do print(g); od;
 

(Typesetting:-mprintslash)([R := differential_polynomial_ring], [differential_polynomial_ring]) 

(Typesetting:-mprintslash)([CC := [characterisable]], [[characterisable]]) 

y0[0, 0] = -1/Pk[0, 0] 

y1[0, 0] = Pk[0, 0] 

c[0, 0] = S2[0, 0]^2+3/2*Pk[0, 0]*S2[0, 0]-S1[0, 0]^2+1/2*Pk[1, 0]-S1[0, 1]-S2[1, 0] 

d[0, 0] = -1/2*Pk[0, 1]-9/2*Pk[0, 0]*S1[0, 0] 

e[0, 0] = 4*Pk[1, 0] 

Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
Pk[2, 1] = -1/3*Pk[0, 3]+3*Pk[0, 0]^2*Pk[0, 1]+8/3*Pk[1, 0]*S2[0, 1]+4*Pk[0, 0]^2*S2[0, 1]-2*Pk[0, 2]*S1[0, 0]-8/3*Pk[0, 1]*S1[0, 1]-11/3*Pk[0, 1]*S1[0, 0]^2-3*Pk[0, 0]*S1[0, 2]-S1[0, 0]*Pk[2, 0]+5*S2...
 

S2[2, 0] = -5/2*Pk[1, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+5/2*Pk[0, 1]*S1[0, 0]+3*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*S2[0, 0]^2+3/2*Pk[0, 0]*Pk[1, 0]+3/2*Pk[0, 0]*S1[0, 1]-3/2*Pk[0, 0]*S2[1, 0]+1/2*Pk[0, ...
S2[2, 0] = -5/2*Pk[1, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+5/2*Pk[0, 1]*S1[0, 0]+3*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*S2[0, 0]^2+3/2*Pk[0, 0]*Pk[1, 0]+3/2*Pk[0, 0]*S1[0, 1]-3/2*Pk[0, 0]*S2[1, 0]+1/2*Pk[0, ...
S2[2, 0] = -5/2*Pk[1, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+5/2*Pk[0, 1]*S1[0, 0]+3*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*S2[0, 0]^2+3/2*Pk[0, 0]*Pk[1, 0]+3/2*Pk[0, 0]*S1[0, 1]-3/2*Pk[0, 0]*S2[1, 0]+1/2*Pk[0, ...
S2[2, 0] = -5/2*Pk[1, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+5/2*Pk[0, 1]*S1[0, 0]+3*Pk[0, 0]*S1[0, 0]^2+3*Pk[0, 0]*S2[0, 0]^2+3/2*Pk[0, 0]*Pk[1, 0]+3/2*Pk[0, 0]*S1[0, 1]-3/2*Pk[0, 0]*S2[1, 0]+1/2*Pk[0, ...
 

S2[1, 1] = 3*Pk[0, 1]*S2[0, 0]+6*Pk[0, 0]*S1[0, 0]*S2[0, 0]+3/2*Pk[0, 0]^2*S1[0, 0]+3/2*Pk[0, 0]*Pk[0, 1]+3*S2[0, 0]*S2[0, 1]+3/2*Pk[0, 0]*S2[0, 1]-2*S1[0, 0]*S1[0, 1]-Pk[1, 1]-2*Pk[1, 0]*S1[0, 0]-S1[...
S2[1, 1] = 3*Pk[0, 1]*S2[0, 0]+6*Pk[0, 0]*S1[0, 0]*S2[0, 0]+3/2*Pk[0, 0]^2*S1[0, 0]+3/2*Pk[0, 0]*Pk[0, 1]+3*S2[0, 0]*S2[0, 1]+3/2*Pk[0, 0]*S2[0, 1]-2*S1[0, 0]*S1[0, 1]-Pk[1, 1]-2*Pk[1, 0]*S1[0, 0]-S1[...
S2[1, 1] = 3*Pk[0, 1]*S2[0, 0]+6*Pk[0, 0]*S1[0, 0]*S2[0, 0]+3/2*Pk[0, 0]^2*S1[0, 0]+3/2*Pk[0, 0]*Pk[0, 1]+3*S2[0, 0]*S2[0, 1]+3/2*Pk[0, 0]*S2[0, 1]-2*S1[0, 0]*S1[0, 1]-Pk[1, 1]-2*Pk[1, 0]*S1[0, 0]-S1[...
S2[1, 1] = 3*Pk[0, 1]*S2[0, 0]+6*Pk[0, 0]*S1[0, 0]*S2[0, 0]+3/2*Pk[0, 0]^2*S1[0, 0]+3/2*Pk[0, 0]*Pk[0, 1]+3*S2[0, 0]*S2[0, 1]+3/2*Pk[0, 0]*S2[0, 1]-2*S1[0, 0]*S1[0, 1]-Pk[1, 1]-2*Pk[1, 0]*S1[0, 0]-S1[...
 

We can use an alternative choice of 4th order differential invariants 

fr > [u[3,0]=Pk, R1=u[4,0]+3*u[2,2],R2=u[3,1]+3*u[1,3], Pk*S1=3*u[3,1]+u[1,3], Pk*S2=3*u[2,2]+u[0,4], u[4,0]=Q0];
 

[u[3, 0] = Pk, R1 = u[4, 0]+3*u[2, 2], R2 = u[3, 1]+3*u[1, 3], Pk*S1 = 3*u[3, 1]+u[1, 3], Pk*S2 = 3*u[2, 2]+u[0, 4], u[4, 0] = Q0]
[u[3, 0] = Pk, R1 = u[4, 0]+3*u[2, 2], R2 = u[3, 1]+3*u[1, 3], Pk*S1 = 3*u[3, 1]+u[1, 3], Pk*S2 = 3*u[2, 2]+u[0, 4], u[4, 0] = Q0]
 

fr > C, Morphism :=dmf(LieAlg,Section, fr, [op(1..10,Section), u[3,0]=Pk, u[4, 0] = Q0, u[1, 3] = -1/8*Pk*S1+3/8*R2, u[2, 2] = 1/3*R1-1/3*Q0, u[3, 1] = -1/8*R2+3/8*Pk*S1, u[0, 4] = Pk*S2-R1+Q0], 'COM'):
 

"Section :" 

[u[0, 3] = 0, u[1, 2] = u[3, 0], u[2, 1] = 0, u[0, 2] = -1, u[1, 1] = 0, u[2, 0] = 1, u[0, 1] = 0, u[1, 0] = 0, u[0, 0] = 0, y = 0, x = 0]
[u[0, 3] = 0, u[1, 2] = u[3, 0], u[2, 1] = 0, u[0, 2] = -1, u[1, 1] = 0, u[2, 0] = 1, u[0, 1] = 0, u[1, 0] = 0, u[0, 0] = 0, y = 0, x = 0]
 

"Transversality condition :"   -u[3   0]
"Invariantizations:"
 

x = 0, y = 0, u[0, 0] = 0, u[1, 0] = 0, u[0, 1] = 0, u[1, 1] = 0, u[2, 0] = 1, u[0, 2] = -1, u[0, 3] = 0, u[2, 1] = 0, u[3, 0] = Pk, u[4, 0] = Q0, u[1, 3] = -1/8*Pk*S1+3/8*R2, u[2, 2] = 1/3*R1-1/3*Q0,...
x = 0, y = 0, u[0, 0] = 0, u[1, 0] = 0, u[0, 1] = 0, u[1, 1] = 0, u[2, 0] = 1, u[0, 2] = -1, u[0, 3] = 0, u[2, 1] = 0, u[3, 0] = Pk, u[4, 0] = Q0, u[1, 3] = -1/8*Pk*S1+3/8*R2, u[2, 2] = 1/3*R1-1/3*Q0,...
x = 0, y = 0, u[0, 0] = 0, u[1, 0] = 0, u[0, 1] = 0, u[1, 1] = 0, u[2, 0] = 1, u[0, 2] = -1, u[0, 3] = 0, u[2, 1] = 0, u[3, 0] = Pk, u[4, 0] = Q0, u[1, 3] = -1/8*Pk*S1+3/8*R2, u[2, 2] = 1/3*R1-1/3*Q0,...
x = 0, y = 0, u[0, 0] = 0, u[1, 0] = 0, u[0, 1] = 0, u[1, 1] = 0, u[2, 0] = 1, u[0, 2] = -1, u[0, 3] = 0, u[2, 1] = 0, u[3, 0] = Pk, u[4, 0] = Q0, u[1, 3] = -1/8*Pk*S1+3/8*R2, u[2, 2] = 1/3*R1-1/3*Q0,...
 

"Commutation rules" 

[D1, D2] = -1/12*(3*u[3, 1]+u[1, 3])*D1/u[3, 0]-1/12*(-u[0, 4]-3*u[2, 2])*D2/u[3, 0] 

"Syzygies" 

[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
[-R2[0, 0]*S2[0, 0]+S1[0, 0]*R1[0, 0]+12*R2[1, 0]-12*R1[0, 1], 8*S2[0, 0]*R1[0, 0]+48*Pk[0, 0]^2*S2[0, 0]-96*Pk[0, 0]*R1[0, 0]+5*R2[0, 0]*S1[0, 0]-4*Pk[0, 0]*S1[0, 0]^2+96*Pk[0, 0]*Q0[0, 0]-12*R2[0, 1...
 

"Syzygies with attempt to eliminate the y's" 

y0[0, 0] = -1/Pk[0, 0] 

y1[0, 0] = Pk[0, 0] 

R2[1, 0] = 1/12*R2[0, 0]*S2[0, 0]-1/12*S1[0, 0]*R1[0, 0]+R1[0, 1] 

Q0[1, 0] = -1/4*S2[0, 0]*R1[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+3*Pk[0, 0]*R1[0, 0]-5/32*R2[0, 0]*S1[0, 0]+1/8*Pk[0, 0]*S1[0, 0]^2-3*Pk[0, 0]*Q0[0, 0]+3/8*R2[0, 1]-9/8*Pk[0, 0]*S1[0, 1]+R1[1, 0]
Q0[1, 0] = -1/4*S2[0, 0]*R1[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+3*Pk[0, 0]*R1[0, 0]-5/32*R2[0, 0]*S1[0, 0]+1/8*Pk[0, 0]*S1[0, 0]^2-3*Pk[0, 0]*Q0[0, 0]+3/8*R2[0, 1]-9/8*Pk[0, 0]*S1[0, 1]+R1[1, 0]
Q0[1, 0] = -1/4*S2[0, 0]*R1[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+3*Pk[0, 0]*R1[0, 0]-5/32*R2[0, 0]*S1[0, 0]+1/8*Pk[0, 0]*S1[0, 0]^2-3*Pk[0, 0]*Q0[0, 0]+3/8*R2[0, 1]-9/8*Pk[0, 0]*S1[0, 1]+R1[1, 0]
 

Pk[1, 0] = 1/4*R1[0, 0] 

S1[1, 0] = 1/12*(-R2[0, 0]*S2[0, 0]+4*Pk[0, 0]*S1[0, 0]*S2[0, 0]+12*Pk[0, 0]^2*S1[0, 0]-36*Pk[0, 0]*R2[0, 0]+4*R1[0, 1]+32*Q0[0, 1]-6*S1[0, 0]*R1[0, 0])/Pk[0, 0]
S1[1, 0] = 1/12*(-R2[0, 0]*S2[0, 0]+4*Pk[0, 0]*S1[0, 0]*S2[0, 0]+12*Pk[0, 0]^2*S1[0, 0]-36*Pk[0, 0]*R2[0, 0]+4*R1[0, 1]+32*Q0[0, 1]-6*S1[0, 0]*R1[0, 0])/Pk[0, 0]
S1[1, 0] = 1/12*(-R2[0, 0]*S2[0, 0]+4*Pk[0, 0]*S1[0, 0]*S2[0, 0]+12*Pk[0, 0]^2*S1[0, 0]-36*Pk[0, 0]*R2[0, 0]+4*R1[0, 1]+32*Q0[0, 1]-6*S1[0, 0]*R1[0, 0])/Pk[0, 0]
 

S2[1, 0] = 3/2*Pk[0, 0]*S2[0, 0]-5/2*R1[0, 0]+4*Q0[0, 0]-1/12*S1[0, 0]^2+1/12*S2[0, 0]^2+S1[0, 1]
S2[1, 0] = 3/2*Pk[0, 0]*S2[0, 0]-5/2*R1[0, 0]+4*Q0[0, 0]-1/12*S1[0, 0]^2+1/12*S2[0, 0]^2+S1[0, 1]
 

Pk[0, 1] = 1/4*R2[0, 0] 

and obtain the following rewriting: 

fr > COM;
 

[[y, x] = [1/12*S1[0, 0], -1/12*S2[0, 0]]] 

fr > R := differential_ring(ranking=[[y0,y1],[Q0,R1,R2],[S1,S2,Pk]], derivations=[x,y],
          commutations=COM,notation=vjet);
CC := Rosenfeld_Groebner(equations(C[1]), [Pk[0,0]], R);
for g in rewrite_rules(CC[1]) do print(g); od;
 

(Typesetting:-mprintslash)([R := differential_polynomial_ring], [differential_polynomial_ring]) 

(Typesetting:-mprintslash)([CC := [characterisable]], [[characterisable]]) 

y0[0, 0] = -1/Pk[0, 0] 

y1[0, 0] = Pk[0, 0] 

Q0[0, 0] = 1/4*S2[1, 0]-3/8*Pk[0, 0]*S2[0, 0]+5/2*Pk[1, 0]+1/48*S1[0, 0]^2-1/48*S2[0, 0]^2-1/4*S1[0, 1]
Q0[0, 0] = 1/4*S2[1, 0]-3/8*Pk[0, 0]*S2[0, 0]+5/2*Pk[1, 0]+1/48*S1[0, 0]^2-1/48*S2[0, 0]^2-1/4*S1[0, 1]
 

R1[0, 0] = 4*Pk[1, 0] 

R2[0, 0] = 4*Pk[0, 1] 

Pk[2, 1] = 3*Pk[0, 0]^2*Pk[0, 1]+2/9*Pk[1, 0]*S2[0, 1]+1/3*Pk[0, 0]^2*S2[0, 1]+1/6*Pk[0, 2]*S1[0, 0]+2/9*Pk[0, 1]*S1[0, 1]-11/432*Pk[0, 1]*S1[0, 0]^2+1/4*Pk[0, 0]*S1[0, 2]+1/12*S1[0, 0]*Pk[2, 0]+5/12*...
Pk[2, 1] = 3*Pk[0, 0]^2*Pk[0, 1]+2/9*Pk[1, 0]*S2[0, 1]+1/3*Pk[0, 0]^2*S2[0, 1]+1/6*Pk[0, 2]*S1[0, 0]+2/9*Pk[0, 1]*S1[0, 1]-11/432*Pk[0, 1]*S1[0, 0]^2+1/4*Pk[0, 0]*S1[0, 2]+1/12*S1[0, 0]*Pk[2, 0]+5/12*...
Pk[2, 1] = 3*Pk[0, 0]^2*Pk[0, 1]+2/9*Pk[1, 0]*S2[0, 1]+1/3*Pk[0, 0]^2*S2[0, 1]+1/6*Pk[0, 2]*S1[0, 0]+2/9*Pk[0, 1]*S1[0, 1]-11/432*Pk[0, 1]*S1[0, 0]^2+1/4*Pk[0, 0]*S1[0, 2]+1/12*S1[0, 0]*Pk[2, 0]+5/12*...
Pk[2, 1] = 3*Pk[0, 0]^2*Pk[0, 1]+2/9*Pk[1, 0]*S2[0, 1]+1/3*Pk[0, 0]^2*S2[0, 1]+1/6*Pk[0, 2]*S1[0, 0]+2/9*Pk[0, 1]*S1[0, 1]-11/432*Pk[0, 1]*S1[0, 0]^2+1/4*Pk[0, 0]*S1[0, 2]+1/12*S1[0, 0]*Pk[2, 0]+5/12*...
Pk[2, 1] = 3*Pk[0, 0]^2*Pk[0, 1]+2/9*Pk[1, 0]*S2[0, 1]+1/3*Pk[0, 0]^2*S2[0, 1]+1/6*Pk[0, 2]*S1[0, 0]+2/9*Pk[0, 1]*S1[0, 1]-11/432*Pk[0, 1]*S1[0, 0]^2+1/4*Pk[0, 0]*S1[0, 2]+1/12*S1[0, 0]*Pk[2, 0]+5/12*...
Pk[2, 1] = 3*Pk[0, 0]^2*Pk[0, 1]+2/9*Pk[1, 0]*S2[0, 1]+1/3*Pk[0, 0]^2*S2[0, 1]+1/6*Pk[0, 2]*S1[0, 0]+2/9*Pk[0, 1]*S1[0, 1]-11/432*Pk[0, 1]*S1[0, 0]^2+1/4*Pk[0, 0]*S1[0, 2]+1/12*S1[0, 0]*Pk[2, 0]+5/12*...
Pk[2, 1] = 3*Pk[0, 0]^2*Pk[0, 1]+2/9*Pk[1, 0]*S2[0, 1]+1/3*Pk[0, 0]^2*S2[0, 1]+1/6*Pk[0, 2]*S1[0, 0]+2/9*Pk[0, 1]*S1[0, 1]-11/432*Pk[0, 1]*S1[0, 0]^2+1/4*Pk[0, 0]*S1[0, 2]+1/12*S1[0, 0]*Pk[2, 0]+5/12*...
Pk[2, 1] = 3*Pk[0, 0]^2*Pk[0, 1]+2/9*Pk[1, 0]*S2[0, 1]+1/3*Pk[0, 0]^2*S2[0, 1]+1/6*Pk[0, 2]*S1[0, 0]+2/9*Pk[0, 1]*S1[0, 1]-11/432*Pk[0, 1]*S1[0, 0]^2+1/4*Pk[0, 0]*S1[0, 2]+1/12*S1[0, 0]*Pk[2, 0]+5/12*...
 

S2[2, 0] = -5/2*Pk[1, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+18*Pk[0, 0]*Pk[1, 0]-5/2*Pk[0, 1]*S1[0, 0]+1/4*Pk[0, 0]*S1[0, 0]^2-3/2*Pk[0, 0]*S2[1, 0]+1/4*Pk[0, 0]*S2[0, 0]^2-3/2*Pk[0, 0]*S1[0, 1]+6*Pk[0,...
S2[2, 0] = -5/2*Pk[1, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+18*Pk[0, 0]*Pk[1, 0]-5/2*Pk[0, 1]*S1[0, 0]+1/4*Pk[0, 0]*S1[0, 0]^2-3/2*Pk[0, 0]*S2[1, 0]+1/4*Pk[0, 0]*S2[0, 0]^2-3/2*Pk[0, 0]*S1[0, 1]+6*Pk[0,...
S2[2, 0] = -5/2*Pk[1, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+18*Pk[0, 0]*Pk[1, 0]-5/2*Pk[0, 1]*S1[0, 0]+1/4*Pk[0, 0]*S1[0, 0]^2-3/2*Pk[0, 0]*S2[1, 0]+1/4*Pk[0, 0]*S2[0, 0]^2-3/2*Pk[0, 0]*S1[0, 1]+6*Pk[0,...
S2[2, 0] = -5/2*Pk[1, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S2[0, 0]+18*Pk[0, 0]*Pk[1, 0]-5/2*Pk[0, 1]*S1[0, 0]+1/4*Pk[0, 0]*S1[0, 0]^2-3/2*Pk[0, 0]*S2[1, 0]+1/4*Pk[0, 0]*S2[0, 0]^2-3/2*Pk[0, 0]*S1[0, 1]+6*Pk[0,...
 

S2[1, 1] = 3*Pk[0, 1]*S2[0, 0]-1/2*Pk[0, 0]*S1[0, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S1[0, 0]+18*Pk[0, 0]*Pk[0, 1]-12*Pk[1, 1]+2*Pk[1, 0]*S1[0, 0]-1/6*S1[0, 0]*S1[0, 1]-1/12*S1[0, 0]*S2[1, 0]+1/4*S2[0, 0]*S2[...
S2[1, 1] = 3*Pk[0, 1]*S2[0, 0]-1/2*Pk[0, 0]*S1[0, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S1[0, 0]+18*Pk[0, 0]*Pk[0, 1]-12*Pk[1, 1]+2*Pk[1, 0]*S1[0, 0]-1/6*S1[0, 0]*S1[0, 1]-1/12*S1[0, 0]*S2[1, 0]+1/4*S2[0, 0]*S2[...
S2[1, 1] = 3*Pk[0, 1]*S2[0, 0]-1/2*Pk[0, 0]*S1[0, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S1[0, 0]+18*Pk[0, 0]*Pk[0, 1]-12*Pk[1, 1]+2*Pk[1, 0]*S1[0, 0]-1/6*S1[0, 0]*S1[0, 1]-1/12*S1[0, 0]*S2[1, 0]+1/4*S2[0, 0]*S2[...
S2[1, 1] = 3*Pk[0, 1]*S2[0, 0]-1/2*Pk[0, 0]*S1[0, 0]*S2[0, 0]-3/2*Pk[0, 0]^2*S1[0, 0]+18*Pk[0, 0]*Pk[0, 1]-12*Pk[1, 1]+2*Pk[1, 0]*S1[0, 0]-1/6*S1[0, 0]*S1[0, 1]-1/12*S1[0, 0]*S2[1, 0]+1/4*S2[0, 0]*S2[...