A variant of Milan's approach: sol: solve( Gcl=Gact*Gth/(Gact*Gth+1) , Gact); -- or solve for Gth factor(subst(sol,myvar)) => Gcl*(Ga*s^2*xbeam+thdhat) -s