diffop(E*X,[X,E],[1,X*E],10)
diffop(Mod(S/C,S^2+C^2-1),[C],[-S],10)
E4(n)=1+240*sum(i=1,n,sigma(i,3)*q^i)+O(q^(n+1));
diffop(E4(15),[q],[q])
