考虑到此类问题和“Mathematica已经装好了，但不能用”即“Mathematica一开始就不能用”有一定区别，特开专帖讨论。今后此类问题原则上请不要单开主题。 这里的“不能用”，包括不限于 1. Mathematica突然无法正常打开 2. Mathematica突然变得无法正常工作（1+1也算不了等） 3. Mathematica突然在执行正确的代码时出现莫名其妙的警告。 如果你的“Mathematica一开始就不能用”，请看置顶。 闲话少叙，说一说可能的修复方法： 此类问题似乎和Mathematica的Paclet更新不

已知 AB^2 == b^2 + h^2 AC^2 == c^2 + h^2 AD^2 == d^2 + h^2 BC^2 == b^2 + c^2 - 2 bc Cos[lbk][lbk]Alpha[rbk][rbk] CD^2 == c^2 + d^2 - 2 cd Cos[lbk][lbk]Beta[rbk][rbk] BD^2 == b^2 + d^2 + 2 bd Cos[lbk][lbk]Alpha[rbk] + [lbk]Beta[rbk][rbk] 我想用AB,AC,AD,BC,BD,CD来表示h，应该怎么做（本人刚用mathematica，希望各位大佬帮忙指导一下，谢谢！）

版本为12.0，Integrate函数给出错误的复数结果，但NIntegrate函数结果正确。 R[m_, n_] := Integrate[(1 - Cos[m x + n y])/(2 - Cos[x] - Cos[y]), {x, -Pi, Pi}, {y, -Pi, Pi}]/(4 Pi^2) (*R[1,0]*) Integrate[(1 - Cos[x])/(2 - Cos[x] - Cos[y]), {x, -Pi, Pi}, {y, -Pi, Pi}]/(4 Pi^2) (*(2 I Log[2])/\[Pi]*) NIntegrate[(1 - Cos[x])/(2 - Cos[x] - Cos[y]), {x, -Pi, Pi}, {y, -Pi, Pi}]/(4 Pi^2) (*0.5*)

NDSolve:无法求得一个导数公式，考虑使用选项设置：...... 大致看一下这是什么问题？我在知乎上问说没有见过。 这个代码有些长，我在评论区也补上了一部分。 l = 1 g = 2 w = 2 b = 2 a = 3 \[Theta] = 2 m = 1 L10 = 1/6 (-6 g w (l Cos[\[Theta]] Cos[\[Alpha]1[ t]] Cos[\[Theta] - \[Beta]1[t]] + Sin[\[Theta]] (b + l Sin[\[Theta] - \[Beta]1[t]])) - 3 m (3 b Sin[\[Theta]] + a Cos[\[Theta]] Sin[\[Alpha]1[t]] - 2 l Cos[\[Beta]1[t]] Sqrt[ Cos[\[Alpha]1[t]]^2 Cos[\[Theta] - \[Beta]1[t]]^2 + Sin[\[Theta] - \[Beta]1[t]]^2]) - 6 g w (Cos[\[Theta]] (-l Cos

也算是对这帖（http://tieba.baidu.com/p/2964416898）的2，3，4条的解说吧。有什么意见或建议欢迎提出： http://note.youdao.com/share/?id=abd51087f44c0b6a41ff6022d549dc4

症状： 顶楼中的 [ 和 ] 会被自动转换为 [lbk] 和 [rbk]。此现象在回复帖中不存在。 影响范围： 此bug至少在去年10月16日即已出现，至今（安卓12.58.1.0）未修。IPhone受影响情况不明，欢迎实验补充。 对策： 对于提问者，最简单的对策就是改用电脑访问贴吧网页版，不要使用使用手机客户端发帖。如果确有相关需求（所以说你们到底为什么会有这种需求？电脑没连网？那也可以开手机热点啊，访问网页又花不了什么流量），避免在顶楼张贴代码也是一个

kb = 1.38*10^-23; T = 300;(*charge*) muG = 0.1; muT = 2*kb*T*Log[2*Cosh[muG/(2*kb*T)]]

先贴代码 dm3=If[d>c,If[0<=(d-c-2)/(d-2)<=1,1-(d-c-2)/(d-2),Infinity],Infinity]; DM3=Table[dm3,{c,0,8,1},{d,0,50,1}]; ListDensityPlot[DM3,Mesh->Full] 这里我明明取值的时候是按间隔1地取的，出来的不应该是一个个色块吗为什么在每个单位1的格子里颜色还有深浅过渡变化，看起来变成了解析的DensityPlot如果想让它每个色块里都是单一深浅的颜色该怎么办呢

下面这个积分，我想求解析解，计算了两个小时还没有出结果来，一直显示在运行。不懂就问，这是什么原因，有可能是算力不够吗？cpu核心全开，新款笔记本机子。如图。 q*k*Q/L/f*Simplify[lbk] Integrate[lbk] Integrate[lbk]((l1*Sin (\[lbk]Theta[rbk]) + d1)^2 + (l1*Cos (\[lbk]Theta[rbk]) + y)^2 + (l2*Sin (\[lbk]Alpha[rbk]) - z)^2)^(-3/2)*(-l1* Cos (\[lbk]Theta[rbk]) - y), {y, L1 - L, L1}[rbk], {z, -f1, f - f1}[rbk], Assumptions -> {Element[lbk]f, Reals[rbk], Element[lbk]f1, Reals[rbk]}[rbk]

链接不发了。楼下简述更新内容

本文系自此帖（http://tieba.baidu.com/p/2736346260）5楼的回复修订增补而来。吧里的老人们阅读本文时可能会觉得各种眼熟：毕竟重要的经验我都是不断强调的。 我最初会得知Mathematica的存在，是因为《数学建模》老师在课堂上的推荐。我第一眼就被它那极度接近传统数学式书写的使用方法给吸引了——相信这个吧里的许多人也是一样。《数学建模》结课后不久校内就搞了数学建模大赛（目的是为了选拔参加全国大赛的选手），我就参加了，用的当然是Mathem

此帖是本吧初代吧主 @mm_酱 发在果壳小组的一篇教程。如今果壳小组已死，mm_酱近期又不太活跃，这里姑且重发一份： https://note.youdao.com/share/?id=058e6037396d925af1f4abe4d54a52 mm_酱你要是想整理下出个v2就吱一声。

https://1drv.ms/u/s!Amd-uedNeUsttDBBDZWPGllTqawk?e=IohgLF

贴吧新弄了个贴条，贴条只能同时存在10个并且只能实际只能显示26字，于是我就在想着是不是该总结些Mathematica新手最常见的问题贴到那上面滚动播出（虽然很多人或许不会看，但有总比没有好嘛……）。试着弄了个草稿（二十戒），大家可以看看，哪些需要修改，哪些需要删节或补充： 二十戒 1　不要为了节约硬盘选择过老的版本。至少要用版本7 2 　自带帮助是最好的教材，提问前先打开软件按下F1仔细查查 3　将光标停在不认识的函数前/中/后再

拉格朗日微分方程组(常系数二阶，四元)一直显示“正在运行”，已经2天了，看内核计算状态确实也是一直在算，电脑一直在出声，前面也没有报错。难道这个只能通过用数值解来慢慢搞出来了吗？对于这样的拉格朗日微分方程组数值解大致怎么用？ \[lbk]Theta[rbk] =ArcSin (-(y[lbk]t[rbk] + l*(Cos (\[lbk]Alpha[rbk][lbk]t[rbk]) - Cos (\[lbk]Beta[rbk][lbk]t[rbk])))/L)d\[lbk]Theta[rbk]dt = (l*(Sin (\[lbk]Alpha[rbk][lbk]t[rbk])*\[lbk]Alpha[rbk]'[lbk]t[rbk] - Sin (\[lbk]Beta[rbk][lbk]t[rbk])*\[lbk]Beta[rbk]'[lbk]t[rbk

嗯……这个问题的复杂度我心里还是有个数的，所以这帖不指望得出答案之类，重在交流讨论。 先说一下我大致的了解吧。因为最近有机会接触到了多核的计算机，所以我也稍微关心了一下。总的感觉是并行似乎真不那么容易，在Stackexchange的未解答率（无人回答或有人答却未被采纳或没有Upvote的问题数比上总问题数）只有35/175.==0.2，这算是相当高了，要知道全站总体未解答率是1325/14474.==0.0915435，一些流行的Tag，比如plot，未解答率只有63./2069==0.0304495

但是，请做好没人答你的问题的心理准备。 决定暂时不做处罚的原因，主要是目前会这么做的人很少。（一个月也不见得有一个。）“司法”实践不足。 顺便，自己没有做任何尝试的帖子原则上都可以判断为求码帖，诸位吧务应该没忘吧？ 最后表个态，我不会回答任何设置了“你需要关注TA才能发表评论”的主题。

调节自变量，控制某个目标量接近最大值最小值或设定值？

我暂时无法将这个错误缩减为一段简单的代码, 完全的代码在二楼. 问题完全出在被命名为problem的函数上, 所以只看problem这一块就好了. problem函数计算了两个积分w0和w1, 第一个积分w1在数学上的结果是1. 尝试运行测试函数 test[1] 第一个积分结果为1, 很好, test[2] 第一个积分结果就不正确了. 但! 如果退出内核, 先运行 test[2] 第一个积分结果为1, 反而后运行的 test[1] 不正确了, 也就是只能正确计算首次接受的参数. 错误可能是Δ[index_] := Δ[index] = … 这样的“有

代码：Integrate[E^(-I*(m*x + n*y))*(x+y)*E^(-Sqrt[x^2 + y^2]), {x, -Infinity,Infinity}, {y, -Infinity,Infinity}] 不是很复杂，但是软件得不到结果，不知咋回事： 这个确定是能积的，结果类似于以下解析式：

第一张图是求解模型所用的方程组，我想得到在那个矩形区域内的温度分布图（ContourPlot）。 第二张图是我在mma上的代码，但是它提示我说方程个数大于了未知量数，删掉1个条件又出现错误，求同志们帮忙。 代码在最下方（可复制）。 谢谢。 a = 0.5*10^-3; b = 20*10^-3; \[Rho] = 11.3*10^3; c = 130; \[Lambda] = 300; i = 10; R = 0.01; H = 15; fun = \[Rho]*c*D[T[r, z, t], t] - \[Lambda]* Laplacian[T[r, z, t], {r, z}] - (i^2 R)/(\[Pi]*a^2*b); ic = T[r, z, 0] == 300.15; bc1 = fun == NeumannValue[H*(T[r, z, t] - T[r, z,

请教一下大佬，有谁知道PrimeQ用的什么算法判定质数？结果是100%正确吗？

Mathematica在本地有个质数库嘛？ Prime[n] 取得很大的质数，速度非常快。 PrimeQ[x] 判断质数速度也很快。 我用python，java写代码判断效率和它差了十条街。 难道它的算法这么优秀？ 怀疑它是不是本地有个质数表，存了一些信息，然后在这个基础上在计算出来的？

现在我有一个Table{x,y}，大概几千个点，怎么直接交换两个坐标的值生成新的Table啊，就把这几千个数据的{x,y}变成{y,x}，F1没找到相关的指令

f[x_ ] := x^2 + 1; g = Compile[{x }, x^2 + 1]; h = #^2 + 1 &; a = f /@ Range[10^5]; b = g /@ Range[10^5]; c = h /@ Range[10^5]; a == b == c ByteCount /@ {a, b, c} (* True {2400080, 800200, 800200} *) 表面上a,b,c结果是相同的，但如果需要将a,b,c用于后续的计算，a与b,c计算时间会有明显的差别 Min[a]; // AbsoluteTiming Min[b]; // AbsoluteTiming Min[c]; // AbsoluteTiming (* {0.0005798, Null} {0.0000589, Null} {0.0000559, Null} *) 这是什么原因

已经把公式中其他内容带入数值，运行不出结果，而且单独计算里面的积分也运行不出来 DSolve[P[\[Omega]] == \!\( \*SubsuperscriptBox[\(\[Integral]\), \(0\), \(Pi/2 - 10^\((\(-3\))\)\)]\(\((\(-\((\((66416.66137691677`\ \ \((\(-3.6551383444493544`*^14\) + \*SuperscriptBox[\(E\), \(38.683223086147535`\ \*SqrtBox[\(0.7514187629` + 0.005167396856733549`\ Cos[\[Theta]]\)]\)])\)\ \ \((\(-2\)\ \(( \*FractionBox[\(9\), \(1600\)] + \*FractionBox[\(5.618`*^21\), SuperscriptBox[\(\[Omega]\), \(2\)]])\)\ \*SuperscriptBox[\(Cos[\[Theta]]\), \(2\)] + 3.0056750516`\ Cos[ 2\ \[T

首先两台设备的软件安装和使用都是我操作的，确认为相同版本。 有区别的状况是在我的个人pc上发生了报错，在远程使用的工作站上没有报错，因为要调试所以两边的设备上都有试运行。但最终输出的结果一致。这种图是报错类型，代码我放在后面。请问这种情况下结果精度问题大吗？

最近在重看有限元相关的东西，结果在看到等参元这部分内容的时候发现了一个数学上的小问题，姑且拿出来水一帖。考虑下面这个参数方程： base = DeleteCases[{-(1/4) (-1 + s) (-1 + t) (1 + s + t), -(1/ 4) (1 + s) (-1 + s - t) (-1 + t), 1/4 (1 + s) (1 + t) (-1 + s + t), 1/4 (-1 + s) (1 + s - t) (1 + t), 1/2 (-1 + s^2) (-1 + t), -(1/2) (1 + s) (-1 + t^2), -(1/2) (-1 + s^2) (1 + t), 1/2 (-1 + s) (-1 + t^2)} /. t -> -1, 0][[{1, 3, 2}]] (* {1/2 (-1 + s) s, 1 - s^2, 1/2 s (1 + s)} *) 这个方程有如下性质：s为-1时，第1项为1，其余

rt，楼主发现在做指数带变量的解析计算的时候，mma给出的最终结果指数部分都没能化简，甚至连x^y·x^(1-y)这样的式子也不会化简成x，有什么办法解决吗

我有一个这样的积分： Integrate[(-t Cos[t] + Sin[t])^2/t^5*Erf[2*a*t], {t, 0, +\[Infinity]}, Assumptions -> a > 0] 在mma 13.3版中，上述表达式能给出正确积分结果： 我最近将mma升级到了14.0版，结果同样的式子计算不出来了，直接把输入原封不动地输出： 有大佬知道这个是什么原因吗？我查了13.3到14.0的更新日志，里面没有显示对这个Integrate函数进行了更新。以及有没有除了回退旧版本以外的解决方法？万分感谢

需要绘制P'/P为函数，Pt/P为自变量的函数图。看到函数后无法下手，在此求助各位大佬了！

请问这种图mathematica怎么画，不同区域有不同的颜色和填充 是用的纹理吗？ 如果只用网格和regionplot能画吗，就是不同区域有不同的颜色和形状填充

函数为 aM[k_, \[Tau]_] := (2/ NSolve[w - 1 + k*Sin[w*\[Tau]] == 0, w, Reals][[1, 1, 2]])* Cos[NSolve[w - 1 + k*Sin[w*\[Tau]] == 0, w, Reals][[1, 1, 2]]*\[Tau]]* Sin[NSolve[w - 1 + k*Sin[w*\[Tau]] == 0, w, Reals][[1, 1, 2]]/2] 但是将其中Nsolve部分单独拿出来设定参数明明是有解的

Reduce[-(( 4 (-1 + r) \[Alpha]^2 \[Beta]^2)/(-4 \[Alpha] \[Beta] + r (\[Beta] (-1 + \[Theta])^2 + \[Alpha] \[Theta]^2))^2) + ( 4 (-1 + r) \[Alpha]^2)/(-4 \[Alpha] + r (-1 + \[Theta])^2)^2 >= 0 && 0 < \[Theta] < 1 && 0 < r < 1 && 0 < \[Alpha] && 0 r < 4 \[Alpha] < 4 \[Beta] && -4 \[Alpha] \[Beta] + r (\[Beta] (-1 + \[Theta])^2 + \[Alpha] \[Theta]^2) < 0 && -4 \[Alpha] + r (-1 + \[Theta])^2 < 0, \[Theta]] // FullSimplify

孩子之前是windows，现在换mac了，一打开软件懵了，界面太干净了，想问问是不是我安装时哪里出问题了

ClearAll 的正确用途及语法请参看帮助，不知道怎么查帮助的人，以及不知道我现在在说啥的人，请参看《你查过自带帮助了吗？你真的会查自带帮助吗？》：https://tieba.baidu.com/p/4392855072 顺便，clc，clear all之类的东西在Mathematica里也是没有任何意义的。 再顺便，如果你模模糊糊地觉得MATLAB、Mathematica、Maple、C、Pascal……其实是一个东西，那么，你是错的。 （这条其实已经收在《十戒》里了，但想想还是单开个标题好懂的帖强调一下。）

我的代码如下： ClearAll; n = 2; Bernstein[i_, j_, k_, u_, v_, w_] := Factorial[ n]/(Factorial[i] Factorial[j] Factorial[k] Factorial[ n - i - j - k]) u^i v^j w^k (1 - u - v - w)^(n - i - j - k); P = Array[x, {n + 1, n + 1, n + 1, 3}]; deltaP = Array[deltax, {n + 1, n + 1, n + 1, 3}]; BezierTet[u_, v_, w_, t_] := Sum[Bernstein[i, j, k, u, v, w] (P[[i + 1, j + 1, k + 1]] + deltaP[[i + 1, j + 1, k + 1]] t), {i, 0, n}, {j, 0, n - i}, {k, 0, n - i - j}]; BTJacobi[u_, v_, w_, t_] := Dot[Cross[D[BezierTet[u, v, w, t], u], D[BezierTet[u, v, w, t], v]], D[BezierTet[u, v, w, t], w]]; Collect[BTJa

Solve[lbk]{40*492 \!\(\*SuperscriptBox[lbk]\(x\), \(\[lbk]And[rbk]\)[rbk]\) 4 - (492*238145 + 492*16021 + 40*16021 + 142*2800*Cos[lbk]xy[rbk]) \!\(\*SuperscriptBox[lbk]\(x\), \(\[lbk]And[rbk]\)[rbk]\) 2 + (-40*2800 \!\(\*SuperscriptBox[lbk]\(x\), \(\[lbk]And[rbk]\)[rbk]\) 3 + 238145*2800 x)* Sin[lbk]xy[rbk] + 238145*16021 == 0 && -40*142 \!\(\*SuperscriptBox[lbk]\(x\), \(\[lbk]And[rbk]\)[rbk]\) 3 - 492*142 \!\(\*SuperscriptBox[lbk]\(x\), \(\[lbk]And[rbk]\)[rbk]\) 3 + 142*2800 \!\(\*SuperscriptBox[lbk]\(x\), \(\[lbk]And[rbk]\)[rbk]\) 2*Sin[lbk]xy[rbk] + 238145*142 x + (-40*2800 \!\(\*Su

Manipulate[lbk]Plot[lbk]aSin[lbk]bx + c[rbk], {x, -6, 6}, PlotRange -> {-5, 5}[rbk], {a, 1, 5, 0.2}, {b, 1, 5, 0.2}, {c, 10, 2 Pi, Pi/100}[rbk]

开这帖没别的目的，就为了把十戒（http://tieba.baidu.com/p/2964416898）里的9到11条 9　没赋值的变量是蓝色的，注意这点能有效避免低级错误 10　赋过值的变量会变黑，注意这点能有效避免低级错误 11　内置函数全是黑色的，注意这点能有效避免低级错误 以一个更显眼的方式强调一遍（话说能理解上面三句的诸位可以出帖了，下面没新东西），顺便补充少许的例子。 为了便于修订，照例开了个笔记，请大家多点一下鼠标： note.youdao.com/noteshare?id=06adaac6aaf85c8

一开始以为是太复杂，因为我知道代码整的比较粗糙计算比较笨。但是少一个变量的大约6×10×100个类似的式子算得还是比较快的，而且换到工作站用96核上还是跑不动。 另外分这么多变量事实上是另一步定积分的数值计算我给手动拆区间取点取均值了，我目前只是单纯想要相关计算结果，看了下精度和之前耗时还可以接受所以暂时先将就用着在，准备后面再慢慢看怎么做合理。

我是做理论建模的，会用到很多符号计算，因此常常会遇到复杂函数化简的问题，还请各位大神讲解！

Integrate[(a b - b x Cos[\[Phi]] - a y Sin[\[Phi]])/ Abs[z^2 + (x - a Cos[\[Phi]])^2 + (y - b Sin[\[Phi]])^2]^(3/ 2), {\[Phi], 0, 2 \[Pi]}] 为什么求不出具体结果呢

一个数用两种类型表示算出的结果为何差距这么大： N[Integrate[ Exp[(10^(-4))*t]*ArcTan[((10^(-4)) - t)^2], {t, 0, 10^(-4)}]/ Integrate[2*t*((10^(-4)) - t), {t, 0, 10^(-4)}], 10] 1.0000000025+ 0.*10^-11i N[Integrate[Exp[0.0001*t]*ArcTan[(0.0001 -t)^2], {t, 0, 0.0001}]/ Integrate[2*t*(0.0001 - t), {t, 0, 0.0001}], 10] -32.7418 请教各位这是怎么回事啊

他的帖子也跟着隐藏了，目前吧里损失不小，希望有空的朋友能一起来帮忙申诉下。具体地说，就是访问 tieba.baidu.com/pmc/feedback 然后选“账号被封禁”，把他的账号 asdasd1dsadsa 贴进去，再回车，如下图：

f[lbk]\[lbk]Omega[rbk]_, n_[rbk] := Module[lbk]{m = 1, m1 = 1, m2 = 1, k1 = 1, k2 = 0.105, k3 = 1, \[lbk]Alpha[rbk]1 = 0.5, \[lbk]Alpha[rbk]2 = 0.02, \[lbk]Alpha[rbk]3 = 0.4, p = 1, \[lbk]Zeta[rbk] = 0.005, \[lbk]Omega[rbk]n1 = 1.`, \[lbk]Omega[rbk]n2 = 1.1`, \[lbk]Phi[rbk]11 = 0.7071067811865476`, \[lbk]Phi[rbk]12 = 0.7071067811865468`, \[lbk]Phi[rbk]21 = 0.7071067811865476`, \[lbk]Phi[rbk]22 = -0.7071067811865482`, A1 = 0.2250000000000001`, A2 = 0.07499999999999998`, A3 = 0.6749999999999999`, A4 = 0.02499999999999962`, A5 = 0.02499999999999998`, A6 = 0.675`, A7 = 0.07499999999999862`, A8 = 0

(\[Kappa] = 1)*；* 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*Sqrt[2]*\[Eta]*\[Kappa]^3*\[Lambda]