请问这是怎么回事？如何理解给出的这个解？ 代码：RSolve[lbk]{a[lbk]n[rbk] == COS (a[lbk]n - 1[rbk]) + 1, a[lbk]0[rbk] == e}, a, n[rbk] Solves a recurrence equation for a[lbk]n[rbk] .

x = 3 low = 1 high = 2 guess[low_, high_] := (low + high)/2 e = 0.01 f[x_] := If[guess[low, high]^2 > x, high = guess[low, high], low = guess[low, high]] While[ Abs[high - low] > e, f[x]] guess[low, high] // N

现在是手动二分法，输入两个自变量x,y，寻找输出的两个量为特定的值时与之对应的xy，现在想能不能在mma里或者python中二分法求解x,y。已知函数单调，和xy大致范围，不会出现不同xy对应同一个值

c = 0.5 \[Alpha] = 0.9 d = 4 k = 4 \[Beta] = ( d (4 - \[Alpha]) (\[Alpha] - 2) + 4 [3 w + c (\[Alpha] - k)] - \[Alpha] (c k + w - \[Alpha]))/( c k (4 - 3 \[Alpha])) Plot[(2 w + \[Alpha] (1 + c) + 2 k c (\[Beta] - 1))/(\[Alpha] (4 - \[Alpha])), {w, 3.7, 4.8}]

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

各位大佬，这是怎么回事？ 代码：RSolve[lbk]a[lbk]n[rbk] = a[lbk]n - 1[rbk] + 1, a[lbk]n[rbk], n[rbk] Solves a recurrence equation for a[lbk]n[rbk] .

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

下图是运算结果，我希望更改频率和线径来得出一个量来进行比较，假如如下图输入所示fs为3＊1的矩阵，线径为1＊3的矩阵，最后不应该α是一个3*3的矩阵，为啥是现在这种情况，请各位大佬帮忙看看 代码如下 fs = {70 Quantity[lbk]1, "Kilohertz"[rbk], 140 Quantity[lbk]1, "Kilohertz"[rbk], 150 Quantity[lbk]1, "Kilohertz"[rbk]};(*开关频率*) WireDs = {{0.1 Quantity[lbk]1, "Millimeters"[rbk], 0.07 Quantity[lbk]1, "Millimeters"[rbk], 0.05 Quantity[lbk]1, "Millimeters"[

重积分我想用那个boole函数表示结果很多简单的都积不出来，第二形曲面积分完全不知道该怎么用mma来算

你好，我在学习mma的时候和朋友一起算一些积分公式。 发现在他的电脑上可以计算出结果，但是在我的电脑上只能得出还是一个定积分的式子。 其他步骤代码也完全相同，之前计算出来的结果也相同。就是Z2的不同，甚至相同的计算步骤出来的W2结果也是一样的。朋友说他也遇到过类似的情况，是关掉再重新运行就可以，但是我尝试了清除重新计算、关掉、换成英文版、重启电脑等都不行。 Derivative[1][z2][t] == 1/6 (Sqrt[3] w1[t]^2 - 2 Sqrt[3] w1[t] z1[t] - 2 Sqrt[6]

类似这种，怎么分别按照上下两个部分分别进行拟合？

代码如下：FullSimplify[ ImplicitD[x - y, (x - y)/(2 z) == n[x]/p[k[x], z] - n[y]/q[k[y], z], x - y, z]] 想要求d(x-y)/dz，其中p和q是k和z的函数，n和k是x或y的函数，求大佬帮忙看看是不是有什么问题，谢谢

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

Limit[lbk]Sum[lbk]n^k,{k,1,n}[rbk]/Sum[lbk]k^n,{k,1,n}[rbk],n->Infinity[rbk] 这样得不到结果，应该怎样输入才能求解呢

现在有如图所示的数据点，z方向每隔0.05有一组数据分立分布的，每一组都形成了一条二维封闭曲线（原先数据量z方向上范围比较大，数据点已经够密了，为了减少贴出来的数据量截取了一小部分） 通过对数据的辐角和z坐标值排序，可以用ListLinePlot3D画出连成线的版本，效果还不错（部分地方变化较快有折线，可能也是后面画曲面有问题的原因） 想请问一下怎么让这些点（或二维线）生成一个三维面，类似一个筒的效果？ 用ListSurfacePlot3D画出来是很

求助有关ListVectorPlot的一些问题，下面的123分别对应着三幅图产生的问题： 1.首先Rainbow最亮的颜色是红色，下面这串代码将大于0.5的部分全部设置为最亮的颜色红色，但其实对于EVPA中的最大值是3.6，这样导致形成的图像大面积的全是红色（存在的问题可能是数据中存在较大的值，而大部分值相对较小，那么图例的颜色范围和实际图像中的颜色会不一致，个人猜测，不一定。）我想要的结果是只在最大值附近才是红色，让红色面积减少。 2.我通过将Vecto

如何使得f[g[x_]]的定义始终优先于g[x] 如上，我期望f[g[2]]始终输出4 代码如下 ClearAll["Global`*"] f[g[x_]]^:=2x; f[g[2]] g[x_]:=x f[g[2]]

，求结微分方程出现这种结果，怎么处理，这是代码，{{\[lbk]Alpha[rbk][lbk]t[rbk] -> InverseFunction[lbk]-(( g^2 \[lbk]CapitalDelta[rbk] ArcTan[lbk](2 Sqrt[lbk]k[rbk] Sqrt[lbk]#1[rbk])/ Sqrt[lbk]-2 \[lbk]Delta[rbk] + g^2 \[lbk]CapitalDelta[rbk][rbk][rbk])/( 2 Sqrt[lbk]k[rbk] Sqrt[lbk]-2 \[lbk]Delta[rbk] + g^2 \[lbk]CapitalDelta[rbk][rbk])) + Sqrt[lbk]#1[rbk] &[rbk][lbk]-t \[lbk]CapitalOmega[rbk][rbk]}}

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

分别画两幅图，一幅图是实线另一幅是虚线，为什么虚线图例线段长度比实线长，之前我都是直接调整图例框让虚线看起来和实线线段长度一样

已经解了两天了没解出来，有没有大佬帮忙看看怎么回事

我想将nb文件另存为pdf，但是这个pdf的字体怎么选呢？以前都自动是黑体，但是有一次突然变成别的体了一直没变回来，我想要黑体，应该怎么设置呢

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]就好了

代码和运行结果如下： 在吧内尝试着找了一下没有找到相关的帖子，我自己觉得应该是与替换规则的使用或者代码的运算顺序有关系，但是不知道应该去找哪些资料来看，所以如果存在现有的资料的话就拜托大家帮忙指个路了 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

明明是收敛的结果我的MMA非要说它不收敛

\[CapitalOmega] = ImplicitRegion[True, {{z, 0, 200}}]; \[CapitalGamma]2 = NeumannValue[0., x == 0 || x == 200] s = NDSolveValue[{D[n[t, z], t] - D[n[t, z], {z, 2}] - \[Tau]/ nb \[Alpha] (1 - R) Pdensity/(Sqrt[2 \[Pi] ] \[Tau]L hv ) gauss[t, 0, \[Tau]L/\[Tau]] Exp[-\[Alpha] Ldi z] == -Ldi e^2 nb Ldi/(k \ T \[Epsilon] \[Epsilon]b) (0 + eF0O[z]) D[n[t, z] , z] - Ldi e^2 nb Ldi/(k T \[Epsilon] \[Epsilon]b) (p0[z] - n0[z] 0 + 1 - E^(-2 u)) n[t, z] + \[CapitalGamma]2, n[-10 \[Tau]L/\[Tau], z] == 0 }, n, {t, -10 \[Tau]L/\[Tau], 5}, {z} \[Element] \[CapitalOmega], PrecisionGoal -> 35, Method ->

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

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

这个是代码 Plot[0.0172[1 + 3.93 (10^-3)[x - 20]]/10^6, {x, 20, 320}] 我尝试加入PlotRange给因变量范围在1.72*10^-8到4*10^8也不行，查了个帖子说用plot[re]指令，也不行

类似上面所示的雷达图。我手中目前已经有了R、G、B三种颜色的视角数据，请问如何才能将这些数据绘制成图片那样的雷达图呢？尝试过ListDensityPlot和Raster等命令，似乎都不太行。各位大佬有没有知道的，感谢！！！

mma 同样作为 lisp-like 语言有 Throw/Catch 有 Dialog 甚至有 Goto 却没有 call/cc 的官方实现也是够奇怪的

昨天@cckinger 同学跟我提到了遗传算法,我觉得很有意思,便尝试用MMTC实现了一下.第一次玩,题目简单,高手莫笑.

请问刚接触mathematica应该如何学习应用，看书还是看网课，本人大一想接触数学建模相关的比赛，求大佬指点🙏

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

已经把公式中其他内容带入数值，运行不出结果，而且单独计算里面的积分也运行不出来 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

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

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

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

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+

我用的是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

下面是代码: 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)左