Reduce[{\[Alpha]*(( A^2 b - 2 A b (-1 + cr) \[Eta] + \[Eta] ((1 + b (-2 + f))^2 + b (-1 + cr)^2 \[Eta] - 4 b \[Beta]1^2 \[Theta]))/( 4 b \[Eta])) - (( A^2*b + \[Eta] + \[Eta]*b^2*(-2 + f) + \[Eta]^2* b*(-1 + cr)^2 + \[Beta]1^2*\[Theta]^2*b + 2 A*b*(\[Eta] - \[Eta]*cr + \[Beta]1*\[Theta]) + 2*\[Eta]* b*(-2 + f + \[Beta]1*\[Theta] - \[Beta]1*\[Theta]*cr))/( 16 \[Eta]*b)) >= 0; (1 - \[Alpha])*(( A^2 b - 2 A b (-1 + cr) \[Eta] + \[Eta] ((1 + b (-2 + f))^2 + b (-1 + cr)^2 \[Eta] - 4 b \[Beta]1^2 \[Theta]))/( 4 b \[Eta])) - (( A^2*b + \[Eta] + \[Eta]*b^2*(-2 + f) + \[Eta]^2*b*(-1 + cr)^2 + 9 \[Be

NDSolve怎么加入一个限制条件，使得函数值低于某特定值时停止计算

新手刚刚上手mathematica，在计算的时候发现矩阵求逆或者求特征值时，结果就和给的输入一样，没有计算，这是我哪里编错了还是个什么情况。 类似这几个图： 程序如下： nl = 13 Nl = EYmod*h Dl = EYmod*h^3/12 A2 = k22*G*h*(1 + h^2/(12*R^2)) A1 = Nl/R + Dl/(R^3) u = 0 lK2 = 12.3 kl2 = K2*Dl*b/(2*R^3) kl1 = kl2 kl3 = kl2 lk2 = lK2*Dl*b/(2*R^3) lk1 = lk2 y1 = 2*Pi/nl y2 = 0 y3 = y1 y0 = 0 nnd = 3 en = 1 m1c = 0 m2c = 0 m1s = 0 m2s = 0 cosm1 = 0 sinm1 = 0 cosm2 = 0 sinm2 = 0 m1c0 = 0 m1s0 = 0 m2c0 = 0 m2s0 = 0 Cn1 = -(nnd^2*A1 + A2/R) Cn2 =

如何用mathematica像图中那样，把一个3D球形矢量图（顶部图）展开到2维平面（底部图）

大佬们，能不能让y轴刻度-10^｛-6｝，-10^｛-5｝，-10^｛-4｝等间隔分布，我只能让正的等间隔分布，那负的怎么办呀，求助

例如：10^（-6）*sin（x），10^（-5）*sin（x），10^（-4）*sin（x），10^（-3）*sin（x），把这几个放到一个图中，而且要显示出振荡来，求助

soltest是得到的待定系数的解，其中第一个参数k的值先代入一个常微分非常求解，之后得到的函数和剩下的解代入到设的函数里面，我尝试一块代入，要不出错，要不有些没被替换 soltest = {{k -> c/q, a0 -> 0, a1 -> -1, b1 -> 0, c0 -> c/q, c1 -> 0, c2 -> 1, d1 -> 0, d2 -> 0}, {k -> (4 c)/q, a0 -> 0, a1 -> -(1/2), b1 -> 0, c0 -> (2 c)/q, c1 -> 0, c2 -> 1/2, d1 -> 0, d2 -> 0}}; odesol = Assuming[ {k /. soltest[[1]], c > 0, q > 0}, sol = DSolve[{y'[x] == k + y[x]^2, y[0] == 0}

比如我想定义 y[x_] := c^2 + x 其中c是关于x的一个表达式，c=x^2 实际上最后的表达式是y=x^4+x，只是通过两段

求大佬解惑 微分方程或者微分方程组中含有定义的已知函数（如一元二次函数）时，求解的结果存在问题 1、微分方程 In[1]:= g[x_] := x - x^2 In[2]:= DSolve[f[x] + g[x]*D[f[x], {x, 2}] == 0, f[x], x] Out[2]= {{f[x] -> -x C[1] Hypergeometric2F1[1/2 - Sqrt[5]/2, 1/2 + Sqrt[5]/2, 2, x] + C[2] MeijerG[{{}, {1/2 (3 - Sqrt[5]), 1/2 (3 + Sqrt[5])}}, {{0, 1}, {}}, x]}} 由于f[x]的两次导前面乘了g[x]结果中出现了很不常见的项（在资料库里说是两个函数），我不太清楚这是否算求得的结果（没接触过这样的通解） 2

我想画：一个图中，两条函数曲线，两条函数曲线用两个纵坐标分布显示。网上的方法好像不行，右侧的坐标只是个摆设，没有实际的用处

如何变为这种形式呢

其实经常在SE答题的盆友应该都有所察觉，那就是Mathematica对LaTeX的导入导出一直不完美。上面提及的这个问题很可能是继承自下面这个问题： ToExpression["\\mathbb{A}", TeXForm] 输出的是普通的A嗯。 \[DoubleStruckCapitalA]右键复制为LaTeX倒是可以正确导出为\mathbb{A}，不过这个步骤是由FrontEnd`CopyAsTeX[]函数控制的。参看 mathematica.stackexchange.com/q/102928/1871 个人感觉这回这个LaTeX直接输入功能的上线有点仓促，把 Underscript 的快捷键占了，结果 Underscript 的自带帮

(* 清除之前的定义 *)ClearAll[myClear, clc] (* 定义myClear命令用于清除全局变量的值，但保留变量名 *)myClear[] := (Unprotect["`*"]; Clear["`*"];) (* 定义clc命令用于清除当前笔记本中的所有单元格 *) clc[] := Module[{obj}, obj = EvaluationCell[]; SelectionMove[obj, All, Notebook]; NotebookDelete[]; "clc"]; 。。。。。。。。。。。。。。。。。。。。。。。。。。。。。 In[1]:= a = 1 Out[1]= 1 In[2]:= myClear[] In[3]:= a Out[3]= a In[4]:= a = 4 Out[4]= 4 In[5]:= myClear[] Out[5]= myClear[]

想要使用mathematica做一些数值求解偏微分方程组的东西，因此先随便使用一个方程去做测试。 方程如下：n[i]·Grad[u[i]]==Σ_{j} M[i][j]*u[j] n为单位矢量矩阵，我将其分量限制为全正数；M为复常数矩阵，将其设置为模值为1。i，j∈{1,2,...,num} 但有一些问题： 1，有结果不稳定的问题，这个原因不知为何，是否会影响正确性，我应该怎么检验其影响大小，又应该如何优化之，使此警报消失？ 2，运行的时候，CPU占用还比较多，但是我运行之时有时会突然停止运

我找到Mathprogramming-intro.org这个网站，但并没有在其中找到下载链接，该网站是已经关闭了吗？

这是一个已经报给官方并被证实的bug。在核心语言里发现bug的情况（相对）不太多见，姑且开帖说一下。总之，大家可以试一下下面的代码： ReplaceAt[Unevaluated@a[[0]], a -> {1}, 1] (*Language`ReplacementsDump`replace*)

用上新版了，高兴中。

一点偶然的发现，姑且记录一下： ClearAttributes[Plus, Listable] {1} + {1} {1} + {2} {1} + 1 SetAttributes[Plus, Listable] Times函数也有类似行为。当然了，这种行为不难模拟，增加一条类似f[a_List, a_List] := {f[a, a]}的规则就能实现。 最后，这个故事再次告诉我们，随便改内置函数是危险的，结果是不可预料的，能不这么干就不这么干。

并且这是没有任何通知或提示的。大家在发帖时千万注意。 不确定除外链外是否有其他内容会受此影响，大家若不放心，可在发帖后通过查看浏览器隐身模式加以检查。

CreateDialog[ Grid[{{B3, C1, C2, C3, C4, C5}, {C1, B3C1C1 = PopupMenu[1, {0, 1, 3, 5, 7, 9, 1/3, 1/5, 1/7, 1/9}], B3C1C2 = PopupMenu[1, {0,

如图，Reduce为何报错？请各位大佬指教

求出的fnpdf是两个瑞利分布和的概率密度，想用概率密度的特征积分为一进行验证，但是可求单点数值，一积分就得到非数值，代码如下 Clear["Global`*"] f[y_] = PDF[RayleighDistribution[0.5], y] te[t_] = (2*Pi)^0.5*FourierTransform[f[y], y, t] fnpdf[x_] = 1/(2*Pi)^0.5*InverseFourierTransform[te[t]^2, t, x](*N个特殊和概率密度*) N[fnpdf[0.5]] NIntegrate[fnpdf[x], {x, 0, Infinity}]

我尝试了下面的代码用 mma 来绘制函数图像 r = 2^64; Plot[r - Floor[r/b] (b - 1) + b, {b, 1, 30}, PlotRange -> All] mma 能成功绘制图像，此时执行其他代码也没有问题，但点击上述代码绘制的图像后 mma 会卡死甚至崩溃 想请问一下大家能复现这个问题吗？还是说这只是我系统上的问题 mma版本：14.0 中文版 + 本地文本包 系统版本：Windows 11

我用show命令画了三条曲线，并且给三条曲线标注了图例。但是发现两个问题：（1）图例之间的间距过大，请问用什么代码可以调整图例之间的间隔？（2）给图例设置字体大小的时候，发现对上标不起作用。比如 y^AB [x] 中 y 和 x 都变成了10号字体，但是上标 AB 的大小却没有改变，请问如何调整上标的字体大小？ 我写的代码如下： Show[Plot[{2^x + 1}, {x, 1, 3}, PlotLegends -> LineLegend[{TraditionalForm[(y^AB)[x]]}, LabelStyle -> 10], PlotStyle -> Red], Plot[{x + 3}, {x, 1/2, 3},

如图，ContourPlot3D的结果按理说至少有一条线吧

请问大家能否给一些初学者建议呢？我对函数变量都不太了解怎么打上去，应该如何学习呢？

就是不用根据网格线来涂色 比如 Plot3D[x^3+y^3-2x*y^2-2x^2*y+x^2-2y^2-10,{x,-10,10},{y,-10,10},RegionFunction→Function[{x,y,z}, (x^2)/100+(y^2)/81<=1.03],PlotStyle→Opacity[0.8],PlotRange→All,Axes→False,Boxed→False,BoxRatios→{1,1,1.2},PlotPoints→100,Mesh→6] 图像如下： 我需要 z ≥0 时作一种色；z ≤0 时作另一种色，请问可以怎么做？

"ABC" "XYZ" (*已知两个字符串可以拼接起来*)"ABC" <> "XYZ" Style["ABC", Red] Style["XYZ", Blue] (*那么两个带格式的字符串要怎么拼接？*)Style["ABC", Red] <> Style["XYZ", Blue] (*不能是 TableForm 这种网格形式的，他会有奇怪的对齐效果*)

Solve[Cosh[24 x]* Sin[24 x]*(Cosh[30 x]*Sin[24 x]*Sin[30 x]* Sinh[24 x] + (Cosh[24 x]*Sin[24 x]*Sin[30 x] - Sin[24 x + 30 x]*Sinh[24 x])* Sinh[30 x]) + (Sinh[ 24 x]*(2*Sin[24 x]*Sin[24 x] - Cosh[30 x]*(2*Cos[30 x]*Sin[24 x]*Sin[24 x] + Sin[30 x]*Sin[48 x]) + Cos[24 x]*Cos[24 x]*Sin[30 x]*Sinh[30 x] + Cos[30 x]*Sin[48 x]*Sinh[30 x]) + Cosh[24 x]* Sin[24 x]*(Cosh[30 x]*Sin[30 x]*Sin[24 x] - Sin[30 x + 24 x]*Sinh[30 x]))*Sinh[24 x] == 0, x]

做了一个MMA的视频教学，争取一周一更新，但不保证。视频上传到腾讯视频了，优酷有写地方不好，所以被我放弃了。 我的腾讯视频主页：

请问各位大佬，这俩积分咋算呀

我发现Solve只能求解单个变量如x,或者两个{x,y}，那么如何求解一个表达式呢（比如x*y） 对于结果是一个常量的表达式我通过FullSimplify来求解，如： 但是如果结果也是一个代数式呢，我发现文档的Reduce，Refine函数也都是只能针对变量进行操作，那么如何才能针对一个代数式进行求解呢，比如： 已知 求： 我发现如果直接使用Solve还是FullSimplify都没法解决，我目前只能手动的去一步步做： verticalEnergy := v1 Sin[theta1] - v2 Sin[theta2] totalEnergy := (m v1^2 - m v2^2)/2 + (

版本是12.2，被积函数是x^2*Sin[a*x+b*x^2+c*x^3]这样的形式，是一个震荡的函数，所以可能数值积分上确实会存在些问题。 我和另一个同学分别写的，我用的是NIntegrate，他用的是Integrate，剩下的都一样。然后发现计算结果不一致。 首先是他用Integrate计算的结果，带有一个相对实部很小的虚部；而且实部和NIntegrate也差几个量级。 然后就是我为了验证这个问题，在同一个notebook里，先后用Integrate和NIntegrate计算，发现和我们分别算的结果是一样的，依然存在着

为选项OutputSizeLimit指定的设置无法被使用.

In[12]:= v[t] = D[x[t], t] Out[12]= Derivative[1][x][t] In[19]:= DSolve[(m g + a x[t] g - a v[t] v[t]) == (m + a x[t]) D[v[t], t], t, x[t]] \:6B63\:5728\:8BA1\:7B97In[19]:= DSolve::dvnoarg: 函数 t 出现时没有参数. Out[19]= DSolve[ g m + a g x[t] - a Derivative[1][x][t]^2 == (m + a x[t]) (x^\[Prime]\[Prime])[t], t, x[t]]

如图所示，在输出的结果中有小括号代表什么意思呢？能否输出一个确定的值呢？源代码如下所示： In[1]:= n = 2;(*稳定器个数*) D1 = 24.4;(*钻头直径Db*) D2 = 24.4; (*井眼直径*) l1 = 1800; (*第一跨钻铤长度*) l2 = 900; (*第二跨钻铤长度*) Dl1 = 20.32; (*第一跨钻铤外径 8寸钻铤*) DL1IN = 7.14;(*第一跨钻铤内径 8寸钻铤*) Dl2 = 17.8;(*第二跨钻铤外径 7寸钻铤*) DL2IN = 7.14;(*第二跨钻铤内径 7寸钻铤*) Dl3 = 17.8;(*第三跨钻铤外径 7寸钻铤*) DL3IN = 7.14;(*第三跨钻铤内径 7寸钻铤*) DS1 =

我想用Mathematica复现一篇论文中的结论，解超定方程组，问题和结论如图,超定方程组我放到最后了 我的代码 (* 这部分代码用于推导方程式，所以你可以忽略它。 *) eqn = c^2*(n - 1)^2*S^3 + a*(n - 1)^2*S^2 + b*n*(2*n - 1)*(1 - Y^2)^2*(DS)^2 - b*n*(n - 1)*S*(1 - Y^2)*(-2*Y*DS + (1 - Y^2)*D2S); S[Y_] := a0 + a1*Y + a2*Y^2; DS = D[S[Y], Y]; D2S = D[DS, Y]; (* 这里，我们可以得到上述提到的公式 *) seqn = eqn /. {S -> a0 + a1*Y + a2*Y^2, DS -> DS, D2S -> D2S} (* 得到超定方程组 *) coeList = CoefficientList[Expand[seqn]

我想求解这样一个关于x1,x2的二元一次方程，可以把mx,ny,theta都认为是参数 Solve[mx == x1 Cos[\[Theta]/2] + x2 Cos[\[Pi]/3 - \[Theta]/2] && my == x2 Sin[\[Pi]/3 - \[Theta]/2] - x1 Sin[\[Theta]/2], {x1, x2}] 如下图所示，如果我单独开一个文档，每次都是直接打开然后开始算，就可以解出来结果。 但有时候我不是这样单独开一个文档，而是放在我的其他脚本中，就算上一步才清理完全部变量，也会报错说 Solve::incs: Warning: Solve was unable to prove that the solution set found is complete. 研究

大键盘中Enter是换行，可是小键盘中的Enter居然是执行，怎么把小键盘的Enter也改成换行啊！！！！！！

无论是内置的Eigenvalues，还是使用SchurDecomposition，总会莫名的改变原矩阵特征值的顺序。比如一个4阶方阵，test ={{600, 20, 20, 20}, {20, 400, 20, 20}, {20, 20, 800, 20}, {20, 20, 20,1000}}，对角原远大于非对角元，显然特征值的顺序应该是600左右的一个，400左右的一个，800左右的一个，1000左右的一个。可是Eigenvalues一定会从大到小牌序，用Schur分解的方法就更离谱了： test 1=SchurDecomposition[N@test][[1]]; Transpose@test1.test.test1 给出的顺序是400,600,1000,800，并且还会随着非对角元

已知有关5个变量的3个隐式方程：f1(x, y, z, a, b), f2(x, y, z, a, b), 以及f3(x, y, z, a, b)。我们可以得到xyz坐标系下的一个曲面S。该如何来画S呢？ ContourPlot3D不行，因为它只能接受一个方程作为输入，这里有3个；ParametricPlot3D也不行，因为它要求输入的方程为显式表达，这里为隐式；我也想通过直接Solve的方式把隐式变为显式，但方程有点复杂也做不到。求求各位大神教教我！有没有什么办法可以画出S？

请问下大伙儿们，Mathematica的NDSolve命令有提供梯形法求解吗？或者其他求解微分方程的命令有没有提供梯形法的？

Abs[lbk]-4 \[lbk]Pi[rbk]^2 - 4 I \[lbk]Pi[rbk] x + x^2[rbk] Assuming[lbk]x \[lbk]Element[rbk] Reals, Simplify[lbk]Abs[lbk]-4 Pi^2 - 4 I Pi x + x^2[rbk][rbk][rbk] 这是我尝试的两个公式，但是都无法直接计算虚数的模，请问加什么限制条件可以化简成最下面的结果呢

In[27]:= Maximize[Sin[2 x - 11 \[Pi]/6], x] Out[27]= {1, {x -> -((143 \[Pi])/6)}} 我希望他输出带C的所有解析解，但是这个命令只给我了其中一个x所对应的值。 用偏导的话感觉很麻烦要一个一个排除，有没有命令可以一次性解决所有的极（大，小）值哇 （不用D【f】==0）

我的目标是用mma帮忙推公式，减去手算的麻烦。目标公式如下，论文里随便找的。 这是我一顿乱写的代码，得出结果是一致的。想问问特别最后一步，单对分母按需求进行合并同类项，有没有更简便的写法。 代码： << Notation` Notation[ NotationTemplateTag[x_ || y_] \[DoubleLongLeftRightArrow] NotationTemplateTag[((x_ y_)/(x_ + y_))]](*定义并联算符*) Subscript[g, m] (1/(j w Subscript[m, 2 ] Subscript[C, gd]) || 1/( j w Subscript[C, ds]) || Subscript[R, L] || Subscript[R, ds]) // Simplify ((Subscript[g, m] Subscr

如图，看了帮助，但是实在找不到和这种形式类似的例子了，希望大佬帮助，万分感谢！