V=newvect(17,[1,a1,a2,a3,a4,a5,a6,a7,a8,a9,a10,a11])$

S=V[1]^2-4*V[2]$

for(N=3;N<12;N++){

  T=V[N-1]^2*S$
  S=T$

  for(K=1;K<N;K++){
    U=(-N)*diff(T,V[1])$
    for(J=2;J<N;J++){
      U=U+((J-(N+1))*V[J-1])*diff(T,V[J])$
    }  
    T=sdiv(U,K*V[N-1])$
    S=S+T*V[N]^K$
  }

  print(nmono(S))$
}

end$