y = Integrate[Exp[x/2 (1 + x^2 - 2 x)], x]; Table[{x, y}, {x, 0, 1, 0.2}] 如例子所示，想求这个函数y(x)的数值Table，但是y是x的一个积分形式没法解析求解，我想生成一个{x,y}的table画它的函数图像但是用Integrate的话它不能直接取值会变成图里这个样子，求问一下怎样才能让它积分后再取值啊用NIntegrate的话它又不支持不定积分一定要输入上下限，如果能直接y=NIntegrate[y(x),x]就好了

我的数据点都是0到1.2跟0到5按0.01间隔进行采样的，用ListDensityPlot画出来的图它的坐标就放大了100倍，有没有什么方法把它改回来？用DensityPlot的话是正常的但是精度又不够，帮助里没找到有调这个的例子

题一： 已知 f(x)+f(1/x)+3f(-x)=x , 求f(x)的解析式 题二： 已知 f(f(x))=x^2-3x+4 , 求f(2)的值和f(x)的解析式

代码和运行结果如下： 在吧内尝试着找了一下没有找到相关的帖子，我自己觉得应该是与替换规则的使用或者代码的运算顺序有关系，但是不知道应该去找哪些资料来看，所以如果存在现有的资料的话就拜托大家帮忙指个路了 mKc = 0.493677; mpi0 = 0.1349768; f = 0.093; \[Epsilon] = 0.001; omega[m_, q_] := Sqrt[m^2 + q^2]; fcut[p_, cut_] := Exp[-(p/cut)^2]; g[s_, cut_] := NIntegrate[((4*Pi)/(2*Pi)^3)*((omega[mKc, q] + omega[mpi0, q])/(2* omega[mKc, q]*omega[mpi0, q]))* ((q^2* fcut[q, cut])/(s - (omega[mKc, q] + omega[mpi0

我尝试了Reduce，Simplify等函数，没做出来，看看有没有高手指点一下，谢谢大家！

才发现ListContourPlot的“描边”颜色总比设置的颜色淡，会自动加了层透明度（大概Opacity[0.5]）😂，搞得和别的图总有个色差。 也可能是我火星了，搜了圈google也没有太相关信息。 也没有默认设置能调，帮助里好像也没看到提这个事，不过用Opacity能改，就是有点麻烦

Subscript[\[Pi], m]=(2 a c (t+\[Beta])-c^2 (t+\[Beta])^2+2 t T (-4+2 \[Alpha]^2+(t+\[Beta])^2 \[Gamma]^2)+((-1+\[Theta]) Overscript[a, _]-\[Theta] Q) (2 a+(-1+\[Theta]) Overscript[a, _]-\[Theta] Q))/(-8+4 \[Alpha]^2+2 (t+\[Beta])^2 \[Gamma]^2); Subscript[\[Pi], r]=-(((-2+\[Alpha]^2) (2 a-c (t+\[Beta]))+a (t+\[Beta])^2 \[Gamma]^2+(-2+\[Alpha]^2+(t+\[Beta])^2 \[Gamma]^2) ((-1+\[Theta]) Overscript[a, _]-\[Theta] Q))^2/(2 (-2+\[Alpha]^2) (-4+2 \[Alpha]^2+(t+\[Beta])^2 \[Gamma]^2)^2));A=Q*\[Theta]+(1-\[Theta])*Overscript[a, _];optw=((-2+\[Alpha]^2) (A+c (t-\[Beta]))+A t (t+\[Beta]) \[Gamma]^2)/(-4+

如图所示，提示我Failed to load some of the packages(see logs) 某些包加载失败？叫我看日志，我也不知道日志在哪...... 我装的是offline版本 wolfram engine我已经装好了，启动正常

我用的是14.0中文版的Mathematica，在打开调试打开堆栈窗口的时候会跳出很多堆栈窗口，直到再跳出“格式化笔记本”的对话框。运气好点打开十八九个窗口后会停下来，其中只有一个窗口能显示堆栈，其它都是空白，内存也是占用了一两个G，运气不好的话停下来后直接完全退出，对文件的改动不能保存。 原来用的是14.0是注册机弄的，现在换成试用版下载安装也是这样，软件路径也确定是英文，但是软件本地的存档没有清过。这个情况常见吗？是什么

假设我有一个关系图，以整数之间的因数关系举例： ns = {1, 2, 3, 4, 5, 6, 7, 8, 9, 10}; as = {1 \[DirectedEdge] 2, 1 \[DirectedEdge] 3, 1 \[DirectedEdge] 4, 1 \[DirectedEdge] 5, 1 \[DirectedEdge] 6, 1 \[DirectedEdge] 7, 1 \[DirectedEdge] 8, 1 \[DirectedEdge] 9, 1 \[DirectedEdge] 10, 2 \[DirectedEdge] 4, 2 \[DirectedEdge] 6, 2 \[DirectedEdge] 8, 2 \[DirectedEdge] 10, 3 \[DirectedEdge] 6, 3 \[DirectedEdge] 9, 4 \[DirectedEdge] 8, 5 \[DirectedEdge] 10}; GraphPlot[as, PlotTheme -> "ClassicDiagram", DataRange -> {{0, 1}, {0, 1}}, GraphLayout -> "CircularE

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

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

请教大佬们，如何从4-7和4-8联立得到式子4-9

下面是代码: n = 4; Fig = Table[1, {i, 1, n}]; For[l = 1, l <= n, l++, p1 = 1; p2 = p1 (1 + 10^(-l - 13)); Clear[{s1, s2}]; s1 = y /. NDSolve[{y'[x] == y[x] Cos[p1 x + y[x]], y[0] == 1}, y, {x, 0, 31}, AccuracyGoal -> 20, PrecisionGoal -> 20, WorkingPrecision -> 64][[1]]; s2 = y /. NDSolve[{y'[x] == y[x] Cos[p2 x + y[x]], y[0] == 1}, y, {x, 0, 31}, AccuracyGoal -> 20, PrecisionGoal -> 20, WorkingPrecision -> 64][[1]]; Fig[[l]] = Plot[{s1[x] - s2[x]}, {x, 0, 30}, PlotRange -> Full]; ] Fig 试了很多方法, 但是, 只要p1 和p2 的区别小于10^(-17)左

吧里关心微分方程的人不多，但这个例子确实很独特，姑且提一句。我们知道（好吧，相当多的人不知道，想着自己写经典龙格库塔的人时不时就会冒出来），NDSolve默认的常微分方程求解器（细节不明——甚至可以说是NDSolve中最不透明的部分——似乎一般是LSODA方法，有时似乎也会调用StiffnessSwitching方法）性能是非常好的，在NDSolve求解出错的时候，常微分方程求解器永远应该放在最后检查。而当常微分方程求解器出现性能问题时，“使用MaxStepSize限

已知 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，希望各位大佬帮忙指导一下，谢谢！）

后半截（不许问这个是什么或是怎么用）： s/1ytBfx9aZZLQvhSBvt9YqgA?pwd=wisz 两本书分别是 《Fundamental Finite Element Analysis and Applications: with Mathematica and Matlab Computations》 《Advanced Topics in Finite Element Analysis of Structures: With Mathematica and MATLAB Computations》 书的优点是讲解极为详细，几乎看不到跳步，即便撇开编程的部分也是很不错的有限元教材（尤其是第一册）。此外随书附赠的源代码（上面的链接中已包含）也是不错的学习材料。

例如：运行 Result1 = f11[t] /. t -> 33之后 得到结果 47.8 [33] 请问怎么把中括号及“33”去掉？ 十分感谢！

怎么修改自动生成的横坐标轴刻度之间的间隔，想把间隔缩小一点，就是0到0.2之间的距离近一点，请问如何设置 代码： Plot[{f13, f14, f15}, {\[Lambda], 0, 1}, PlotLegends -> {"s=0.3", "s=0.5", "s=0.7"}, Frame -> True, FrameLabel -> {"\[Lambda]", "\!\(\*SubsuperscriptBox[\(\[CapitalDelta]p\), \(t\), \(DT\)]\)", "b=0.2", None}, PlotStyle -> {{Black, Dotted}, {Black, DotDashed}, {Black, Dashed}}]

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

版本为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*)

如第二张图里下划黄线的部分，为什么会出现这样的不规则表示。 同时也说一下我在其他的mma文件里写一些什么好像都不会出现这，只有这一个会这样。 看看第一张图 d\[Theta]dt 这条输入命令和 k1 这条命令输入和第二张的那几处对应，有没有什么问题呢？我是没有发现。 L = 1.6 H = 4.3 c = 0.7 a = 0.7 b = 3.2 l = 2.5 m1 = 0.9 d = 1.3 p = 0.6 m = 4.2 \[Alpha]0 = 0.6 \[Beta]0 = 0.6 \[Alpha]10 = 0.2 \[Beta]10 = 0.1 g = 9.79 \[Theta] = ArcSin[-(y[t] + l*(Cos[\[Alpha][t]] - Cos[\[Beta][t]]))/L] Sin[\[Theta]] d\[

我有一个非对易乘法的表达式， \[Sigma][y, 3] ** \[Sigma][y, 2] ** \[Sigma][y, 1] 我想利用循环替换指令ReplaceRepeated将上式变为 \[Sigma][y, 1] ** \[Sigma][y, 2] ** \[Sigma][y, 3]， 即每一个Sigma[y,i]按照i的大小，从小到达排列相乘。 以下是我的MMA代码： *************代码区********************* \[Sigma][y, 3] ** \[Sigma][y, 2] ** \[Sigma][y, 1] //. \[Sigma][y, i_] ** \[Sigma][y, j_] :> If[i > j, \[Sigma][y, j] ** \[Sigma][y, i], \[Sigma][y, i] ** \[Sigma][y, j]] ************代码区********************* 然而运行结果是//.指令

把求根公式换成chareqn2/chareqn3 跑的时候会显示计算溢出 搞不明白 有没有老哥救一下 底下放源码 Clear; \[Lambda] = 0.633; nSiO2 = Sqrt[ 1 + x1*\[Lambda]^2/(\[Lambda]^2 - x2^2) + x3*\[Lambda]^2/(\[Lambda]^2 - x4^2) + x5*\[Lambda]^2/(\[Lambda]^2 - x6^2) /. {x1 -> 0.6961663, x2 -> 0.0684043, x3 -> 0.4079426, x4 -> 0.1162414, x5 -> 0.8974794, x6 -> 9.896161}]; nSi = Sqrt[ 11.6858 + 0.939816/\[Lambda]^2 + 0.000993358/(\[Lambda]^2 - 1.22567)]; n1 = nSiO2; n2 = 1; k0 = 2*\[Pi]/\[Lambda] c = 2.99979*10^14; d0 = 0.12; (*0.12um=120nm*) \[CapitalDelta]d =

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

如图所示，我使用的mathematica版本是10.1 已经使用了两个多月，一直没有发现有什么问题，最近使用的时候发现很简单很简单的计算，发出运行命令后，软件一直在running，无法完成计算。计算机也没有发现有发热之类别的故障，其他程序运行十分正常。 请问广大吧友有没有遇到过这种情况，有没有什么办法可以解决？ 非常感谢。

2. \!\(\*SuperscriptBox[\(A\), TagBox[ RowBox[{"(", RowBox[{"0", ",", "1"}], ")"}], Derivative], MultilineFunction->None]\)[\[Alpha], \[Beta]] \!\(\*SuperscriptBox[\(B\), TagBox[ RowBox[{"(", RowBox[{"1", ",", "0"}], ")"}], Derivative], MultilineFunction->None]\)[\[Alpha], \[Beta]] \!\(\*SuperscriptBox[\(\[ScriptW]\), TagBox[ RowBox[{"(", RowBox[{"1", ",", "1"}], ")"}], Derivative], MultilineFunction->None]\)[\[Alpha], \[Beta]] == 2. A[\

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

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

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

先贴代码 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才能发表评论”的主题。

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

代码：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,

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