Thank you for the reply. I realized my mistake, actually it should be val1 = gn /. solvegn1[[1]] for reading first root out of four roots.
Solving first equation for each n gives me four roots of gn and I was trying to choose only real positive ones for each n (let's say n=0 and four gn and trying to choose for only real positive values) so I used Real//Positive.
I will forget nature of the roots and If I take into account all four roots then plugging those in last two equations gives me four y and four z values but journal calculation shows only one positive y and z values so I was trying to choose only real positive gn value.
Could you please solve for y and z values for each gn?
In short I am trying to solve for y and z values from last two equations.
Thank you
Attachments: