题主正在用mma解决物理问题，求解带有角速度的球体的运动轨迹。将方程组输入后，方程“y'[0]==36.53”被替换为True,删除后提示方程约束条件的数目不等于系统的总微分阶数 （6），请大家帮忙看看如何解决这个问题。 R = 0.105; M = 0.45; Cd = 0.24; Cw = 1.7; g = 9.8; S = 0.0346; \[Rho] \ = 1.3; w0 = 88 ans = NDSolve[{M* x''[t] == -(1/2) \[Rho] Sqrt[x'[t]^2 + y'[t]^2 + z'[t]^2] x'[t] S Cd - \[Rho]*R^3 w0*y'[t] Cw, M*y''[t] == -(1/2) \[Rho] Sqrt[x'[t]^2 + y'[t]^2 + z'[t]^2] y

鉴于这个问题在吧里的出现率也有明显的上升趋势，特开这个标题比较好懂的主题，以便今后查找及索引。简单的答案： Clear@Derivative 就行了。 想对有更深入的了解请阅读这个：http://mathematica.stackexchange.com/a/46239/1871

这是原方程， A = Sqrt[(1/(t^2) - w^2) (dVal^2 + dh2) + 1]/ Sqrt[(1/(t^2)) (dVal^2 + dh2) + 1]; z[dVal_] := t*Integrate[1, {A, 0, rh}，PrincipalValue -> true] 这是写的相关程序。学业需要用该方程作图，但是发现积不出来，是代码问题吗？应该如何处理 z[d_] := t \!\( \*SubsuperscriptBox[\(\[Integral]\), \(0\), \(rh\)]\(1 \[DifferentialD] \*FractionBox[\(d* \*SqrtBox[\(\((1/ \*SuperscriptBox[\(t\), \(2\)] - \*SuperscriptBox[\(w\), \(2\)])\) \(( \*SuperscriptBox[\(d\), \(2\)] + dh2)\) + 1\)]\), SqrtBox[\(\(1/ \*SuperscriptBox[\(t\), \(2\)]\)/\(( \*Supe

Solve[1/2 + AB^2 - 1/(8 A^2) - (9 B^2)/(20 A^2) - 9/(80 A^3) == 0, B],为什么显示B=0.595375，然后显示FAlSE，B=0.595375是我上组计算出来的数据，我新建笔记本后还是不行，重新输入公式也不行，这是为什么呢？

Clear["Global`*"]; delta[a_] := Z^2*(1 - 2 a)/(a*(Z + 1)); a11[a_] := (3*Pi/32)*(288 + 604*Sqrt[2]*delta[a] + 217*delta[a]^2)/ 4/(72 + 61*Sqrt[2]*delta[a] + 16*delta[a]^2); dc[a_] := Piecewise[{{1, a <= 0}, {a11[a]/(1 - a*(2 - mi/mI*(Z + 1)))/(a*(Z - 1) + 1), a > 0 && a <= 1/2}, {dc[1/2], a > 1/2}}]; dp[a_] := dc[a]*a*(1 - 2 a)*(Z - 1)^2/(Z - a*(Z - 1)); D0 = 278*4/100*Sqrt[2]*5^(5/2)/(2*1.81)*(Z + 1)/Z^2; (*cm^2/s*) Z = 6; mi = 2; mI = 36; a0[z_] := Piecewise[{{1/2, z < 0}, {0, z >= 0}}]; sol = NDSolve[{D[a[z, t], t] == D[D0*(dc[a[z, t]] + dp[a[z, t]])*D[a[

怎么把表达式里的\[Phi]提出来，将这个表达式表达成\[Phi]和其他部分的乘积 n2 = 1/(2*\[Omega]22) (-(1/\[Tau]2 + ( c \[Sigma] f1 \[Phi]*Abs[U[x, y, z]]^2)/nn) + Sqrt[(1/\[Tau]2 + (c \[Sigma] f2 \[Phi]*Abs[U[x, y, z]]^2)/ nn)^2 + 4*\[Omega]22*(( c \[Sigma] f1 \[Phi]*Abs[U[x, y, z]]^2)/nn*n1 - 1/\[Tau]1*n1 + 2*R[x, y, z])]);

楼主是学物理的，尝试手解在光滑桌面上的小球在一端固定且有原长l的弹簧下做的二维运动，结果呵呵了

没有求出来解，也没有显示报错

这里的内置TeX功能报错了，不知道是什么原因。版本是13.2 windows。我同学的mac是可以使用的，我这个不管输什么都红，如果换成13.1可以解决吗？

Manipulate[ DensityPlot[ 1 + Cos[(4 \[Pi]/(\[Lambda]*10^-9)*((d + (x - r)* Sin[(\[Pi]*\[Theta])/180])*10^-3* f/Sqrt[x^2 + y^2 + f^2]))], {x, -3, 3}, {y, -3, 3}, PlotPoints -> 200, PlotRange -> All, ImageSize -> 1.05 {400, 400}, PerformanceGoal -> "Speed"], Row[{Control@{{\[Lambda], 632.8, "波长\[Lambda](nm)"}, 450, 750, Appearance -> "Labeled", ImageSize -> Small}, Control@{{d, -0.15, "空气膜厚度d(mm)"}, -0.2, 0.2, Appearance -> "Labeled", ImageSize -> Small}，Control@{{f, 20, "透镜焦距"}, 10, 40, A

Manipulate[ DensityPlot[ 1 + Cos[(4 \[Pi]/(\[Lambda]*10^-9)*((d)*10^-3*(Sqrt[ m^2 - n^2*Sqrt[(x^2 + y^2)/(f^2 + x^2 + y^2)]])))], {x, -3, 3}, {y, -3, 3}, PlotPoints -> 200, PlotRange -> All, ImageSize -> 1.05 {400, 400}, PerformanceGoal -> "Speed"], Row[{Control@{{\[Lambda], 632.8, "波长\[Lambda](nm)"}, 450, 750, Appearance -> "Labeled", ImageSize -> Small}, Control@{{d, 0.027, "空气膜厚度d(mm)"}, -0.2, 0.2, Appearance -> "Labeled", ImageSize -> Small}, Control@{{f, 20, "透镜焦距"}, 10, 40, Appearanc

u = 0.33; Er = 1/(12.3*10^-12*(1 - u^2)); v = Sqrt[Er/7500]; a = 0.015; l = 0.005; e = 15.1; w = k*v; x = k*a; Y = (I*w*e*Pi*a^2)/ l*(1 - 0.58^2 + 0.58^2*((1 + u)*BesselJ[1, x])/( x*BesselJ[0, x] - (1 - u)*BesselJ[1, x])); Z = 1/Y; f = x/(2*Pi*a)*Sqrt[1/((1 - u^2)*7500*12.3*10^-12)] NSolve[{x*BesselJ[0, x] == (1 - u)*BesselJ[1, x]}, k] Plot[Norm@(Z), {k, 0, 700}] 555.1031891081523`*137.82424095932583` 555.1031891081523`*359.6737230889798` 555.1031891081523`*571.693401597557`

In[8]:= Abs[(4*k1*k2*Exp[I*k1*d])/( 4*k1*k2*Cos[k2*d] - 2*I*(k1^2 + k2^2)*Sin[k2*d])] Out[8]= 4 E^-Im[d k1] Abs[(k1 k2)/( 4 k1 k2 Cos[d k2] - 2 I (k1^2 + k2^2) Sin[d k2])] In[9]:= ComplexExpand[ 4 E^-Im[d k1] Abs[(k1 k2)/( 4 k1 k2 Cos[d k2] - 2 I (k1^2 + k2^2) Sin[d k2])]] Out[9]= (4 Sqrt[k1^2] Sqrt[k2^2])/Sqrt[ 16 k1^2 k2^2 Cos[d k2]^2 + 4 (k1^2 + k2^2)^2 Sin[d k2]^2] 楼主刚学，看了半天帮助文档也没明白，上面那个指数函数为啥求模量就变成虚部然后没了。

给出五个积分。 Integrate[p*BesselK[2, p] BesselK[2, p], {p, 0, Infinity}] Integrate[p*BesselK[2, p] BesselK[2, p], {p, 0, Infinity}, Generat

为什么有时候“，”和“&&”不能互换，有时候又可以互换？

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

请问各位大神，图中的方程组要怎么解？ NSolve[a == 15/(1 - Exp[-74 b]) && 15 == 74 a b Exp[-74 b] + 550 b, {a, b}]

这应该算是使用 ReplaceAll （ /. ）的过程中必踩的两个坑之一。（另一个请参看《为什么 a^4 /. a^2 -> b 不输出 b^2 ？说说语义式匹配》：https://tieba.baidu.com/p/5728734141 。）为方便修订，正文内容放在了云笔记上，请多点一下鼠标： http://note.youdao.com/s/8cleVUiT

自己输入的代码是这样的 In[6]:= \[Psi][\[Xi]] = E^(-\[Xi]^2/2) H[\[Xi]]; H[\[Xi], Subscript[a, \[Nu]]] = Sum[Subscript[a, \[Nu]] \[Xi]^\[Nu], {\[Nu], 0, \[Infinity]}]; \[Psi][\[Xi]] Out[8]= E^(-(\[Xi]^2/2)) H[\[Xi]] In[9]:= \[Psi]'[\[Xi]] = D[E^(-(\[Xi]^2/2)) \!\( \*UnderoverscriptBox[\(\[Sum]\), \(\[Nu] = 0\), \(\[Infinity]\)]\( \*SuperscriptBox[\(\[Xi]\), \(\[Nu]\)]\ \*SubscriptBox[\(a\), \(\[Nu]\)]\)\), \[Xi]] Out[9]= E^(-(\[Xi]^2/2)) \!\( \*UnderoverscriptBox[\(\[Sum]\), \(\[Nu] = 0\), \(\[Infinity]\)]\(\[Nu]\ \*SuperscriptBox[\(\[Xi]\), \(\(-1\) + \[Nu]\)]\ \*SubscriptBox[\

在使用FindFit函数时，我的模型是一个含有多个参数的12阶矩阵。我希望用{kx,ky,kz}为自变量，矩阵的第6个本征值为因变量Re[Eigenvalues[ham][[6]]]，拟合最佳参数。但是，在求解方程本征值时候因为计算量太大，解不出来。老师提到，可以对矩阵数值化之后再去做拟合，可以降低计算量，不太理解。请问各位老师有没有好的方法能求解本征值？ 实际上这是一个紧束缚模型，考虑一个原子周围的原子对自己的影响。考虑的原子越多，参数越多，但是矩阵的大

我想要使用一个局部放大的图像。以下是代码。正如您所看到的，在放大图像中，“Min”和“RF”等“PlotLabels”未能完全显示。然而，在可以插入“PlotLabels”的位置存在明显的空白空间。如何解决这个问题呢？ ``` ttt = Table[{5, 10, 20, 50, 100, 200, 500, 1000, 2000, 5000, 10000, 50000, 100000, 500000, 1000000}] ttttt = {0.5411080105624226`, 0.4746114896441091`, 0.45560476211921547`, 0.3920375515864554`, 0.3688071154680211`, 0.37303695707239365`, 0.35986703180178964`, 0.3454025968800972`, 0.33737968696907666`, 0.3328931553201

看了帮助文档和网上的经验贴都没解决问题 无奈之下向贴吧各位大佬们求助 谢谢各位 报错截图如下 具体代码文本如下 v0 = 100; Rp = 0.1; R = 0.025; \[Nu] = 0.3; \[Mu] = 8.5; EE = 50*10^9; mp = 0.03; h = 0.005; G = 2*10^10; \[Rho]s = 1200; k1 = 1; k2 = 1; \[Sigma]ys = 50*10^6; n1 = 1; n2 = 1; \[Omega][r_, t_] := \[Omega]0[ t]*((2*r^3)/\[Eta][t]^3 - (3*r^2)/\[Eta][t]^2 + 1); \[Gamma]rz[r_, t_] := \[Gamma]rz0[t]*Sin[(Pi*r)/\[Eta][t]]^2; \[Rho][Z] = 0.8*Z + 0.1; \[Rho]Z0 = Integrate[\[Rho][Z], {Z, 0, 1}]; \[Sigma]0[\[Rho]] = k1*\[Sigma]ys*\[Rho]Z0^n1; CC[\[

搜了一下，发现吧里好像没有主题帖专门说这个。（尽管我觉得这应该是任何一个有心学习Mathematica的人都能很容易地接触到的常识。）那就专开一帖好了。 为了便于修订，照例放到云笔记上： note.youdao.com/noteshare?id=7334ad886b97de3565d59beea50b5350 ……其实不看云笔记都无所谓，概括起来就那么几句话： 1. For循环写起来很冗长——这其实是其次的； 2. For循环在Mathematica里很慢——这其实也是其次的； 3. For循环的表意逻辑，几乎在任何情况下，都无意义地复

E0 = 1; s0 = 0.5; lambda = 0.0006328; r = 0; zr = Pi*s0^2/lambda; z = zr*1; sz = s0*Sqrt[lbk]1 + ((z*lambda)/(Pi*s0^2))^2[rbk]; k0 = 2*Pi/lambda; Rz = z (1 + ((Pi*s0^2)/(z*lambda))^2); phiz = ArcTan[lbk](lambda*z)/(Pi*s0^2)[rbk] E1 = (E0/sz)*Exp[lbk]-1*(r^2/sz^2)[rbk]* Exp[lbk]I*k0*(z + (r^2/(2*Rz))) - I*phiz[rbk]; 我想计算r从0到2，取值间隔0.01，依次计算200个值，要怎么样实现呢

DSolve[lbk]{\[lbk]InvisiblePrefixScriptBase[rbk] Subscript[lbk]a, 0[rbk] u'[lbk]t[rbk] + \[lbk]Lambda[rbk] Subscript[lbk]a, 0[rbk]^2 Subscript[lbk]a, 1[rbk] \[lbk]Alpha[rbk][lbk]t[rbk] + \[lbk]Lambda[rbk] Subscript[lbk]b, 0[rbk]^2 Subscript[lbk]b, 1[rbk] \[lbk]Beta[rbk][lbk]t[rbk] + \[lbk]Lambda[rbk]^3 Subscript[lbk]a, 1[rbk] Subscript[lbk]q, 2[rbk] \[lbk]Gamma[rbk][lbk]t[rbk] == 0, 2 \[lbk]Lambda[rbk] (Subscript[lbk]a, 0[rbk] Subscript[lbk]a, 1[rbk]^2 \[lbk]Alpha[rbk][lbk]t[rbk] + Subscript[lbk]b, 1[rbk]^2 Subscript[lbk]b, 0[rbk] \[lbk]Beta[rbk][lbk]t[rbk]) == 0, \[lbk]Lambda[rbk] (Subscript[

1.电脑硬盘空间不多，就想着清理掉本地文档。 2.已经直接删除了文档所在文件夹，但是不论F1还是主页文档，都没有响应。 右键获取帮助可以打开，但我还是希望能够把F1恢复正常。 3.网速不好的朋友请勿模仿。

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

MMA中的PDE的NDslove是用的有限元来求解的吗？

mathematic积分时，对被积函数原封不动的输出，并没有做任何变化，求各位大佬帮忙看看这是什么问题，已经困扰好久了，下面是代码 \[lbk]Lambda[rbk] = 1;r = 1/2;\[lbk]Theta[rbk] = \[lbk]Gamma[rbk]/2*(t + (E^(-\[lbk]Lambda[rbk]*t) - 1)/\[lbk]Lambda[rbk]);\[lbk]Alpha[rbk] = \[lbk]Gamma[rbk]/2*(1 - E^(-\[lbk]Lambda[rbk]*t));d\[lbk]Rho[rbk]11 = -E^(-2*\[lbk]Theta[rbk]*Cosh[lbk]2*r[rbk])*\[lbk]Alpha[rbk];d\[lbk]Rho[rbk]10 = \[lbk]Alpha[rbk]/ 4 (E^(-2*r - E^(2*r)*\[lbk]Theta[rbk]) (E^(I*\[lbk]Beta[rbk]) - 1) - E^(2*r - E^(2*r)*\[lbk]Theta[rbk]) (E^(I*\[lb

这个也是讲烂了的内容，但吧里好像没有标题足够显眼的帖，姑且专开一帖，方便索引。 比较常用的“超级函数”基本都有与之对应的数值求解器。如标题所说，如果符号解没算出来并且符号解对你来说不是必要的，那么不妨去试试数值解。 只写这么几行的话内容好像太少，那就再稍微谈一下比较常见的成对的符号/数值求解器及其特性。 对于代数方程：Solve（虽然现在已经发展成了超级求解器，但它主要是个多项式方程求解器）/NSolve（引入时间非

求助各位吧友们，请问List1中的第一个元素为什么不能被识别出来呢？感谢~ List1 = {Subscript[\[CurlyPhi], 1], Subscript[\[CurlyPhi], 2], \[CurlyPhi], Subscript[\[CurlyPhi], 1]*\[CurlyPhi], \[Theta], \[CurlyPhi]1} Table[Cases[List1[[i]], Subscript[\[CurlyPhi], 1], Infinity], {i, Length[List1]}] TreeForm[List1]

我的问题需要用到mellin逆变换，但没有解析解，所以想用个简单算例比较下数值积分与解析解的误差。代码如下： r = 2; f[s_] := Gamma[s]*Gamma[1 - s] InverseMellinTransform[f[s], s, r] NIntegrate[1/(2 \[Pi] I) f[s] r^-s, {s, 1 - 5 I, 1 + 5 I} (*Exclusions->{0}*), WorkingPrecision -> 20, MaxRecursion -> 20] // Timing NIntegrate[ Re[1/(2 \[Pi]) f[1 + p I]/ r (Cos[p Log[r]] - I Sin[p Log[r]])], {p, -5, 5} (*Exclusions->{0}*), WorkingPrecision -> 20, AccuracyGoal -> 8, MaxRecursion -> 100, PrecisionGoal -> 100, Method -> {"Glob

y1[lbk]t_[rbk]y2[lbk]t_[rbk] m = 1, c = 1, l = 1, Jc = 1, r = 1 xs[lbk]t[rbk] = -((l - r - c)*sin[lbk]y1[lbk]t[rbk][rbk] + r*y1[lbk]t[rbk]) ys[lbk]t[rbk] = (l - r - c)*cos[lbk]y1[lbk]t[rbk][rbk] + r D[lbk]xs[lbk]t[rbk], t[rbk] 报错显示第二行给的不完整

现在我对这个微分方程的解猜出了一个形式解，如图所示，现在我想要级数展开然后求解前面的常系数，但是使用Series函数只能得到幂级数的展开解，带log z的这些项体现不出来，我要怎么做才能得到我想要的形式展开呢？ 以下是代码： In[4]:= SetAttributes[{u1, u0, zh, v4}, Constant] In[2]:= eq1 = y''[z] + (-3/z + 2*u1^2*z + 4*z^3/(z^4 - zh^4)) y'[z] + 1/z^2*zh^4/(zh^4 - z^4) (3*y[z] - 4*v4*y[z]^3) == 0 Out[2]= (zh^4 (3 y[z] - 4 v4 y[z]^3))/( z^2 (-z^4 + zh^4)) + (-(3/z) + 2 u1^2 z + (4 z^3)/( z^4 - zh^4)

具体来说，用多项分布生成了一个两层的数据集，现在要对生成元素计数。用tally或counts都可以计数，不过我的需求是像{0,9,0}与{9,0,0}、{0,0,9}这样的项认为是等价的，直接给出这些计数之和。我想如果改变List的属性为Oderless就可以实现这个目的，不知道有没有其他的技巧？ data = RandomVariate[MultinomialDistribution[9, {1/3, 1/3, 1/3}], 10^5]; Counts[data]; Tally[data];

二分法查找一直不能正确输出啊，是用Which语句实现递归的过程错了吗 下面是代码：： temp = 10^-6;(*设置一个阈值来判断函数值是否接近零*) Findzero[f_, start_, final_, temp_] := Module[{a = N[start], b = N[final], c}, c = (a + b)/2; fc = f[c]; Which[Abs[a - b] <= temp, Return[c], fc <= temp, Return[c], f[a]*fc > 0, Findzero[f, c, b, temp], f[b]*fc > 0, Findzero[f, a, c, temp]]]; f[x_] := Sin [x]; Findzero[f, -1.5, 0.5, temp] 希望大佬们能够为我点拨一二。

Clear["Global`*"] V[x_, y_] := {Integrate[\[Lambda]/Sqrt[(x - l)^2 + y^2], l] /. {l -> L}} - {Integrate[\[Lambda]/Sqrt[(x - l)^2 + y^2], l] /. {l -> -L}} Ex[x_, y_] = -Div[V[x, y], x] // Simplify Ey[x_, y_] = -Div[V[x, y], y] // Simplify L := 1; \[Lambda] := 1; P=ContourPlot[V[x, y], {x, -5, 5}, {y, -5, 5}, Contours -> 10,AxesLabel->Automatic] G=VectorPlot[{Ex[x, y],Ey[x, y]},{x, -5, 5}, {y, -5, 5}] Show[P,G]

大佬们，这个我算不出来，能帮我看看吗？eqn = D[lbk]u[lbk]t, x, y[rbk], t[rbk] == D[lbk]u[lbk]t, x, y[rbk], x, x[rbk] + D[lbk]u[lbk]t, x, y[rbk], y, y[rbk] - u[lbk]t, x, y[rbk] + u[lbk]t, x, y[rbk]^3;ic = {u[lbk]0, x, y[rbk] == Exp[lbk]-(x^2 + y^2)[rbk], Derivative[lbk]1, 0, 0[rbk][lbk]u[rbk][lbk]0, x, y[rbk] == 0};sol = NDSolve[lbk]{eqn, ic}, u, {t, 0, 1}, {x, -5, 5}, {y, -5, 5}[rbk];Plot3D[lbk]Evaluate[lbk]u[lbk]t, x, y[rbk] /. sol[rbk], {t, 0, 1}, {x, -5, 5}, {y, -5, 5}, PlotRange -> All[rbk]

um=4; list=Table[ParametricPlot3D[{Cos[u]+n/2,Sin[u]+n/2,v}, {v,0,1},{u,0,um},Axes->False,Boxed->False, PlotStyle->Directive[FaceForm[Green,Texture["ABCDEFG"]]], TextureCoordinateFunction->({2#5,#4}&),Lighting->"Neutral", Mesh->None,PlotRange->All,ViewPoint->{-30,10,20}, ViewAngle->.025],{n,0,2}] Show[list[[1 ;; 3]]]

各位好！我用网上一个生成mma3 spikey的代码玩mma，发现生成图形的缺省图，把理论上看不见的线也给展现了出来。如下图所示，杂线很多，尤其是最外边。 转角度后，上述问题有所减轻，但还是有点线头。 想问下这是由什么原因造成的？有什么解决办法吗？谢谢大家！ 附：我用的mma是11.3破解版，代码（从网上扒的）： Graphics3D[{EdgeForm[{GrayLevel[0.25], Thickness[0.001]}], SurfaceColor[GrayLevel[0.8], RGBColor[1, 0.4, 1], 3], {Polygon /@ #1, Polygon /@ #2, MapThread[ Polygon[Join[#1, #2]]

请问这种公式是怎么在mma中输入的？

phi1关于theta1的表达式

本帖是果壳网友 马甲与小号 发在果壳小组的一篇入门教程。如今果壳小组已关，故在吧里发一份备份。 note.youdao.com/noteshare?id=c5114ccda3270199fd801952cf785bff

例如u为x和y的函数，s等于u对x的偏导加z乘以u对y的偏导，麻烦告知一下谢谢

在后缀为.nb的Mathematica文档中,插入内嵌TeX输入,出现十分严重的乱码,格式不正确,字体不协调,汉字无法显示的问题. 对于代码： \begin{array}{cl}当r < R时, &\vec {E} = 0 \\ 当r > R时, &\vec {E} = \frac {Q} {4\pi\varepsilon_ 0 \r^2}\vec{e_r}\end{array} 对比图如下: 此外，\[CurlyEpsilon],也就是\Varepsilon符号也有非常明显的问题,上两图也有显示,在下面附带一个更明显的: 而且,\mathrm{}这种很常见的把字体变成罗马体的指令,也会给人一种排版上的非常难受的感觉,仿佛额外加粗

问题总结：给定几个点，绘制从原点到这几个点的垂直平分线，然后几个向量的垂直平分线相交围成的最小面积（三维情况下就是体积）的图形就是需要绘制的图形。如何能绘制出这个图形 背景：打算绘制一个布里渊区示意图；布里渊区图绘制方法：给定两个基矢量b1,b2；绘制所有G=n1*b1+n2*b2，其中n1,n2为整数，比如我们给定n从-10取值到10，一个晶格常数为1的简单格子的布里渊区就是{{-10pi,-10pi},{-10pi,-9pi},...,{-10pi,10pi},{-9pi,-10pi},{-9pi,-9pi},...,{-9pi,10pi},......

方程不是原本的方程，我为了debug删得只剩这么点了，居然还报错，删到只剩一边了也报错，到底为什么啊？我看这个报错的官方解释就是自变量的顺序要一致，比如f[x, t]不能写成f[t, x]，但我这没有这个错误啊 DSolve[{D[c[t, \[Zeta][t]], t] == D[c[t, \[Zeta][t]], \[Zeta][t]]}, c, {t, \[Zeta][t]}] DSolve[{D[c[t, \[Zeta][t]], t] == 1}, c, {t, \[Zeta][t]}] DSolve[{1 == D[c[t, \[Zeta][t]], \[Zeta][t]]}, c, {t, \[Zeta][t]}]