| 读帖时,帖子不存在 |
|
莫回避苹果方程——我的方程可以是2个半苹果相加=1个苹果,这才是奇迹。
曾云海 |
|
>> syms x y r
>> [x,y]=dsolve('Dx+2*x/r=Dy/2','3*Dy+2*y*r^2=2*Dx','r') %f=r^2的解 x = exp(-1/6*r^3)*(C1*WhittakerM(1/6,1,1/3*r^3)+C2*WhittakerW(1/6,1,1/3*r^3))/r^(3/2) y = 2/3*exp(-1/6*r^3)*(3*C1*WhittakerM(1/6,1,1/3*r^3)*r^3-6*C1*WhittakerM(1/6,1,1/3*r^3)-12*C1*WhittakerM(-5/6,1,1/3*r^3)+3*C2*WhittakerW(1/6,1,1/3*r^3)*r^3-6*C2*WhittakerW(1/6,1,1/3*r^3)+8*C2*WhittakerW(-5/6,1,1/3*r^3))/r^(9/2) 一堆的惠泰克函数 >> [x,y]=dsolve('Dx+2*x/r=Dy/2','3*Dy+2*y*(r^3+r)=2*Dx','r') %f=r^3+r x = DESol({(r^4+r^2)/r^2/(r^2+1)*diff(_Y(r),$(r,2))+(r^7+2*r^5+r^3+2*r)/r^2/(r^2+1)*diff(_Y(r),r)+(2*r^6+4*r^4-10*r^2-6)/r^2/(r^2+1)*_Y(r)},{_Y(r)}) y = -2*(diff(DESol({(r^4+r^2)/r^2/(r^2+1)*diff(_Y(r),$(r,2))+(r^7+2*r^5+r^3+2*r)/r^2/(r^2+1)*diff(_Y(r),r)+(2*r^6+4*r^4-10*r^2-6)/r^2/(r^2+1)*_Y(r)},{_Y(r)}),r)*r+3*DESol({(r^4+r^2)/r^2/(r^2+1)*diff(_Y(r),$(r,2))+(r^7+2*r^5+r^3+2*r)/r^2/(r^2+1)*diff(_Y(r),r)+(2*r^6+4*r^4-10*r^2-6)/r^2/(r^2+1)*_Y(r)},{_Y(r)}))/r^2/(r^2+1) DESol是无解析解 >> [x,y]=dsolve('Dx+2*x/r=Dy/2','3*Dy+2*y*r=2*Dx','r') %f=r的解 x = exp(-1/2*r^2)*(r^2+3)*(C1+C2*Int(1/(r^2+3)^2*r^4*exp(1/2*r^2),r))/r^3 y = 2*(4*r^2+3+r^4)/r^3/(r^2+3)*C2*exp(-1/2*r^2)*Int(1/(r^2+3)^2*r^4*exp(1/2*r^2),r)+2*(4*r^2+3+r^4)/r^3/(r^2+3)*C1*exp(-1/2*r^2)-2/(r^2+3)*C2 朱顶余:你还差远了,别吹牛了,全部都是错误答案。 沈建其:我怀疑你不会使用数学软件。 ※※※※※※ 天地之道,以阴阳二气造化万物。是故易有太极,是生两仪。两仪生四象,四象生八卦。 |
|
对【96楼】说: 梁斌斌, X'+2X/r=Y'/2, 这个方程组真的没有解析解么?当f=k/r时,该方程组具有系列解析解 你用数学软件处理的结果如何? |
|
>> [x,y]=dsolve('Dx+2*x/r=Dy/2','3*Dy+2*y*k/r=2*Dx','r')
x = C1*r^(-1/2*k-3/2+1/2*(k^2-2*k+9)^(1/2))+C2*r^(-1/2*k-3/2-1/2*(k^2-2*k+9)^(1/2)) y = (C1*r^(-1/2*k-3/2+1/2*(k^2-2*k+9)^(1/2))*k-3*C1*r^(-1/2*k-3/2+1/2*(k^2-2*k+9)^(1/2))-C1*r^(-1/2*k-3/2+1/2*(k^2-2*k+9)^(1/2))*(k^2-2*k+9)^(1/2)+C2*r^(-1/2*k-3/2-1/2*(k^2-2*k+9)^(1/2))*k-3*C2*r^(-1/2*k-3/2-1/2*(k^2-2*k+9)^(1/2))+C2*r^(-1/2*k-3/2-1/2*(k^2-2*k+9)^(1/2))*(k^2-2*k+9)^(1/2))/k 给定初值条件能解出特解,这是通解。 ※※※※※※ 天地之道,以阴阳二气造化万物。是故易有太极,是生两仪。两仪生四象,四象生八卦。 |
|
对【101楼】说: 如果 我的破题技巧也能被数学软件所取替,那还算是什么绝技呢? x^(1/2)y"=y^(3/2) 这道著名的 汤马斯-费米(Thomas-Fermi)方程;你用数学软件所获取的解析解的形式 是什么样子的?无能为力了吧? |
|
梁彬彬、朱顶余,
你们是赌了谁输就封号1个月的,别人不好为你俩做裁判,自裁吧! 如果是朱顶余输了,就封朱顶余一个月完事,我也不说什么了! 如果是梁彬彬输了,就封梁彬彬一个月、且梁辞去版主。然后就裁定103楼这两个方程,无庸说,是我赢了,版主是我的,我爱封谁就可以封谁了。 曾云海 |
|
曾云海:
我问他意下如何,他没说啊,就是他不敢赌。跟我赌,朱早就输了N回了。 原来是冷落你了,你想插句话是不是?那你可以随便插,我可不奉陪。 ※※※※※※ 天地之道,以阴阳二气造化万物。是故易有太极,是生两仪。两仪生四象,四象生八卦。 |
|
对【109楼】说: 沈博士,你先阅读[94楼]: [楼主] [94楼] 作者:541218 发表时间: 2011/02/19 19:22 [加为好友][发送消息][个人空间]回复 修改 来源 删除 …………………………………………………………………………………………………………………………………………… 这一环很重要!必须首先将其试探函数的类型明确下来!否则,下面的环节就无法进行…… |
|
对【96楼】说: xhzjzs 用f=r^2和f=k/r为例用计算程序算了解,这个我也会的(Maple, Matlab我都会)。 但我要的是f待定的情形。这个“f待定”对我很重要。计算机恐怕无法给出通解。我是寄希望于手工。我相信手工解是存在的。朱先生改动了一个系数(将第一个X'改为2X'),得到的解析解好像是对的。如果不改动系数呢?难道真的就影响很大,以致无解析解? |