Symmetry group for orthogonal wave equations 

Study case suggested by G. Metivier, Universite de Rennes. 

> coord_frame([x,y],[s,phi,psi],fr):
LieAlg:= [ vect(x),vect(y),
      v_zip([-y,x],[x,y],plus, vect),            
      v_zip([x,y,2*s[0,0]],[x,y,s[0,0]],plus,vect),
      v_zip([phi[0,0]],[phi[0,0]],plus,vect),
      v_zip([psi[0,0]],[psi[0,0]],plus,vect)]:
 

> Section:=[x=0,y=0,s[0,0]=1,phi[0,0]=1, psi[0,0]=1, s[1,0]=0]: stair(Section);
 

Plot 

Plot 

Plot 

> C, Invariantizations :=dmf(LieAlg,Section, fr, [op(Section), s[0,1]=2*s1, s[2,0] = 2*s1*(s3+s1),s[1,1] = 2*s1*s2, phi[1,0]=phi1, phi[0,1]=phi2, psi[1,0]=psi1, psi[0,1]=psi2],'COM'):
 

"Section :" 

[s[1, 0] = 0, psi[0, 0] = 1, phi[0, 0] = 1, s[0, 0] = 1, y = 0, x = 0] 

"Transversality condition :"   -s[0   1]
"Invariantizations:"
 

x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1, s[2, 0] = 2*s1*(s3+s1), s[1, 1] = 2*s1*s2, phi[1, 0] = phi1, phi[0, 1] = phi2, psi[1, 0] = psi1, psi[0, 1] = psi2,...
x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1, s[2, 0] = 2*s1*(s3+s1), s[1, 1] = 2*s1*s2, phi[1, 0] = phi1, phi[0, 1] = phi2, psi[1, 0] = psi1, psi[0, 1] = psi2,...
x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1, s[2, 0] = 2*s1*(s3+s1), s[1, 1] = 2*s1*s2, phi[1, 0] = phi1, phi[0, 1] = phi2, psi[1, 0] = psi1, psi[0, 1] = psi2,...
x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1, s[2, 0] = 2*s1*(s3+s1), s[1, 1] = 2*s1*s2, phi[1, 0] = phi1, phi[0, 1] = phi2, psi[1, 0] = psi1, psi[0, 1] = psi2,...
 

"Commutation rules" 

[D1, D2] = -1/2*(-2*s[2, 0]+s[0, 1]^2)*D1/s[0, 1]+s[1, 1]*D2/s[0, 1] 

"Syzygies" 

[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
[s2[0, 0]*y5[0, 0]-s2[0, 0]*y7[0, 0]-psi2[0, 0]*y6[0, 0]-2*y6[0, 0]*s3[0, 0]+psi1[0, 0]*y7[0, 0]-y6[0, 1]+y7[1, 0], -2*s2[0, 0]*y6[0, 0]+y5[0, 0]*s1[0, 0]-psi2[0, 0]*y5[0, 0]-y5[0, 0]*s3[0, 0]+y7[0, 0...
 

"Syzygies with attempt to eliminate the y's" 

y0[0, 0] = -1/2/s1[0, 0] 

y7[0, 0] = psi2[0, 1]-s2[0, 0]*psi1[0, 0]-psi2[0, 0]*s1[0, 0]+psi2[0, 0]^2 

y6[0, 0] = psi1[0, 1]+s2[0, 0]*psi2[0, 0]-s1[0, 0]*psi1[0, 0]+psi2[0, 0]*psi1[0, 0] 

y5[0, 0] = psi1[1, 0]+psi2[0, 0]*s3[0, 0]+psi2[0, 0]*s1[0, 0]+psi1[0, 0]^2 

y4[0, 0] = phi2[0, 1]-s2[0, 0]*phi1[0, 0]-phi2[0, 0]*s1[0, 0]+phi2[0, 0]^2 

y3[0, 0] = phi1[0, 1]+s2[0, 0]*phi2[0, 0]-s1[0, 0]*phi1[0, 0]+phi2[0, 0]*phi1[0, 0] 

y2[0, 0] = phi1[1, 0]+phi2[0, 0]*s3[0, 0]+phi2[0, 0]*s1[0, 0]+phi1[0, 0]^2 

y1[0, 0] = 2*s1[0, 1]+2*s1[0, 0]^2 

s2[1, 0] = s1[0, 1]+s2[0, 0]^2+s3[0, 0]^2+s1[0, 0]*s3[0, 0]+s3[0, 1] 

s1[1, 0] = s2[0, 0]*s1[0, 0] 

psi2[1, 0] = psi1[0, 1]+s2[0, 0]*psi2[0, 0]+psi1[0, 0]*s3[0, 0] 

phi2[1, 0] = phi1[0, 1]+s2[0, 0]*phi2[0, 0]+phi1[0, 0]*s3[0, 0] 

fr > COM;
 

[[y, x] = [-s3[0, 0], -s2[0, 0]]] 

fr > Invariantizations;
 

[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
[x = 0, y = 0, s[0, 0] = 1, phi[0, 0] = 1, psi[0, 0] = 1, s[1, 0] = 0, s[0, 1] = 2*s1[0, 0], s[2, 0] = 2*s1[0, 0]*(s3[0, 0]+s1[0, 0]), s[1, 1] = 2*s2[0, 0]*s1[0, 0], phi[1, 0] = phi1[0, 0], phi[0, 1] ...
 

fr >