Computation of dparz(z)

This section is basically the same as the previous one in execution.

> dparz1:=collect(subs(d(u)=0,dg),p);

The next three steps are a bit of 'sleight of hand' to get Maple to remove the appropriate terms... namely those with a p in the numerator, and no p's in the denominator.

> dparzromeo:=simplify(p*dparz1);

> dparzlist:=convert(dparzromeo,list);

> dparzjuliet:=remove(has,dparzlist,{p^2,p^3,p^4});

> dparz2:=convert(dparzjuliet,`+`)/p;

Now we simplify dparz2, set it equal to zero, and solve for d(z). This is the final step which gives us dparz(z).

> dparzz:=solve(dparz2=0,d(z));

>