1 为什么在数学助手里输入的(d/dx)无法解析 2 高阶导除了用D[f（x），x，n]还有什么办法 3 怎么让它自动化简计算结果 4 求解时复数结果为什么是诸如（-1）^(1/4)这种形式，怎么让它表示成a+ib形

这问题其实很有趣啊。Dialog这个函数，实际上就是中断计算（Alt+,）后点击“使用子会话”时被调用的函数，它会形成一个对话（Dialog），在这个对话里还真可以干些有意思的事。 假设我们有个很费时间的Table： Table[Pause[1]; i, {i, 10}](*我就用Pause意思一下了*) 我们可以Alt+，使用子对话，然后，查看Table算到哪里了！： 全体使用了Block进行变量局部化的函数，都可以用这个方法来查局部变量的值。 Module也是可以查的，不过麻烦一点： 在对话中对变量进行

各位大佬，我想要利用mathematica求方程的特解。可是程序跑完之后没有output，连空集都没有。然后左侧的行号直接变成了In[·]，有大佬知道为啥嘛 ClearAll["Global`*"] gin@2@2 = s11@l2@2 + ( s21@l2@2^2*((54.548 - 16.455*I - 50)/(54.548 - 16.455*I + 50)))/( 1 - s22@l2@2*((54.548 - 16.455*I - 50)/(54.548 - 16.455*I + 50))); zin@2@2 = 50*(1 + gin@2@2)/(1 - gin@2@2); yin@2@2 = 1/zin@2@2; FindInstance[ Re[yin@2@2] == 0.02, {s11@l2@2, s21@l2@2, s22@l2@2}, Complexes, 4]

有没有大佬知道，怎么通过mathematics把传递函数转换为 首一标准型 和 尾一标准型 例如下面的 (4 s - 4)/(2 s + 3 s^2 + s^3) 传递函数 我图里面已经给了首一型的方法，但不知道尾一型的方法，求

se\[Gamma] = \[Gamma] -> 0.9; seg = g -> 9.8; Subscript[se\[Rho], 1] = Subscript[\[Rho], 1] -> 1024; sev = v -> 0.3; Subscript[se\[Rho], e] = Subscript[\[Rho], e] -> 0.9; seE = E -> 6000; seH = H -> 100; fValues = Range[0.01, 100.01, (100.01 - 0.01)/100]; s22 = ((g*\[Omega]^2)/ H - (g k/ H ((1 + (d^3 E (k)^4)/(12 (1 - v^2))) Tanh[ 0.2` k] + (1 + (d^3 E (k)^4 \[Gamma])/( 12 (1 - v^2) )) Tanh[0.8` k] + k (1 - \[Gamma]) Subscript[\[Rho], e] Tanh[0.2` k] Tanh[ 0.8` k]) (1 + \[Sqrt](1 - (4 (1 - \[Gamma]) (1 + ( E (k)^4)/(12 (1 - v^2))) Tanh[0.2` k] Tanh[ 0.8` k] (1 + \[Gamma] T

简单来说就是 {{1, 2}, {1, 3}, {1, 4}, {1, 5}, {2, 3}, {2, 4}, {2, 5}}， 我是想将这些第二个位置的数一样的坐标放一起， 我用的是 Sort[{{1, 2}, {1, 3}, {1, 4}, {1, 5}, {2, 3}, {2, 4}, {2, 5}}, #1[[2]] == #2[[2]] &] 结果 {{2, 3}, {1, 5}, {2, 5}, {1, 4}, {2, 4}, {1, 3}, {1, 2}} 这是为什么呢？按道理不应该是 {{1, 5}, {2, 5}, {1, 4}, {2, 4}, {1, 3}, {2, 3}, {1, 2}}吗

之前看到过一次，找不到了，感觉很有意思

Clear["Global`*"] mu=SetPrecision[0.005,20]; ms=SetPrecision[0.1283,20]; G=SetPrecision[3.672,20]; K=SetPrecision[59.628,20]; Λ=SetPrecision[0.6816,20]; r=SetPrecision[0.1,20]; T=SetPrecision[0.01,20]; uu=SetPrecision[0,20]; us=SetPrecision[0,20]; w=SetPrecision[0.6,20]; Mu[SS_,UU_]:=mu+(2*K*SS-4*G)*UU; Ms[SS_,UU_]:=ms-4*G*SS+2*K*UU^2; f1[Pt_,Pz_,SS_,UU_,n_]:=((BesselJ[n+1,Pt*r])^2+(BesselJ[n,Pt*r])^2)*{4*2*Exp[uu/T]*(Exp[(2*n*w+w)/T]-Exp[2*Sqrt[Mu^2+Pt^2+Pz^2]/T])*(2*G-K*SS)*Mu/(Sqrt[Mu^2+Pt^2+Pz^2]*(Exp[Sqrt[Mu^2+Pt^2+Pz^2]/T]+Exp[(2*n*w+2*uu+w)/(2*T)])*(Exp[(Sqrt[Mu^2+Pt^2+Pz^2]+

Solve[{8 D1 h + M - 2 D1 M - 2 D0 (8 h + M) + D2 (1 + 8 (-1 + 4 D0) h + M) == 0, - 2 D2 + 8 D2 h + 2 D0 (1 + 8 (-1 + 2 D1) h) - M + D1 (1 - 8 h + M) == 0, 2 (D2 + 16 D0 h - 8 D1 h + 8 (-1 + 2 D1) D2 h + M - D1 M) == 0, 4 (D1 - D2)^2 - 16 D0 (-D0 + D1 + D2 - 2 D1 D2) + P == 0}, {D0, D1, D2, h}] 显示正在运行一个多小时了，请问是正常计算吗，还是方程组存在问题？

准新手，不知道怎么描述问题才能有效搜索，所以只好提问！ 已知一个函数 f(x,y,z)的表达式，比如 f(x,y,z)=1/x + y^2 + z 另外推得一个公式 g(a,b,c,x,y,z)= (a+b)/x + (a+b) y^2 + (a+b) z + c^7+, 内部肯定包含了f(x,y,z)的表达式，其中(a+b)是个复杂表达式，请问有什么办法 将g 代换成下列形式： g(a,b,c, x, y, z) = (a+b) f(x,y,z) + c^7 说明： 我是期望能够保持 f(x,y,z), 做了这样的反向代换，以简化整个表达式。 （本人准新手，基本操作了解，对mathematica缺少整体了解） 请赐教！

现在在弄RNN神经网络（用来解微分方程），然后定义损失函数里需要输出端对输入端求2阶导数，mma自带了网络端口梯度函数netportgradient[]，也就是1阶导，现在2阶导怎么办呢？我看见论文说是用自动求导技术，这个部分mma能实现嘛不然还得花大把时间学py

如题，我想要画一个正方体沿着正中间的正六边形斜着切开的图形，但是在网上没找到，只找到了 ConvexHullMesh 用来画凸包，但是我不想要多面体面内多余的线，请问有什么方法解决？ ConvexHullMesh[{{0, 0, 0}, {1, 0, 0}, {0, 1, 0}, {0, 0, 1}, {0, 0.5, 1}, {1, 0.5, 0}, {0.5, 0, 1}, {0.5, 1, 0}, {0, 1, 0.5}, {1, 0, 0.5}}] （鄙人小白，用的是 Mathematica 12）

为什么绘图这么慢，是哪个环节出问题了？ Clear["Global`*"] s0 = PauliMatrix[1]; s1 = PauliMatrix[2]; s2 = PauliMatrix[3]; s3 = PauliMatrix[4]; Ham = kx s1 + ky s2 ; DD = Eigenvalues[Ham]; VV = Eigenvectors[Ham]; region = RegionUnion[Disk[{0, 0}, 0.1], RegionDifference[Disk[{0, 0}, 0.3], Disk[{0, 0}, 0.2]], RegionDifference[Disk[{0, 0}, 0.5], Disk[{0, 0}, 0.4]], RegionDifference[Disk[{0, 0}, 0.7], Disk[{0, 0}, 0.6]], RegionDifference[Disk[{0, 0}, 0.9], Disk[{0, 0}, 0.8]], RegionDifference[Disk[{0, 0}, 1.1], Disk[{0, 0}, 1]], RegionDifference[Disk[{0, 0}, 1.3], Disk[{0

如下一个程序包函数，用于求解数独的。在.nb和.wl中使用按普通方式或者Compile来打印是确认可行的。但是用来确认的语句是Compile[{}, Do[Print@Grid@#, ##2]] &[array, Sequence @@ circulate][] 会打印所有解。因此，想改成带选项函数，返回第一个解或返回解的个数。 因为有的数度可能有上百万个解。 还是用Compile语句去改写(程序段最后用Which限定的那部分语句，之前调试时，用(**)当作注释括起来了，要使用请去掉)。 但是，去掉(**)进行调试，没有预期效果，也不

比方说一个非常简单的线性方程组： x+y=2 x-y=3 解出来之后，我想引用x和y进行下一步计算，应该怎么编程序呢？我刚入门，真诚求解各位大佬

默认都是六位，怎么增加或者减少啊？

clear["Global'*"] L == 0.1; tm == 60; k == 0.5; m1 == 0.02; m2 == 0.02; m3 == 0.02; r1 = Sqrt[(x2[t] - x1[t]) \.b2 + (y2[t] - y1[t]) \.b2]; f1 = k*(r1 - L); r2 = Sqrt[(x2[t] - x3[t]) \.b2 + (y2[t] - y3[t]) \.b2]; f2 = k*(r2 - L); equ = { m1*x1''[t] = f1*(x2[t] - x1[t])/r1, m1*y1''[t] = f1*(y2[t] - y1[t])/r1, m2*x2''[t] = -f1*(x2[t] - x1[t])/r1 + f2 *(x3[t] - y2[t])/r2, m2*y2''[t] = -f1*(y2[t] - y1[t])/r1 + f2*(y3[t] - y2[t])/r2, m3*x3''[t] = -f2*(x3[t] - x2[t])/r2, m3*y3''[t] = -f2*(y3[t] - y2[t])/r2, x1[0] == -0.1, x1'[0] == 0, y1[0] == 0, y1'[0] == 0, x2[0] == 0, x2'[0] == 0, y2[0]

K1是一阶的第二类修正贝塞尔函数，积分大典上好像没找到合适的公式

我试了assuming函数： Assuming[0 < a < b , Integrate[Cos[4 pi x]/(-a Cos[2 pi x] + b), x]] 但答案仍然以a>b为准，积出了这个： (-2 b x + ((a^2 - 2 b^2) ArcTanh[((a + b) Tan[pi x])/Sqrt[ a^2 - b^2]])/(Sqrt[a^2 - b^2] pi) - (a Sin[2 pi x])/pi)/a^2 求助各位大佬，怎么让它认为b>a呢？

R=100.;a=5.;n1=n3=1;n2=3;yb=Sqrt[a*(2R-a)]; light={};\[Theta]0=0;k1={Cos[\[Theta]0],-Sin[\[Theta]0]}; Do[ p1={-R,y0}; equl={y==p1[[2]]+k1[[2]]/k1[[1]]*(z-p1[[1]]),y^2+(z-R+a)^2==R^2}; s=Solve[equl,{z,y}] m=z/.s s=Select[s,(z/.#)=Min[m]&]; p2={z,y}/.s[[1]]; k= (n1*k1[[2]]* (p2[[2]]*k1[[2]]+(p2[[1]]-R+a)*k1[[1]])+(n2-n1)*p2[[2]])/(n1*k1[[1]]* (p2[[2]]*k1[[2]]+(p2[[1]]-R+a)*k1[[1]])+(n2-n1)*(p2[[1]]-R+a)); kz=1/Sqrt[1+k^2];ky=k*kz;k2={kz,ky}; equ2={y==p2[[2]]+k2[[2]]/k2[[1]]*(z-p2[[1]]), y^2+(z+R-a)^2==R^2}; s=Solve[equ2,{z,y}]; m=z/.s;s=Select[s,(z/.#)=Max[m]&]; p3={z,y}/.s[[1]]; k=(n2*k

FinancialData函数可以获取股票相关的数据，但是它不能获取实时数据，只能获取昨天或者最近一天的收盘价。 有什么办法可以获取实时股价？

鼠标滑过一次就卡死一次, 服了

我的目标是将一个矩阵元为关于x的函数的3*3矩阵里的函数提取出来分别绘图，在最后一步我尝试用Do循环解决，但是代码可以运行但没有输出，不知道为什么会这样。 下面是我的代码： L = {{x, Sqrt[x], x^2}, {1 + 3 x, 2 + x, E^x}, {Log[x], Sinh[x], Cosh[x]}}; func[11] = L[[1, 1]]; func[12] = L[[1, 2]]; func[13] = L[[1, 3]]; func[21] = L[[2, 1]]; func[22] = L[[2, 2]]; func[23] = L[[2, 3]]; func[31] = L[[3, 1]]; func[32] = L[[3, 2]]; func[33] = L[[3, 3]]; Do[{Plot[Tooltip[func[ij]], {x, 1, 6}]}, {i, 1, 3}, {j, 1, 3}]

我有一个老师写的cdf文件，打开可自动运行，输入参数，点击按键输出结果，我想提取源码，应该怎么操作

请问各位大佬，这行代码"listplot[datalist[[ ;; ;; 10 ]]]"中的;; ;; 是从第一个元素取到最后一个元素吗？10是每10个数据画一个点？求教

Export:没有足够的 Java 堆栈空间，以完成此操作. 尝试增加 Java Virtual Machine 的堆栈大小 请问这个堆栈大小怎么调？谢谢。

如图 我知道CreateASectionHeading，但是后面的折叠符号是怎么来的？没在帮助文档里找到 谢谢！

以下是Mathematica 8.0.4帮助文档中讲Range的页面，其中inverse[[perm]]这里双方括号是什么意思没有看懂，请哪位大神来解释解释~ 因为这里inverse、perm都不是内置函数，而且已经被赋值为一组数，所以顺便也说说一个方括号在这种情况下是什么意思。 求逆排列: In[1]:= perm = RandomSample[Range[10]] Out[1]= {4, 2, 8, 1, 10, 5, 3, 6, 9, 7} In[2]:= inverse = perm; inverse[[perm]] = Range[Length[perm]]; inverse Out[2]= {4, 2, 7, 1, 6, 8, 10, 3, 9, 5}

方程中自变量量是一个未知函数怎么定义？

如下图，按照说明文档 配置好 Vscode 的 Wolfram Language Server之后，没有出现文档中的浮窗帮助和自动补全，请问可能还有什么问题呀，在vscode 里服务似乎也有报错

根据点集散点图生成了一个面，如何根据x、y轴求z轴的值呢？ date = \!\(\* TagBox[ RowBox[{"(", "", GridBox[{ {"1.`", RowBox[{"-", "10.`"}], "0.3`"}, {"1.`", RowBox[{"-", "5.`"}], "0.4`"}, {"1.`", "0.`", "0.6`"}, {"1.`", "5.`", "0.9`"}, {"1.`", "10.`", "1.2`"}, {"1.`", "15.`", "1.7`"}, {"1.`", "20.`", "2.3`"}, {"1.`",

具体情况可去官网了解。

如图所示。 这种中文提示实在没有必要，看着太碍眼了，从哪里设置可以去除呢？ （我用的是Mathematica 13版本）

大佬们，帮忙看一下，这个用mathematica来推导的话，是用Solve还是Eliminate，还是其他的。 在这里，Up1是已知量，还有剩下的电阻都是已知量，我的目标是通过这三个方程推出Us1、I1、Ip1、Is1的表达式（用前面的已知量表示）。

用mathematica11我这样作显示不出来 ContourPlot3D[3Sin[x^2+y^2+z^2]-2x*y*z==1,{x,-4.8,4.8},{y,-4.8,4.8},{z,-4.8,4.8},{x,y,z}∈Ball[{0,0,0},2],ContourStyle→{Opacity[0.8],Yellow},Mesh→None,Boxed→False,Axes→None,PlotPoints→30]//AbsoluteTimin 提示：ContourPlot3D……中位置 4 外应该是选项（而不是 {x,y,z}\[Element]Ball[{0,0,0},2]）. 选项必须是一个规则或者规则列表.

新人！想求解的问题是在热源项的激发下在金属薄膜中电子温度Te(z,t)与晶格温度Tl(z,t)随时间的变化。金属薄膜为100nm（100*10^-9m），时域为0-20ps(20*10^-12s)。方程组如下所示： 参数的话在这里我以金为例，具体见下，说一下初始条件与边界条件： 初始条件：Te(x,0)=Tl(x,0)=300K(室温） 边界条件： 参数：Ce(Te)=67.7*Te(z,t); ke=318; g=2.2*10^16; Cl=2.49*10^6; Qpump=0.30375*10^-9; rpump=15*10^-6; tp=2*10^-14s , =16.3*10^-9m 请大佬帮忙看看！

致的有趣问题。 “别随便拿内置函数/符号当变量名”——Mathematica教程里一般都会带一句类似的话，这句话是这么的理所应当，似乎没什么好多说的，可是，如果稍加思考，就会发现事情似乎没这么简单：为什么我们不能拿内置函数/符号当变量名？如果这么干了，会导致什么问题？ 比较容易想到的一种情况是：内置符号基本上都有Protected属性，因此诸如C=1的赋值默认情况下会失败。 可是除此之外呢？ 我标题里所说的，就是其中之一： 好了，废话了

我在想换元积分，两个积分分别是： NIntegrate[{(*换积分元*) g[z, rm] D[ z, {r, 1}]},(*积分上下限z=0\[Rule]r=rm,z=1\[Rule]r=\[Infinity]*){r, rm, \[Infinity]}] 和 NIntegrate[{(*换积分元*) g[z, rm]},(*积分上下限z=0\[Rule]r=rm,z=1\[Rule]r=\[Infinity]*){z, 0, 1}] 其中z = 1 - rm/r; 两个积分结果在a = 0时相同，但是趋势却完全不一样，求解答，以下是全部代码: ClearAll[ff, r, M, A, B, c, ddc, ddA, dc, \[Alpha], rm, g, z, r0, h, \ p, q, a, b, br, bd, d, um, R, hh, RR1, RR, rr1, rr, ff, dff, L] Clear[Subscript] Clear["Global`*"] \[

如图，新人解微分方程组时遇到两种报错，没明白报错原因，请问分别该怎么处理？ 方程组一： solve1 = NDSolve[{-3 x[t] y[t]^2 x'[t]^2 + 2 x[t]^2 y[t] x'[t] y'[t] - 4 y[t]^3 x'[t] y'[t] + 2 x[t]^3 y'[t]^2 - x[t] y[t]^2 y'[t]^2 + x[t]^2 y[t]^2 x''[t] + y[t]^4 x''[t] == 0, -x[t]^2 y[t] x'[t]^2 + 2 y[t]^3 x'[t]^2 - 4 x[t]^3 x'[t] y'[t] + 2 x[t] y[t]^2 x'[t] y'[t] - 3 x[t]^2 y[t] y'[t]^2 + x[t]^4 y''[t] + x[t]^2 y[t]^2 y''[t] == 0, x[0] == 1, x[1] == 0, y[0] == 0, y[1] == 1},

func1 = {D[u[x, y], {x, 2}] + \[Mu]*D[v[x, y], {x}, {y}] - \[Beta]^2* u[x, y] == 0, D[v[x, y], {y, 2}] + \[Mu]*D[u[x, y], {x}, {y}] - \[Beta]^2* v[x, y] == 0} cons = {u[0, y] == 0, v[x, 0] == 0, D[u[0, y], {y}] + D[v[0, y], {x}] == 0, D[u[x, 0], {y}] + D[v[x, 0], {x}] == 0，D[u[L/2, y], {x}] + \[Mu]* D[v[L/2, y], {y}] - (1 + \[Mu])*\[Alpha]*T == 0, D[v[x, B/2], {y}] + \[Mu]*D[u[x, B/2], {x}] - (1 + \[Mu])*\[Alpha]* T == 0} DSolveValue[{func1, cons}, {u[x, y], v[x, y]}, {x, y}] 求解之后没有反应，中间有一部分的边界条件没有显示 输出的红框部分与我的代码不一

NIntegrate去计算柯西主值积分;kz11[w]是已经过内插的在自变量w域1.5-25对应的函数；当例如计算x=1.6的时候，出现第一个提示，然后我根据指引用了MaxRecursion -> 24，虽然没有出现了提示，但例如又在计算x=2的时候，出现第二个提示，根据指引用了PrecisionGoal，然而又会出现新的提示，直到设置PrecisionGoal -> 1，但计算结果相差很大，应该如何解决。 主值积分背景：https://baike.baidu.com/item/柯西主值积分/18885835?fr=ge_ala NIntegrate[kz11[w]/(1.6 - w), {w, 1.6 + 10^-6, 25}]

1.kx11为经过内插（Interpolation）的函数； 2.http://note.youdao.com/share/?id=bf2a031f7d0d69064941032c7857dfd。这篇关于Compile的笔记当中，《原则三，模式匹配不能编译》还是不能很好理解，可以再解释一下吗 Jwij = Compile[{{x, _Real}}, NIntegrate[kx11[w]/(x - w), {w, 0.01, x - 10^-6}] + NIntegrate[kx11[w]/(x - w), {w, x + 10^-6, 25}], RuntimeAttributes -> {Listable}, Parallelization -> True]; Jwij[{0.011 , 2., 3., 6., 8., 10., 11., 12., 15., 16.}]

我的目的是想导入一个excel数据，但却发现了一个非常玄学的事情，前后两行代码明明完全一样，但就是有一个死活运行不出来。奇怪的是，如果我仅仅把第一行的路径覆盖到第二行里，第二行也还是运行不出来。但如果把路径两行的引号也一起复制和覆盖，那么第二行就可以运行了。有大佬知道是什么原因吗，感谢。代码如下： data1 = Import["C:\\Users\\DELL\\Desktop\\SC.xlsx"] data2 = Import["C:\\Users\\DELL\\Desktop\\SC.xlsx"] 这软件我真的是服了，每次

问题是这样的：已知 x 和 y 满足 (x - 2)^2 + y^2 <= 1 （即圆所围成的区域）。我想在这个圆围成的区域内画出满足不等式 Log[x] + x/2 > y + 2^y 的区域。 我想到的是用 RegionPlot 命令， 但是我查了MMA 的代码说明，好像 RegionPlot 只能在长方形区域内（即 a<=x<=b, c<=y<=d）求某个不等式代表的区域。但是我的问题是要在一个圆内画出一个不等式代表的区域，所以不知道该怎么写代码。 请高手指教，谢谢！

比如获得BesselJ[0，x]=0的所有x的具体数值

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

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