想求解第一类一介贝塞尔函数的前10个根，根用于xi表示，i从1到10，请问要怎么做呢，我是初学者，这个代码解决不了Solve[BesselJ[1,x]==0,x]，希望能帮

比如2x=3y，或者2x<3y，可以通过某个命令调整成2x—3y=0和2x—3y<0嘛

近期发现的bug一个（已返馈给WRI），行为太诡异了，姑且发一帖说说： f1 = Function[f, If[f > 1/2, Red, Blue]]; f2 = (If[# > 1/2, Red, Blue] &); SliceDensityPlot3D[Sin[π x], "BackPlanes", {z, -2, 2}, {y, -2, 2}, {x, -2, 2}, ColorFunction -> #] & /@ {f1, f2} 是的你没看错，SliceDensityPlot3D仅识别出了以 # 号表示的纯函数，而对等效的“有名字”的纯函数无动于衷……

因为Cambria Math字体对字符采用了特殊编码。关于这个可以参考 mathematica.stackexchange.com/q/288212/1871 mathematica.stackexchange.com/q/193668/1871

RootSearch::err: The Private RootSearch function zBent had an error in branch 4 when searching for a root near Ersek`RootSearch`Private`newX. Please email a notebook containing this use of the RootSearch function to RootSearch author at (ted.ersek@navy.mil). 这是我的代码 Needs["Ersek`RootSearch`"] seg = g -> 9.8; seh = h -> 4000; Subscript[se\[Rho], 1] = Subscript[\[Rho], 1] -> 1024;(*水的密度*) Subscript[se\[Rho], e] = Subscript[\[Rho], e] -> 917; sec = c -> 1500; seE = capse -> 5*10^9; sev = v -> 0.3; seT = T -> 0; se\[CurlyEpsilon] = \[CurlyEpsilo

背景：已知钢材本构，求矩形截面弯矩与曲率关系（平截面假定） 代码描述：定义钢材本构sigma[e]，为四折线函数，用k*y代替e(k为曲率，y为截面坐标)，求sigma*y对y取[-h/2,h/2]的积分 问题： 1.常数符号较多，进行全局假设后，积分仍无法进行，提示无法确定[-h/2,h/2]为实数 2.放弃全局假设，在积分中引入假设，结果中出现很多不符合假设的分段条件 代码： 问题1： SetAttributes[{ey, eh, em, eu, Es, Eh, Eu, h}, Constant}] $Assumptions = {0 < ey < eh < em < eu, Es > 0, E

主要是想得到步长为0.1的csv数据，然后得到结果是这样的，请问是程序出错还是结果就是这样 y = {\[Kappa] = 1; \[Epsilon] = 0.02*\[Kappa]; \[Eta] = 0.88; \[Mu]2 = 0.99*\[Kappa]; \[Mu]1 = 0.96*\[Kappa]; \[Lambda] = 7/10^5; Log[(((2*(\[CapitalDelta] + \[Eta]) - I*\[Kappa])^2 - 4*\[Mu]1*\[Mu]2)*((2*(\[CapitalDelta] + \[Eta]) + I*\[Kappa])^2 - 4*\[Mu]1*\[Mu]2)* (16*\[CapitalDelta]^4*\[Lambda] - 16*I*\[CapitalDelta]^3*(\[Epsilon]^2 + 5*I*\[Eta]*\[Lambda] + 2*\[Kappa]*\[Lambda]) - 2*I*\[Epsilon]^2*((2*\[Eta] - I*\[Kappa])^2*(4*\[Eta] - I*\[Kappa]) + 4*\[Eta]*\[Mu]1*\[

R = 0.045;L = 14;r = 0.009;e = 0.2;n = L/(\[Pi]*R*R);v = 4/3*Pi*r^3;g = (1 - 7/16*v)/(1 - v)^2; G = g*v;H = 1 + 1/(2*G);M = 1 + (9*\[Pi])/16*(1 + 2/(3*G))^2;k = 2/Sqrt[\[Pi]]*n*2*r*G*M*Sqrt[T[x, y]];\[CapitalGamma] = 4/Sqrt[\[Pi]]*(1 - e)*n*G*T[x, y]^1.5*1/(2*r);CirR = ImplicitRegion[x^2 + y^2 <= R^2, {x, y}];(*求解出n和T的分布，自变量为x和y，其中边界条件为：T在边界上沿径向偏导为0，n在区域内求积分等于L*)eqn1 = Simplify[{ Grad[T[x, y], {x, y}] == {GTx[x, y], GTy[x, y]}, (*Inactive[Div][k[x,y]*Inactive[Grad][T[x,y],{x,y}],{x,y}]*) k[x, y]*Div[{GTx[x,

8名学生参加一次数学测验，测验共有8道题，已知每道题都有至少4名学生答对，问是否一定有两个学生共同答对了所有的问题？

kcube = Graphics3D[{EdgeForm[{Thickness[0.005], Blue}], FaceForm[], Cube[{0, 0, 0}, N[2*persat*Sqrt[1/3]]]}]; 请问Cube[{0, 0, 0}, N[2*persat*Sqrt[1/3]]]}]是什么意思？

我现在只要打开mma，即便啥也不干，电脑就会偶尔卡一下。我发现只要我重新打开mma，第一次右键桌面都会卡一下。但平时不开mma就不会有这种问题。而且这个卡顿是最近出现的。

我定义了一个分段函数Lens[x_,y_],其中包含了一个Subscript[Y, i ],i取不同值时Y[]是不同的函数。但是某些Y[]的计算数值会产生一个很小的虚部 0.000000i，在运行我的分段函数时会提示我 “Less::nord: 尝试使用 37344.8 -3240.47 I 进行无效的比较.” 之后我将我的分段函数Lens[x_,y_]化为了数组Lensarray然后用MatrixPlot函数画图。画的图是正常的。 虚部按理来说是个及其小的值，如果我将Y[]前面加上Abs[]或者用Chop[]函数最后画出来的图应该是没有影响的，但是实际上画出来

\[Kappa] = 1; \[Epsilon] = 2/100; \[Lambda] = 7/10^5; Solve[{-16*Sqrt[2]*\[CapitalDelta]^3*\[Epsilon]^2 - 64*Sqrt[2]*\[CapitalDelta]^2*\[Epsilon]^2*\[Eta] - 80*Sqrt[2]*\[CapitalDelta]*\[Epsilon]^2*\[Eta]^2 - 32*Sqrt[2]*\[Epsilon]^2*\[Eta]^3 + 12*Sqrt[2]*\[CapitalDelta]*\[Epsilon]^2*\[Kappa]^2 + 16*Sqrt[2]*\[Epsilon]^2*\[Eta]*\[Kappa]^2 - 32*Sqrt[2]*\[CapitalDelta]^3*\[Kappa]*\[Lambda] - 120*Sqrt[2]*\[CapitalDelta]^2*\[Eta]*\[Kappa]*\[Lambda] - 144*Sqrt[2]*\[CapitalDelta]*\[Eta]^2*\[Kappa]*\[Lambda] - 56*Sqrt[2]*\[Eta]^3*\[Kappa]*\[Lambda] + 8*Sqrt[2]*\[CapitalDelta]*\[Kappa]^3*\[Lambda] + 10*S

k = Cos[aa3]^2 + 2 Sin[aa3]^2 Simplify[k] 结果还是Cos[aa3]^2 + 2 Sin[aa3]^2 而不是1+Sin[aa3]^2 即使用FullSimplify　结果还是一样，哪位大神能帮帮我啊。。。

当CellContext为默认的Global时，变量名可以自动补全，但是当我使用 SetOptions[EvaluationNotebook[], CellContext -> Notebook] 将上下文设置为当前笔记本时，自定义变量的输入自动补全功能就失效了。 我又试了通过Context[]查看当前上下文，然后手动输入该上下文然后再接着输入变量名，这时自动补全功能又起作用了。 请问这是新版本的bug吗？

SetAttributes[{a, b, q}, Constant]; eqn = {Laplacian[w[x, y], {x, y}] == v[x, y], Laplacian[v[x, y], {x, y}] == 7}; bcs = {w[0, y] == 0, w[x, 20] == 0, w[10, y] == 0, w[x, 0] == 0, Derivative[1, 0][w][0, y] == 0, v[x, 20] == 0, Derivative[1, 0][w][10, y] == 0, v[x, 0] == 0}; ufun = NDSolveValue[{eqn, bcs}, w, {x, 0, 10}, {y, 0, 20}] NDSolveValue::fembdnl: The dependent variable in (w^(1,0))[0,y]==0 in the boundary condition DirichletCondition[(w^(1,0))[0,y]==0,x==0.] needs to be linear. 怎么搞成线性的啊

如题，想问问bug出在哪了，代码如下： Clear["Global`*"]; Remove["Global`*"]; Get["VariationalMethods`"]; x0[t_] := Cos[omega*t]; y0 = 0; xx1[t_] := x0[t] + l*Cos[theta1[t]]; yy1[t_] := (-l)*Sin[theta1[t]]; ek[t_] := 0.5*m*(Derivative[1][xx1][t]^2 + Derivative[1][yy1][t]^2); ep[t_] := m*9.81*yy1[t]; lag[t_] := ek[t] - ep[t]; equ = EulerEquations[lag[t], theta1[t], t]; m = 1; l = 2; omega = 2; ini = {theta1[0] == 0., Derivative[1][theta1][0] == 0}; sol = NDSolve[{equ, ini}, theta1, {t, 0, 50}, Method -> {"IndexReduction" -> "Pantel

对于一个方程我用ContourPlot作出了2D图 nmax = 40; c[n_] := v (1 - (-1)^n) *2/(n*Pi*Sinh[n*Pi*b/a]); \[Phi] = Sum[c[n]*2/a*Sin[n*P

用的是manipulate函数画的plot3D。不仅是这个图，还有一个图也是这样显示不全

尝试分别用ArcLength函数或者用原本的弧长积分公式计算弧长，结果是能出来的，但是耗时起码3分钟。 结果与三维软件UG的曲线对照，计算出的弧长是一致的。电脑是内存8G的。 版本是mathematica10.3 代码如下，以及截图 绘制曲线也是立刻出来的，就是弧长算半天。。。 曾尝试用标准圆来测试ArcLength函数，也是立刻得出弧长的。。。 ClearAll; a = 87; amp = 1; i21 = -(1/2 + amp*(Sin[t/2])^2); rv[t_] := a*i21/(i21 - 1); rw[t_] := a - rv[t]; xv[t_] := rv[t]*Cos[t]; yv[t_] := rv[t]*Sin[t]; xw[t_]

是这样的，现在算出来了一个巨大的矩阵。它内存大小大概是10000乘30000 我看了一下占了8g内存 电脑运行内存16g然后后面下一步我还要对这个矩阵求转置。现在是一求转置电脑就卡死了 内存占用100%。 想问问吧里的大佬这个问题怎么解决啊换个电脑吗

真是被气的血压飙升 SetAttributes[B, Constant] SetAttributes[n, Constant] SetAttributes[C, Constant] SetAttributes[A, Constant] FindMaximum[{-(1/( 4 (-1 + n))) (8 + 6 B^2 + 8 C + 8 C^2 - 5 n - 4 B^2 n - 4 C n - 5 C^2 n + (8 + 4 B^2 + 8 C - 4 n - 2 B^2 n - 4 C n) Cos[2 x] + n Cos[4 x] + (4 B C - 2 B C n) Cos[ x - 3 y] + (2 B^2 + 8 C - B^2 n - 6 C n) Cos[ 2 x - 2 y] + (16 B + 16 B C - 10 B n - 10 B C n) Cos[ x - y] + (4 B - 2 B n) Cos[ 3 x - y] + (4 B^2 + 8 C + 8 C^2 - 2 B^2 n - 4 C n - 4 C^2 n) Cos[2 y] + C^2 n Cos[4 y] + (12 B + 12 B C - 6 B n - 6 B C n) Cos[x + y] + 2 B n Cos[3 x + y]

(*第二题*) Clear[l, m, a, b, \[Theta]1, \[Theta]2, L, funcs, sol]; Clear["Global`*"]; {g, m, A, \[Omega]0} = {9.81, 1, 2, 1}; l = g/(\[Omega]0^2)； Manipulate[ x = A*Cos[\[Omega]*t] + l*Sin[\[Theta][t]]； y = -l*Cos[\[Theta][t]]; funcs = {D[D[x, t], t]*Cos[\[Theta][t]] + D[D[y, t], t]*Sin[\[Theta][t]] == g*Sin[\[Theta][t]], \[Theta][0] == 0, \[Theta]'[0] == 0}; sol = Quiet@ NDSolve[funcs, \[Theta], {t, 0, 10}, Method -> {"IndexReduction" -> {True, "ConstraintMethod" -> "Projection"}}];, {\[Omega], 0.5, 3}] p[t_] = {A*Cos[\[Omega]*t] + l

Clear[l, m, a, b, \[Theta]1, \[Theta]2, L, funcs, sol]; Clear["Global`*"]; {g, m, A, \[Omega]0} = {9.81, 1, 2, 1}; l = g/(\[Omega]0^2)； Manipulate[\[Omega], {\[Omega], 0.5, 3}] x = A*Cos[\[Omega]*t] + l*Sin[\[Theta][t]]； y = -l*Cos[\[Theta][t]]; funcs = {D[D[x, t], t]*Cos[\[Theta][t]] + D[D[y, t], t]*Sin[\[Theta][t]] == g*Sin[\[Theta][t]], \[Theta][0] == 0, \[Theta]'[0] == 0}; sol = Quiet@ NDSolve[funcs, \[Theta], {t, 0, 10}, Method -> {"IndexReduction" -> {True, "ConstraintMethod" -> "Projection"}}]; p[t_] = {A*Cos[\[Omega]*t] + l*Sin[\

代码的文本是Solve[Cos[x]＝＝x,x]

请问为什么Mathematica中对同一个函数级数展开的一次项的系数和求偏导后取相同x0得到的值会不同In[907]:= ClearAll[ff, r, M, A, B, c, ddc, ddA, \[Alpha], rm, g, z, \r0, h, p, q, a, b, br, bd, d, um, R, hh, RR1, RR, rr1, rr, fff]M = 1/2;ff[r_] = 1 - 1/r + \[Alpha]/r^2;A[r_] = ff[r];B[r_] = 1/ff[r];c[r_] = r^2;ddc[r_] = D[c[r], {r, 2}];ddA[r_] = D[A[r], {r, 2}];dA[r_] = D[A[r], {r, 1}];rm = Solve[D[c[r], {r, 1}]/c[r] == D[A[r], {r, 1}]/A[r], r][[2, 1, 2]]; // Simplifyrm = 3/4 (1 + Sqrt[1 - 32 /9 \[Alpha]]);r = (r0^2 + Sqrt[ r0^4 - 4 r0^2 \[Alpha] (r0 - r0 z - \[Alp

我们知道在复平面中解析区域内取任何一个回路进行积分，当回路中没有奇点时积分结果应该为0。我们取被积函数func[\[Beta]_, i_] := \[Beta]^i/(2 \[Pi] I )，其中含一个指数项。我们选取的路径保证这个函数始终没有任何奇点。当选取单位圆积分时，结果如图1，没有问题；当选取一个平移后的半径为1的圆时（图2），积分结果在i较大时出现了一个拐点，随后保持相同的斜率增大，造成了明显的数值问题（图3）。根据个人经验，原因应该是积分路径一部分在

Eigensystem求解本征值默认按照绝对值大小排序，请问如何让本征值按照实部从小到大排序，对应的本征矢量也和本征值顺序对应。 （因为要求解一个随连续参数变化的矩阵的本征值问题，这样做是想得到一个类似随参数解析的表达式，每个不同值都解出来排序一次不能满足连续性）

DSolve[{e*i*(v'''')[x] == p0, v[0] == 0, v[l] == 0, (v'')[0] == 0, (v'')[l] == 0}, v[x], x]

这个表达式在Sin[2 n2 pi/M]≠0的时候应该可以化简成1，但是不知道为什么我化简不成功，所以想请教一下大家有没有什么解决办法。 另外还想问问assume和assumption除了格式不一样以外，有没有其他的区别？ 谢谢！ Assuming[Sin[(2 n2 \[Pi])/M] != 0 && n2 \[Element] PositiveIntegers && M \[Element] PositiveIntegers, Csc[(2 n2 \[Pi])/M] Sin[2 (2 - 1/M) n2 \[Pi]]] Simplify [Csc[(2 n2 \[Pi])/M] Sin[2 (2 - 1/M) n2 \[Pi]], Assumptions -> {Sin[(2 n2 \[Pi])/M] != 0, n2 \[Element] PositiveIntegers, M \[Element] PositiveIn

如题。我想知道在Mathematica中，如何统计一段代码在运行时，总共进行了多少次浮点运算。 不过我也不知道我的描述是否有问题，或者Mathematica是不是真的做到这件事，所以这里简单说一下来龙去脉。 事情是这样的，我的一篇论文中，算法部分的代码是用Mathematica实现的。现在主编对于我关于计算机这一部分的讨论不满意，希望我能用roofline模型，对代码进行定量的评估，以此来判断代码性能是否已达到我计算机的极限。 我之前并没有接触过什么是roof

b = 10.0; n = 2.0; m = 20; a = {}; n0 = n; n1 = b; Do[ak = 0; n2 = n0; While[n2 >= n1, n2 = n2/n1; ak++], AppendTo[a, ak]; n0 = n1; n1 = n2; If[N[Abs[n2 - 1]] < N[10^(-m)], Break[]], {k, 1, m}]; len = Length[a]; y = a[[len]]; Do[y = a[[i]] + 1/y, {i, len - 1, 1, -1}]; Print[a]; Print[{m, y, N[y, m]}]; N[Log[10, 2], m] Do::nliter: 位置 2 处的非列表迭代器(non-list iterator) AppendTo[a,ak];n0=n1;n1=n2;If[N[Abs[n2-1]]<N[10^-m],Break[]] 没有算出实数值.

最近在做一个方程参数拟合的东西，是通过相平衡来实现的，简单来说就是输入温度和压力，输出两个组

不求代码，只求大神能够给个简单的思路，非常感谢。题目原文在此，Mathematica中，FindRoot命令能够求出方程的根，但每次使用FindRoot，只能求出一个根，试利用FindRoot函数，编写新的函数FRoots，它能够求出单变量方程在某个范围的所有根，并要求函数FRoots返回的根不重复，最后，利用函数FRoots求出下列方程的根：

用findroot命令解方程组，求得的结果画出图像来明显不满足其中的一些方程（具体表现在有几个方程是利用连续性条件得到的，但是最终的结果画出图像来很明显不连续），我自己觉得可能是求解精度不够，所以希望有人能帮忙看看指点指点，刚用mathematica不久，谢谢大家。代码如下（复制过来有些符号显示不出来了）： apation={a1*((c1+c2)/(a1+a3))+c4==0,-a3*((c1+c2)/(a1+a3))+c6==0,c3-c5==0,c4-c6==a2/(a1+a3)*\[Tau]m,\[Beta]1*r*h==ks*((a5+a6)/6*r*h^3-a2/2*\[Tau]m*h^2)+(c9-c11)*h+c10-c12,a1*(c

\[Alpha]=0.5;\[Beta]=1;\[Rho]=1;\[Lambda]=0.01;b=-1;a=-0.18; M=-((a(a^2+b^2))/(b^2 \[Lambda])); f=M+(x-a t)^2+b^2 t^2; g=1/((a \[Alpha]^2+\[Lambda])^2+b^2 \[Alpha]^4) 4(I(a^2 \[Alpha]^3+b^2 \[Alpha]^3+a \[Alpha] \[Lambda])x-I \[Alpha](a^2+b^2)(a \[Alpha]^2+\[Lambda])t-\[Alpha]^2 (a^2+b^2))+f; A=\[Rho] Exp[I \[Alpha] x - (I \[Lambda] t)/\[Alpha]]*g/f; Plot3D[A*(Conjugate[A]),{x,-10,10},{t,-10,10},PlotPoints->60,AxesLabel->{"x "," t"," \[VerticalSeparator]A\[VerticalSeparator]^2"},PlotRange->All,AxesEdge->{{-1,-1},{-1,-1},{-1,-1}},Ticks->{{-5,0,5},{-3

DSolve[ { -y''[t] == (G*m^2)/(x[t]^2 + y[t]^2)*y[t]/Sqrt[x[t]^2 + y[t]^2], -x''[t] == (G*m^2)/(x[t]^2 + y[t]^2)*x[t]/Sqrt[x[t]^2 + y[t]^2] }, {x[t], y[t]}, t]

代码如下： feq1 = {D[u[t, x], t, t] + D[u[t, x], x, x, x, x] == 0}; feq2 = {D[v[t, x], t, t] + D[v[t, x], x, x, x, x] == 0}; bc1 = {v[t, 0] == 0, Derivative[0, 1][v][t, 0] == 0, Derivative[0, 2][v][t, 1] == 0, Derivative[0, 3][v][t, 1] == u[t, 1]}; bc2 = {u[t, 0] == 0, Derivative[0, 1][u][t, 0] == 0, Derivative[0, 2][u][t, 1] == 0, Derivative[0, 3][u][t, 1] == If[t < 1, t, 1]}; st1 = {u[0, x] == 0, Derivative[1, 0][v][0, x] == 0}; st2 = {v[0, x] == 0, Derivative[1, 0][u][0, x] == 0}; sol = NDSolve[{feq1, feq2, st1, st2, bc1, bc2}, {v, u}, {t, 0, 10}, {x, 0, 1}]; Plot[{Evaluate[v[t, 1]

代码如下 e={1.131,1.33,1.594,1.33,1.131,1.33,1.594,1.33,1.131,1.33,1.594,1.33,1.131,1.33,1.594,1.33,0.605,0.9468,2.668,1.222};(*各膜层厚度倍率*) n={2.3,1.63,1.38,1.63,2.3,1.63,1.38,1.63,2.3,1.63,1.38,1.63,2.3,1.63,1.38,1.63,2.3,1.63,1.38,1.63};(*各膜层折射率*) \[Theta]=Table[ArcCos[(1-1.52*1.52*(Sin[Pi/4]/n[[j]])^2)^0.5]/Degree,{j,Length[n]}];(*各膜层入射角*) d=Table[e[[j]]*550*n[[j]]/4,{j,Length[n]}];(*各膜层厚度*) \[Eta]p=Table[n[[j]]/Cos[\[Theta][[j]] Degree],{j,Length[n]}];(*p偏振光的\[Eta]*) \[Eta]s=Table[n[[j]]*Cos[\[Theta][[j]] Degree],{j,Length[n]}];(*s

本人萌新，想用mma做一些光学的仿真

Eliminate[{x^2 + y^2 + z^2 - 2 Rz == 0, x^2 + y^2 + z^2 - R^2 == 0 },z] // Simplify

NetTrain这个函数在我的i7-12700H上面不能跑满CPU(我已经在Windows11和linux上面都进行过测试了) 但是在R7-5800H上面可以跑满整个CPU(测试了windows11) 请问有人知道是什么问题吗 代码如下 dat = Flatten@ Table[{x, y} -> Exp[-Norm[{x, y}]], {x, -3, 3, .005}, {y, -3,3, .005}]; ne = NetChain[{32, Tanh, 1}]; NetTrain[ne, dat, All]

数据从excel导入的

想求 x^5 - 5 x - 2 = 0 的公式解，但是这方程不存在由根式构成的公式解。 而 mathematica 给出了以下解答，5 个

\[Psi] = ({ {a E^(I \[CurlyPhi]1)}, {b E^(I \[CurlyPhi]2)}, {c E^(I \[CurlyPhi]3)}, {d E^(I \[CurlyPhi]4)} }) ConjugateTranspose[\[Psi]]