Greetings Liyaqat,
Before we work on solving this problem, have you checked that the answer provided is correct?
I tried to verify it and assuming everything else is correct, I find that
Vcc == 0 && -1. C L \[Omega] + C^2 L^2 \[Omega]^3 != 0
This means that Vcc must be zero when you say it is assumed to be constant. Is this an expected result?
I attached some the edits I made to produce this result.
qsol = 1.607;
phisol = -41.614;
VRsol = 0.9253*Vcc;
eqn1 = eq4 /. {q -> qsol, phi -> phisol, VR -> VRsol} // FullSimplify
eqn2 = eq5b /. {q -> qsol, phi -> phisol, VR -> VRsol} //
FullSimplify
eqn3 = eq6c /. {q -> qsol, phi -> phisol, VR -> VRsol} // FullSimplify
Reduce[{eqn1, eqn2, eqn3}]
Attachments: