求助二元四阶微分方程求解
本帖最后由 TBE_Legend 于 2009-12-11 09:38 编辑二元四阶微分方程如附件,怎么结果中有积分和复数啊?
后面是代码
omega := f*sin(Pi*m*x/a)*(sin(Pi*y/b))^2;
eq1 := E*((diff(diff(omega, y), x))^2-(diff(diff(omega, x), x))*(diff(diff(omega, y), y)));
eq2:=combine(eq1);
eq3 := diff(diff(diff(diff(F(x, y), x), x), x), x)+2*(diff(diff(diff(diff(F(x, y), y), y), x), x))+diff(diff(diff(diff(F(x, y), y), y), y), y);
pde := eq2 = eq3;
sol := pdsolve(pde,F(x,y)); 最好用代码或者图片说明意图. 第二张图太小,不能说明问题 截图真的很小,无法看清. 顶起来,希望高手帮忙,谢谢 不是所有的pde都有解的,4阶2元出现复数很正常,君不见2元ODE在delta<0的情况。 那后面那部分怎么是积分型式呢? 用build选项尝试获取显式. 各位,我已经把代码上传,希望大家帮帮忙 9# bpzf
pdetest(sol, pde)
解没有问题. 哦,谢谢maplelab 了,我再看看 本帖最后由 bpzf 于 2010-1-20 11:29 编辑
为什么用显式的就计算不出来了,结果用pdetest检验不等于0???
omega := f*sin(Pi*m*x/a)*(sin(Pi*y/b))^2;
eq1 := E*((diff(diff(omega, y), x))^2-(diff(diff(omega, x), x))*(diff(diff(omega, y), y)));
eq2:=combine(eq1);
eq3 := diff(diff(diff(diff(F(x, y), x), x), x), x)+2*(diff(diff(diff(diff(F(x, y), y), y), x), x))+diff(diff(diff(diff(F(x, y), y), y), y), y);
pde := eq2 = eq3;
sol := pdsolve(pde, HINT = `+`, INTEGRATE, build);
页:
[1]