用MATLAB怎么求这个带参数的二元三次方程组的解？（解决问题追加高分）

Question

前面发了个帖子，但是问了多个问题，不好回答。现在把最关键的问题提出来。不知道这样是不是违反了版规，求版主高抬贵手。

我要求一个带参数的二元三次方程。二元三次方程是对一个函数的两个变量分别求偏导构成的。这个函数由以下两个函数构成
Obj=-116*100*y/(2*pi*x*(2.2+0.02+0.02*x^2+5/x^2+(50000/(10*x^2))*(y/(2*pi))^2)) %这个是第一个函数，实际是从一个M文件中出来的
function [c,ceq]=Constr(x,y) %这个是第二个函数，也是一个M文件
syms x y;
c(1)=x/y-1500/(8*100/60);
c(2)=0.2/(2*100/60)-x/y;
c(3)=2.2/(2.2+0.02+0.02*x^2+5/x^2+50000/(10*x^2)*(y/(2*pi))^2)-0.8;
c(4)=0.5-2.2/(2.2+0.02+0.02*x^2+5/x^2+50000/(10*x^2)*(y/(2*pi))^2);
c(5)=2.2+0.02+0.02*x^2+5/x^2+5000/(10*x^2)*(y/(2*pi))^2-9.6*0.08*116*100/1500;
ceq=[];
f=-Obj(x,y)-r0*sum([Constr].\[1 1 1 1 1]); %这个是第三个函数，由上面两个函数复合而成
二元三次方程
dpx=diff(f,'x')=0
dpy=diff(f,'y')=0
这个方程里实际有三个变量 x y r0， 我希望求出x和y 的带有r0的表达式。
代码怎么写啊 （如果要用初值，x=2 y=1）
求大神帮助！！！！！！！！！！！

Answer 1

很简单，首先要把第一个函数中的x转变成左右等式的形式我们把它记作1式，再把1式与第二个函数结合依然写成于1式相同的形式，我们把这个形式记作2式，最后一个函数是最长也是最复杂的一个函式，我们把该式进行左右归类的简单划分，再把2式代入，用复合函数等量代换，等式左右恒等的性质可解出，以上为解题思路，具体步骤不细说了，麻烦楼主自己去解吧。