How to use fsolve to solve a set of nonlinear equations in two variables?

I am attempting to use fsolve to solve a set of equations in X & Y as given below:

Here, the terms B33,B35,B55,C1 and C2 are also functions of X. I tried to use fsolve as explained in the documentation here. How do I make sure that fsolve uses the modified value of B33,B35,B55,C1 and C2 during each iteration? I have attached the .m files. Any help would be much appreciated.

[Matt J]

I see nothing wrong with the way you've done it.

[Sarath S]

Thank you! I was a bit confused since the result I was getting did not match with some previously verified values. I believe the mistake must be in the actual expression itself. Must have typed out a wrong factor or something. Anyways, thank you again.