# 
usk:=seq(u||i, i=1..genus):
# 
SS2 := - proc (i, j) 
         factor(diff(log(si(usk)), u || i, u || j));
         end proc:
SS3 := - proc (i, j, k) 
         factor(diff(log(si(usk)), u || i, u || j, u || k));
         end proc:
SS4 := - proc (i, j, k, l)
         factor(diff(log(si(usk)), u || i, u || j, u || k, u || l));
         end proc:
SS5 := - proc (i, j, k, l, m)
         factor(diff(log(si(usk)), u || i, u || j, u || k, u || l, u || m));
         end proc:

p2subs := seq(seq(P(i,j)=SS2(i,j),j=i..genus),i=1..genus):
p3subs := seq(seq(seq(P(i,j,k)=SS3(i,j,k),k=j..genus),j=i..genus),i=1..genus):
p4subs := seq(seq(seq(seq(P(i,j,k,l)=SS4(i,j,k,l),l=k..genus),k=j..genus), j=i..genus),i=1..genus):
p5subs := seq(seq(seq(seq(seq(P(i,j,k,l,m)=SS5(i,j,k,l,m),m=l..genus),l=k..genus),k=j..genus), j=i..genus),i=1..genus):
psubs  := p5subs,p4subs,p3subs,p2subs: