[求助]VC 调用 fsolve 函数

先表谢意，小弟问题如下，还望高手给与解决.
下面为matlab 中的 .m文件
function Q=fun(T)
x=T(1);y=T(2);Q=zeros(2,1);
a=0.5; //此处数值0.5为变量，需在VC中调用变量P[i][j]。
P=a*a*a*(10-15*a+6*a*a);
D=1050; //此处数值1050为变量，需在VC中调用变量温度 T 。
C=0.02; //此处数值0.02为变量，需在VC中调用变量C[i][j] 。
Q(1)=C-[1-P]*x-P*y;
Q(2)=(322050+75.677*D)/3+8.314*D*log(y/(3-y))+(1/3-2*y/9)*(-190*D)-77207+15.877*D-8.314*D*log(x/(1-x))-(1-2*x)*(-34671);

然后在matlab中调用函数 X=fsolve(@fun,[0.02 0.01],optimset('Display','off')) 所得到的数值 需要赋值给VC中的Ur[i][j]和Ua[i][j].
现在小弟想在VC中实现此函数的中间计算，希望大侠能够给点帮助！谢谢了！