求解非线性方程组（含参数）

编写如下代码可以求解两个方程的数值解

 [x,y]=solve('-24279+(8.372+8.314*(ln(1-y)-ln(1-x)))*2850=0','-16660+(11.72+8.314*(ln(y)-ln(x)))*2850=0');
>> x=vpa(x,4)
y=vpa(y,4)

现在想要将2850换为变参数t（2000-3000的整数），想得到x=f(t); y=f(t)的表达式，用solve函数的话无法求解，在论坛里面发现有人用fsolve求解，请问该如何解决？