k = 1.38/10^23; r = 8.314; model = c*Exp[-(eV1/(k*t))] + a*(1 - Exp[-(eV1/(k*t))])*Exp[-(eV2/(r*t))]; data = {{298, 74.9}, {278, 39.7}, {258, 26.8}, {238, 11.1}, {218,7.}, {198, 3.}, {178, 1.8}, {158, 1.1}, {138, 0.7}}; nlm = NonlinearModelFit[data, model, {c, a, eV1, eV2}, t]; Show[ListPlot[data], Plot[nlm[t], {t, 130, 300}] 代码和图片如上，要拟合的model中指数上方有常数k，有没有k的拟合结果是一样的，但是拟合的方程发生了变化 TraditionalForm[model = nlm[t]] 有k时，结果(98268.67366331698*(1 - E^(-7.246376811594202*^22/t)))/E^(2145.3785836139396/t

Clear["`*"] b = Exp[I (a - I \[Eta]) (\[Epsilon] - C1 (1/(1 - I C2 \[Epsilon]) + 1/(1 + I C2 \[Epsilon])))]; Integrate[b, {\[Epsilon], -Infinity, \[Epsilon]F}, Assumptions -> {a \[Element] Reals, \[Eta] \[Element] Reals, C1 \[Element] Reals, C2 \[Element] Reals}]

请问各位前辈大佬，我有一个方程组，是由8个方程联立，但是未知数是7个，也就是8个方程构成的方程组求解7个未知量，请问用什么指令最好呢？个人试过DSolve、Nsolve、LeastSquares（最小二乘法）、Reduce以及广义逆矩阵的方法等求解，结果要么是零解，要么解不准确。但是我心里知道这个方程在虚数范围内是有解的。请问如何到底怎么求解呢？谢谢。

请问大佬 怎么用mathematica画出跳棋棋盘？

各位大佬，我求解一个简单的微分方程，在Mathematica7.0版本可以直接运算得到结果，但是在14上面，总是这样显示出问题，请问到底是怎么回事？不知道哪里出问题了。谢谢。 AAA = \!\( \*SubscriptBox[\(\[PartialD]\), \( \*SubscriptBox[\(x\), \(1\)], \*SubscriptBox[\(x\), \(1\)]\)]\((Exp[I*\((k* \*SubscriptBox[\(x\), \(1\)] - \[Omega]*t)\)]*f[ \*SubscriptBox[\(x\), \(2\)]])\)\) BBB = \!\( \*SubscriptBox[\(\[PartialD]\), \( \*SubscriptBox[\(x\), \(2\)], \*SubscriptBox[\(x\), \(2\)]\)]\((Exp[I*\((k* \*SubscriptBox[\(x\), \(1\)] - \[Omega]*t)

wholedata = {{0, 394, 372, 11.02`}, {28, 394, 372, 12.7`}, {56, 394, 372, 14.56`}, {84, 394, 372, 6.17`}, {112, 394, 372, 17.25`}, {168, 394, 372, 22.59`}, {224, 394, 372, 21.63`}, {280, 394, 372, 19.34`}, {336, 394, 372, 16.12`}, {392, 394, 372, 14.11`}, {224, 0, 372, 6.39`}, {224, 49, 372, 9.48`}, {224, 98, 372, 12.46`}, {224, 147, 372, 14.33`}, {224, 196, 372, 17.1`}, {224, 294, 372, 21.94`}, {224, 394, 372, 22.64`}, {224, 489, 372, 21.34`}, {224, 587, 372, 22.07`}, {224, 685, 372, 24.53`}, {224, 394, 0, 15.75`}, {224, 394, 47, 16.76`}, {224, 394, 93, 16.89`}, {224, 394, 140, 16.24`}, {224,

有人遇到过这种问题，程序运行后，在重新打开，所有的代码前面都会添加上"In[1] / Out[1]"之类的字样，而且所有的Cell会合并成一个Cell。如图所示：第一张第二张分别是运行后，和关闭重新打开后的

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

和Maple对照很明显，而且用Simplify也无法化简了

如题。自编规则当然也不难，不过这里我想知道的是，能不能用Simplify，TrigExpand一类的化简型指令来得到这个结果。试了半天没试出来……是正好不支持这个化简呢，还是我没有找到？是我漏了什么重要假设吗？（我试验过的假设只有d属于Reals……）

ArcTan[cos[170 D]/(cos[170 D] + 8/3 sin[40 D] sin[80 D])] 怎么得结果 ？

各位大佬，我的代码是这样的 F[q_NumericQ, q0_NumericQ]:=Hold[NIntegrate[f(q,q0,k),{k,0,inf}]] ReleaseHold[NIntegrate[F[q,q0],{q,0,inf},{q0,0,inf}]] 这样得到的结果是0，但它肯定不为0 用Activate和Inactive的方法无法得到答案，或者说计算极慢 如果不用Hold，单纯 F[q_NumericQ, q0_NumericQ]:=NIntegrate[f(q,q0,k),{k,0,inf}] NIntegrate[F[q,q0],{q,0,inf},{q0,0,inf}] 结果得到非数值，这样应该如何解决？

输入：RSolve[{a[n] == \!\( \*UnderoverscriptBox[\(\[Sum]\), \(i = 0\), \(n\)]\((a[i]* a[n - i])\)\), a[0] == a[1] == 1}, a, n] 会报错：RSolve:提供的函数不是给定函数的差分方程 输出：RSolve[{a[n] == \!\( \*UnderoverscriptBox[\(\[Sum]\), \(i = 0\), \(n\)]\(a[ i]\ a[\(-i\) + n]\)\), False == 1/3 == 1}, a, n] 还将a[1]识读成了1/3.将[方程,a,n]的a换成a[n]也没用。求助应该怎么解

想请教一下，为什么不管是式子大于0还是小于0，输出的结果都是false呢？这是什么原因？ Reduce[(-18 m t - (dA kA - dB kB + \[CapitalDelta]q + 3 t z)^2 + (\[CapitalDelta]d kf + \[CapitalDelta]q + 3 t (z + \[Beta]))^2) > 0 && \[CapitalDelta]q > 0 kA > kB > kf && kf > 0 && kA < 1 && dA > dB > 0 && \[CapitalDelta]d > 0 && t > 0 && 0 < z < 1 && 0 < \[Beta] < 1 - z && 0 < z + \[Beta] < 1 && m > 0 && t > 0 && t > (dA kA

拓扑学中的一些拓扑空间，像莫比乌斯带，克莱因瓶这些，它们的图怎么用mathematica比较方便地画？或者说有什么软件能比较方便地画出这些图？

mathematica input 中怎样把平方打在数字或符号的右上角，就像我们平时作业书写一样？不是" ^2 ".

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

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

在尝试求解把vk/vkdig这样的变上限积分携带自变量代入到微分方程中求解时，发现它的计算速度非常的慢。请问该怎么优化呢？

背景：我正在尝试求解一个不稳定性问题，它最终被我归结为含有一个未知参数w（也是我希望求解的值）的两个耦合微分方程组。由于这个问题的复杂性，许多传统方法（Shooting method,本征矩阵法）都不太合适。我目前想到的最好方法是这样的： ClearAll["Global`*"]; k = Sqrt[3/8]*wp/V0; wp = 1; wc = wp*V0/c; c = 1; V0 = 1/5; gamma0 := 1/Sqrt[1 - V0^2/c^2]; w0[k_] := N[Sqrt[wp^2/2/ gamma0^3*(Sqrt[1 + 8*k^2*V0^2/wp^2*gamma0^3] - 1 - 2*k^2*V0^2/wp^2*gamma0^3)]]; Omega[V_] := w - k*V; A[V_] := (w - k*V)^2*(w^2

lz没有系统性的学过积分方程，想用mathematica求解积分方程，查阅帮助文档后发现，mathematica好像没有关于积分方程的文档，百度了一下，一维积分方程可以用Sovle求解，但对于我的问题却失效了，以下是代码： eqn = {Integrate[f[x]/(2\[CenterDot]\[Pi]) 1/(x1 - x), {x, 0, c}] == Subscript[V, \[Infinity]] \[Alpha], f[0] == Infinity, f[c] == 0} Solve[eqn, f[x]] 求问mathematica怎么求解积分方程qvq

为什么会显示“无法求得一个明确的导数公式，考虑使用选项设置SolveDelayed->Ture”.这个选项设置在哪里？我没有找到。我的微分方程编辑逻辑应该没问题，在matlab里面能算出数据，但是不理想，来mathmatical算结果直接出错，大佬们帮我看看^_^。 代码如下： m = 0.15; g = 10; a1 = (2.5^0.5)*10^-1; a2 = (50^0.5)*10^-2; a3 = 10^-3; I1 = 2.019*10^-4; I2 = 4.481*10^-4; I3 = 6.*10^-4; I12 = -2.1706*10^-6; I21 = -2.1706*10^-6; \!\(\* TagBox[ RowBox[{"II", "=", RowBox[{"{", RowBox[{ RowBox[

我算了许多omega，我现在想要提取里面的实部在0~10，虚部在-10~0的，并且按照虚部大小排列，应该怎么写？数据如下 {-178.676 + 3100.58 I, 178.676 + 3100.58 I, -676.972 + 2827.62 I, 676.972 + 2827.62 I, 1091.23 + 2497.01 I, -1091.23 + 2497.01 I, -1373.58 + 2118.98 I, 1373.58 + 2118.98 I, -1544.1 + 1732.84 I, 1544.1 + 1732.84 I, 1444.31 + 1650.16 I, -1444.31 + 1650.16 I, 159.172 + 1645.15 I, -159.172 + 1645.15 I, 1117.67 + 1644.89 I, -1117.67 + 1644.89 I, 477.962 + 1644.62 I, -477.962 + 1644.62 I, 797.587 + 1644.32 I, -797.587 + 1644.32 I, 1717.78 + 1529.86

求！ a-y/x-y 是关于θ的值，要画的也是θ=a-y/x-y 的三维平面 <math xmlns='http://www.w3.org/1998/Math/MathML' mathematica:form='StandardForm' xmlns:mathematica='http://www.wolfram.com/XML/'> <mrow> <mi>Show</mi> <mo>[</mo> <mrow> <mrow> <mi>Plot3D</mi> <mo>[</mo> <mrow> <mrow> <mi>NF</mi> <mo>/.</mo> <mi>vall11</mi> </mrow> <mo>,</mo> <mrow> <mo>{</mo> <mrow> <mi>&#952;</mi> <

如题，想要实现“矩阵直乘再乘积等于对应乘积的直乘”运算，可以通过定义 A_\[CircleTimes]B_ := KroneckerProduct[A, B]; kroneckerProductRule = (a_\[CircleTimes]b_) . (c_\[CircleTimes]d_) :> \ (a . c)\[CircleTimes](b . d); 实现。现在想将这个运算推广到多个矩阵，也即 (A \[CircleTimes] B \[CircleTimes]C) . (D \[CircleTimes] E \[CircleTimes] F) = （A . D) \[CircleTimes] (B . E) \[CircleTimes] (C . F)，尝试模式匹配后发现不成功，甚至两项矩阵都未能实现想要的结果。请问有没有大佬能解答一下。代码如下：

\[CapitalOmega] = ImplicitRegion[True, {{z, 0, 200}}]; \[CapitalGamma]1 = NeumannValue[0 , z == 0 || z == 200]; TraditionalForm[\[CapitalGamma]1] {nO, pO} = 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] + 1 - E^(-2 u) + p[t, z] - n[t, z]) n[t, z] + \[CapitalGamma]1, D[p[t, z], t] == 0.1 D[p[t, z], {z, 2}] + \[Tau]/ nb \[Alpha] (1 - R) Pde

理论上来说,第一类曲面积分对1积分的结果应该是曲面的面积 然而今天计算x^2+y^2==Cos[z]^2的函数时,其积分结果总为0 同时尝试了一下x^2+y^2==z^2,积分结果是正常的 然后又尝试了一下x^2+y^2==(Arcsin[Sin[x]])^2,积分结果仍然是0 即使Region输出的结果貌似是完全相同的 Clear["Global`*"] surfacea = ImplicitRegion[ x^2 + y^2 == (Cos[z])^2 && 0 < z < \[Pi]/2, {x, y, z}]; surfaceb = ImplicitRegion[ x^2 + y^2 == (ArcSin[Sin[z]])^2 && 0 < z < \[Pi]/2, {x, y, z}]; surfacec = ImplicitRegion[x^2 +

如题。

求助，各位大佬，程序中间有没有错误，现在的情况是没有办法出图 代码如下 中间要对特征值进行处理，最后画图，原本的问题是特征值排序的问题，我用了一位大佬解答的程序片段，现在的情况是没有办法出图，不清楚是哪里出了问题，结果是后面的图片 p = K0^2 \[Omega] - Ke^2 \[Omega] + Kf^2 \[Omega] - Kv^2 \[Omega] + \[Epsilon] \[Omega]^2 - \[Omega]^3 - 2 Ke Kf Subscript[l, B] + 2 K0 Kv Subscript[l, B] + \[Epsilon] \!\(\*SubsuperscriptBox[\(l\), \(B\), \(2\)]\) - \[Omega] \!\(\*SubsuperscriptBox[\(l\), \

\[lbk]Omega[rbk]m = \[lbk]Pi[rbk]*4.3*10^5; m = 1.45*10^-7; k1 = k2 = \[lbk]Omega[rbk]m; \[lbk]CapitalDelta[rbk] = 1.2*\[lbk]Omega[rbk]m; \[lbk]CapitalDelta[rbk]2 = 2*\[lbk]Omega[rbk]m; P = 2*10^5; e = 1.6*10^-19; J = 1.5*k1; \[lbk]CurlyTheta[rbk] = 0; G = 0; \[lbk]HBar[rbk] = 1.05*10^-34; L = 0.025; \[lbk]Omega[rbk]L = (6*\[lbk]Pi[rbk]*10^8)/(1.064*10^-6); g = (\[lbk]CapitalDelta[rbk] + \[lbk]Omega[rbk]L)/(L - Q); Q = (\[lbk]HBar[rbk]*g*Abs[lbk]A[rbk]^2)/(m*\[lbk]Omega[rbk]m^2); \[lbk]Omega[rbk]1 = ((\[lbk]CapitalDelta[rbk] + \[lbk]Omega[rbk]L)*L)/(L - Q); \[lbk]CurlyEpsilon[rbk] = Sqrt[lbk](

各位大佬，我在计算特征值时，出现错序 跳变问题，直接画图就能出现。两个都出现了，应该是同一个问题，劳烦各位大佬解答，代码不长，只是占的页数比较大，请见谅。 程序一 \[Gamma] = 1; Subscript[\[Gamma], c] = 0*\[Gamma]; Subscript[\[Gamma], d] = 2*\[Gamma]; Subscript[\[Gamma], a] = 0*\[Gamma]; Subscript[\[Gamma], b] = 2*\[Gamma]; k = 1; J = 1; \[CapitalDelta] = 0.01; \[CapitalPhi] = 2*\[Pi]; n = 0;(*AW=0;*) U = N[BesselJ[n, AW]]; W = N[BesselJ[n, AW]]; TE = 20; e1 = Collect[Eigenvalues[1/I ({ {-I*Subscript[\[Gamma], a]/2, -k

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])) t = (3*\[Alpha] + w*(8 - \[Alpha]) + d*(4 - \[Alpha]) - c*(4*k - 4 - \[Alpha]) - 4)/(4 - 3*\[Alpha]) Plot[{(2*w + \[Alpha]*(1 + c) + 2*k*c*(\[Beta] - 1))/(\[Alpha]*(4 - \[Alpha])), (2* w + \[Alpha] + \[Alpha]*c + \[Alpha]*t - 2*c*k)/(\[Alpha]*(4 - \[Alpha]))}, {w, 3.7, 4.8}] 这段代码有啥错误吗，我的输出（白色底）和原文（绿色底）不一样

Clear["Global`*"] Clear@Derivative DSolve[-g m + k 2 x[t] - 2 k l x[t]/Sqrt[l^2 + x[t]^2] + m (x^\[Prime]\[Prime])[t] == F[t], x[t], t] 求解微分方程的代码，但是运行后却直接把原式还给我了？不知道怎么回事？求大佬教教

我用mathematica13.0计算一个很简单的极限并要求给出step by step solution, 发现答案没错但给的过程有明显错误（划线处），怎么回事？ Limit[(2*Sin[p*Pi] - 2*p*Pi*Cos[p*Pi])/p^2, p -> Plus[Infinity]]

画一个三维的箭头，z方向的范围比x，y方向大很多，默认比例画出来的箭头形状是没问题的 Graphics3D[{Red, Arrowheads[0.1], Arrow[Tube[{{1, 1, -1}, {2, 2, 0}, {3, 3, -1}, {4, 4, 20}}, 0.05]]}] 但是不能展示比例这么奇怪的图，要把z方向压缩一下，这样就改变了箭头的形状，变得很奇怪 Graphics3D[{Red, Arrowheads[0.1], Arrow[Tube[{{1, 1, -1}, {2, 2, 0}, {3, 3, -1}, {4, 4, 20}}, 0.05]]}, BoxRatios -> {1, 1.5, 1}] 好像没有三维箭头设置长宽比的选项，只能改大小，请问怎么让下面的图片中箭头形状变

代码是 Clear["Global`*"] m := 1.320; \[Sigma] := 0.200; \[Alpha]s := 0.35; \[Mu]R := m/4; \[Delta] := 10^-10; V[r_] := -4/3*\[Alpha]s/r + \[Sigma]*r; Veff[r_] := 2*\[Mu]R*V[r]; SolveRadialEq[\[CurlyEpsilon]_] := NDSolve[{u''[r] == (Veff[r] - \[CurlyEpsilon]^2/4 + m^2) *u[r], u[0] == \[Delta], u'[0] == 1}, u, {r, 0, 20}]; CheckBoundaryCondition[sol_, \[CurlyEpsilon]_] := Module[{uSol, value}, uSol := u[r] /. First[sol]; value := uSol[20]; value]; FindEnergyLevel[l_, \[CurlyEpsilon]Min_, \[CurlyEpsilon]Max_, tolerance_] := Module[{\[CurlyEpsilon]Mid, sol, check}, While[

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

把求根公式换成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 =

SetAttributes[{a, b, q}, Constant]; eqn = {Laplacian[w[x, y], {x, y}] == v[x, y], Laplacian[v[x, y], {x, y}] == 7}; bcs = {w[0, y] == 0, w[x, 20] == 0, w[10, y] == 0, w[x, 0] == 0, Derivative[1, 0][w][0, y] == 0, v[x, 20] == 0, Derivative[1, 0][w][10, y] == 0, v[x, 0] == 0}; ufun = NDSolveValue[{eqn, bcs}, w, {x, 0, 10}, {y, 0, 20}] NDSolveValue::fembdnl: The dependent variable in (w^(1,0))[0,y]==0 in the boundary condition DirichletCondition[(w^(1,0))[0,y]==0,x==0.] needs to be linear. 怎么搞成线性的啊

想请问一下吧友我想用mathematica的dsolve解微分方程组，我带入的参数（如Y（T））都是矩阵形式的，但是dsolve一直不输出，想请问下如何解决，是不是必须一条一条把方程列出来啊。

无论是内置的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，并且还会随着非对角元

\[Phi]1 = D[\[Phi][x, y, t], x, x, t] \[Phi]2 = D[\[Phi][x, y, t], x, x, x] \[Phi]3 = D[\[Phi][x, y, t], x, x, y] \[Phi]4 = D[\[Phi][x, y, t], y, y, t] \[Phi]5 = D[\[Phi][x, y, t], y, y, x] \[Phi]6 = D[\[Phi][x, y, t], y, y, y] \[Phi]7 = D[\[Phi][x, y, t], t] \[Phi]8 = D[\[Phi][x, y, t], x] \[Phi]9 = D[\[Phi][x, y, t], y] \[Phi] = \[Psi]0[y, t] + \[Psi] \[Phi]1 + \[Phi]4 - \[Lambda]^2 \[Phi]7 - \[Phi]9*\[Phi]2 + \[Phi]9*\ \[Phi]5 + \[Lambda]^2*\[Phi]9*\[Phi]8 + \[Phi]8*\[Phi]3 + \[Phi]8*\ \[Phi]6 - \[Lambda]^2*\[Phi]8*\[Phi]9 + \[Beta]*\[Phi]8 == 0 TraditionalForm[\[Phi]1 + \[Phi]4 - \[Lambda]

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

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

最近在重看有限元相关的东西，结果在看到等参元这部分内容的时候发现了一个数学上的小问题，姑且拿出来水一帖。考虑下面这个参数方程： 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，其余