【求助】解偏微分方程

jelin

A,B,C,k,h,m均为非零常数，x0>0.

新手，不知如何编写？望大侠救助。

TBE_Legend

http://bb.3861.com.cn/attachment/200911/30/59857_1259572831m44F.jpg
A,B,C,k,h,m均为非零常数，x0>0.

新手，不知如何编写？望大侠救助。
jelin 发表于 2009-12-2 09:12

解不出来。

1. 
   
2. Clear["Global`*"]
   
3. a = b = c = k = h = m = 1;
   
4. NDSolve[{x^2 D[y[x, t], {x, 2}] + 2 x D[y[x, t], x] + a x^2 y[x, t] +
   
5.     b x^2 D[y[x, t], t] == c/x^2, y[x, 0] == 1, y[10^50, t] == 10^-5,
   
6.   Derivative[1, 0][y][1, t] == y[x, t] + m},
   
7. y[x, t], {x, 1, 2}, {t, 0, 1}]
   

复制代码

不带入任何初边值条件，maple能帮你分解成ode，试试这个ode吧。
restart;
pde:=x^2*diff(y(x,t),x,x)+2*x*diff(y(x,t),x)+A*x^2*y(x,t)+B*x^2*diff(y(x,t),t)=C/x^2;
simplify(pdsolve(pde),size,build);

1. restart;
   
2. pde:=x^2*diff(y(x,t),x,x)+2*x*diff(y(x,t),x)+A*x^2*y(x,t)+B*x^2*diff(y(x,t),t)=C/x^2;
   
3. simplify(pdsolve(pde),size,build);
   

复制代码

jelin

太感谢了，向你学习了。
只有寻求其他方式解了， 想试试Green函数法。