数学建模社区-数学中国
标题:
Maple解偏微分方程,解的化简问题 谢谢
[打印本页]
作者:
一汪清水
时间:
2013-11-29 19:52
标题:
Maple解偏微分方程,解的化简问题 谢谢
PDE := x(1-x)*(diff(z(x, y), x))/(1-bx)+y(1-y)*(diff(z(x, y), y))/(1-bx)+(1+ay/(1-bx))*z(x, y)-ay*(1-y+ay+b(b-1)*y^2/(1-bx))/(1-bx)^3 = 0
* d1 B. F* }& q9 Z. O' x. {
% @' d2 x; B l4 U' H; ]
pdsolve(PDE, z(x, y))
* J- S2 G# W6 O( ? \
7 R" }3 A3 }) A$ c" _* D
z(x, y) = exp(Int(-(1-bx+ay)/x(1-x), x))*Intat(-ay*exp((1-bx+ay)*(Int(1/_f(1-_f), _f)))*(ay-ay*bx+b(b-1)*RootOf(Int(1/_f(1-_f), _f)-Intat(1/_a(1-_a), _a = _Z)-(Int(1/x(1-x), x))+Int(1/y(1-y), y))^2-RootOf(Int(1/_f(1-_f), _f)-Intat(1/_a(1-_a), _a = _Z)-(Int(1/x(1-x), x))+Int(1/y(1-y), y))+RootOf(Int(1/_f(1-_f), _f)-Intat(1/_a(1-_a), _a = _Z)-(Int(1/x(1-x), x))+Int(1/y(1-y), y))*bx+1-bx)/(_f(1-_f)*(-1+bx)^3), _f = x)+exp(Int(-(1-bx+ay)/x(1-x), x))*_F1(-(Int(1/x(1-x), x))+Int(1/y(1-y), y))
. x# ]) `0 y$ Q6 I4 k
' f' S' A4 s' N% N3 i. Q8 X
这个解太长了,含有很多积分,而且还有好几个rootof
; |5 ~$ }" s, ]) ^' q9 ~3 z- ^# N
能否将解进行一下简化呢?该如何简化,求指点
$ B4 H0 ~4 B! ~" I# M: i) e8 v$ E
作者:
madio
时间:
2014-3-6 22:31
使用remove_RootOf试试看
欢迎光临 数学建模社区-数学中国 (http://www.madio.net/)
Powered by Discuz! X2.5