[转载] 请maple 或者 mathematica 的高手帮个忙！# Computation - 科学计算
l*t
1 楼
【 以下文字转载自 Physics 讨论区,原文如下 】
发信人: leot (leot), 信区: Physics
标 题: 请maple 或者 mathematica 的高手帮个忙！
发信站: Unknown Space - 未名空间 (Mon Apr 26 15:20:42 2004) WWW-POST
最近看到一片paper,上面给出一个微分方程（以下所列方程均为maple编辑形式）：
(beta[1]/(D+lambda[1])+beta[2]/(D+lambda[2]))*diff(n(D),D)=(beta[1]/(D+lambda[
1])^2+beta[2]/(D+lambda[2])^2-1/tau-alpha)*n(D)+1/tau*m(D);
initial condition is, n(D[0])=m(D[0]);
作者给了一个非常复杂的解析解（m(D)应该是一个已知函数）：
n(D)=(1-alpha*tau)*theta[4]*exp(-theta[1])*(theta[2]^theta[3])*m(D[0])+theta[5
]*exp
(-thet
发信人: leot (leot), 信区: Physics
标 题: 请maple 或者 mathematica 的高手帮个忙！
发信站: Unknown Space - 未名空间 (Mon Apr 26 15:20:42 2004) WWW-POST
最近看到一片paper,上面给出一个微分方程（以下所列方程均为maple编辑形式）：
(beta[1]/(D+lambda[1])+beta[2]/(D+lambda[2]))*diff(n(D),D)=(beta[1]/(D+lambda[
1])^2+beta[2]/(D+lambda[2])^2-1/tau-alpha)*n(D)+1/tau*m(D);
initial condition is, n(D[0])=m(D[0]);
作者给了一个非常复杂的解析解（m(D)应该是一个已知函数）：
n(D)=(1-alpha*tau)*theta[4]*exp(-theta[1])*(theta[2]^theta[3])*m(D[0])+theta[5
]*exp
(-thet