Mathematica求数值解,错在哪?h = 4.13566743*(10^-15)m = 0.1*9.109*(10^-31)Do[NSolve[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] + Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Cos[Sqrt[8*Pi^2

来源:学生作业帮助网 编辑:作业帮 时间:2024/11/29 15:35:25
Mathematica求数值解,错在哪?h = 4.13566743*(10^-15)m = 0.1*9.109*(10^-31)Do[NSolve[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] + Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Cos[Sqrt[8*Pi^2
xTOVW,Vv2&4_H&J6!Ж6c@mmAH$$񵝧 ;Fځ{9xnKҵd[tߢR;n<\Xd8^3(Ïd@L #:J9C٤jڜ^KdX8B ^aqMFVO(>ḛx0 <?}*`fa7yj7a#0Z^Ѐ8NΧL8ORB MB*Nb0}?GJ?ܭ#]?qvewA&,6i:5jzs:?v@~Vvί8ޛv.Fj'h5ڛSRrvypVyuv!S.C 

Mathematica求数值解,错在哪?h = 4.13566743*(10^-15)m = 0.1*9.109*(10^-31)Do[NSolve[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] + Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Cos[Sqrt[8*Pi^2
Mathematica求数值解,错在哪?
h = 4.13566743*(10^-15)
m = 0.1*9.109*(10^-31)
Do[NSolve[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*
Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*
Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] +
Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*
Cos[Sqrt[8*Pi^2*m*y/h^2]*x/2] - Cos[x] == 0,{x,y}],{x,0,10,
1}]

Mathematica求数值解,错在哪?h = 4.13566743*(10^-15)m = 0.1*9.109*(10^-31)Do[NSolve[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] + Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Cos[Sqrt[8*Pi^2
1 你的m的赋值式里面用了中文标点.
2 Do是不会主动显示循环中的计算结果的,要使用Print语句,当然最方便的改法还是把Do改成Table.
3 最关键的一点,这个超越方程恐怕不是Solve和NSolve能解的了的.可以考虑用FindRoot,不过这就需要你对解的范围有大概的了解,要了解解的范围,一个可能的方法,是使用ContourPlot观察解的位置:
h = 4.13566743*(10^-15);
m = 0.1*9.109*(10^-31);
ContourPlot[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] + Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Cos[Sqrt[8*Pi^2*m*y/h^2]*x/2] - Cos[x] == 0, {x, 0, 11}, {y, 0, 20}, PlotPoints -> 30]



然后就根据图上显示的解的范围去寻根吧.

mathematica里面求方程的数值解怎样才可以得到具体的数值解啊 Mathematica求数值解,错在哪?h = 4.13566743*(10^-15)m = 0.1*9.109*(10^-31)Do[NSolve[(0.4 - 2*y)/(2*Sqrt[0.4*y - y^2])*Sinh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Sin[Sqrt[8*Pi^2*m*y/h^2]*x/2] + Cosh[Sqrt[8*Pi^2*m*(0.4 - y)/h^2]*x/2]*Cos[Sqrt[8*Pi^2 mathematica 求微分方程数值解 并画图求微分方程 x'(t)=x(1-0.1y);y'(t)=-0.5y+0.02xy;x(0)=25;y(0)=2在【0,15】的数值解,并画出的图形请把程序代码告诉我.急.求. mathematica 画微分方程数值解求y'(t)=y^2+t^2,y(0)=0的解的图像! mathematica 微分方程数值解的引用我解出了微分方程的数值解,想画出这个数值解得积分,怎么办?不行,,,ganwei1989的答案还是不行.... mathematica与matlab计算比较?各自优点在哪里?仅纯数值计算而言,谁更适合? 求mathematica大神帮助 mathematica怎么求微分方程 MATHEMATICA 8.0 求翻译 问个Mathematica求方程数值解的问题FindRoot[psi[r] == 12,{r,1},WorkingPrecision -> 10]Out[168] {r -> 17.58534817}解是解出来了,但怎么把解自动赋值给其他变量,以供后续运算 使用mathematica解常微分方程(数值解)方程是t为自变量,c、k、h为常数.取正号. mathematica 我需要在二维坐标系中描出100来个数据点,并最佳拟合.同时对一部分线性拟合.我在mathematica中输入之后,总是报错,二维图都画不出来,求大神指导下,怎么写编码,错在哪了? Mathematica求约束函数最大值,哪里错了? 用matlab或者mathematica求积分解微分方程 mathematica怎么求导比如求sinβ导数,其中β是时间t的函数,正确解为cosβ*β(一点),一点在β上边,这里打不出来,在mathematica里怎么表示 我要用mathematica求一个方程组的数值解,可是应该是程序有问题出不来结果,麻烦你帮我看看,是要求这个方程组的x'[t]=-ax[t],y'[t]=x[t],x'[0]=0,y'[0]=0,求它的数值解,可以的话再画出图来,其中a是常数, 如何用Mathematica实现复合函数的反解?比如,知道g[x]的表达式,知道f[g[x]]的最终表达式子,如何求f[x]?在Mathematica上实现 求数值微分的mathematica方法有这样一组数列{1,1,1,2,3,5,7,10,13,15,16,17,18,18,18},请对它求导.只能用mathematica,任意版本.