level 1
junolee8
楼主
想用mathematica求解两个自变量的字符方程组, 有一些复杂,运行后几个小时一直都不出结果,想问一下哪里出了问题。。。下面是代码 自变量是de和phie
Solve[{de ^(-alpha - 1)*(1 +
a exp (-b (Pi/2 - phir - a)))/(grr dr^(-alpha) rr A hr (1 +
a exp (-b (Pi/2 - phie - a)))) == (dr cos (phir - phie) -
de)*ur^2*(de^2 + dr^2 - 2 de dr cos (phir - phie))^(alpha/2 +
1)/(ur^2*(de^2 + dr^2 - 2*de*dr*cos (phir - phie))^(alpha/2) +
rr*grt*ge*A)^2,
de == -alpha/a/b*dr*
sin (phir - phie)*(exp (b*(Pi/2 - phie - a)) + a) +
dr*cos (phir - phie)}, {de, phie}]
2018年06月06日 06点06分
1
Solve[{de ^(-alpha - 1)*(1 +
a exp (-b (Pi/2 - phir - a)))/(grr dr^(-alpha) rr A hr (1 +
a exp (-b (Pi/2 - phie - a)))) == (dr cos (phir - phie) -
de)*ur^2*(de^2 + dr^2 - 2 de dr cos (phir - phie))^(alpha/2 +
1)/(ur^2*(de^2 + dr^2 - 2*de*dr*cos (phir - phie))^(alpha/2) +
rr*grt*ge*A)^2,
de == -alpha/a/b*dr*
sin (phir - phie)*(exp (b*(Pi/2 - phie - a)) + a) +
dr*cos (phir - phie)}, {de, phie}]