mathematica吧
关注: 19,720 贴子: 73,577

自带帮助是最好的教材,提问前请先仔细查查

  • 目录:
  • 程序设计及人工智能
  • 0
    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}]
  • 1
    请问各位前辈大佬,我有一个方程组,是由8个方程联立,但是未知数是7个,也就是8个方程构成的方程组求解7个未知量,请问用什么指令最好呢?个人试过DSolve、Nsolve、LeastSquares(最小二乘法)、Reduce以及广义逆矩阵的方法等求解,结果要么是零解,要么解不准确。但是我心里知道这个方程在虚数范围内是有解的。请问如何到底怎么求解呢?谢谢。
    Repentanze 12-23
  • 47
    请问大佬 怎么用mathematica画出跳棋棋盘?
  • 2
    各位大佬,我求解一个简单的微分方程,在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)
    000000_lp 12-23
  • 1
    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,
    Repentanze 12-23
  • 3
    有人遇到过这种问题,程序运行后,在重新打开,所有的代码前面都会添加上"In[1] / Out[1]"之类的字样,而且所有的Cell会合并成一个Cell。如图所示:第一张第二张分别是运行后,和关闭重新打开后的
  • 0
    症状: 顶楼中的 [ 和 ] 会被自动转换为 [lbk] 和 [rbk]。此现象在回复帖中不存在。 影响范围: 此bug至少在去年10月16日即已出现,至今(安卓12.58.1.0)未修。IPhone受影响情况不明,欢迎实验补充。 对策: 对于提问者,最简单的对策就是改用电脑访问贴吧网页版,不要使用使用手机客户端发帖。如果确有相关需求(所以说你们到底为什么会有这种需求?电脑没连网?那也可以开手机热点啊,访问网页又花不了什么流量),避免在顶楼张贴代码也是一个
    xzcyr 4-6
  • 17
    和Maple对照很明显,而且用Simplify也无法化简了
    草红样 5-3
  • 3
    如题。自编规则当然也不难,不过这里我想知道的是,能不能用Simplify,TrigExpand一类的化简型指令来得到这个结果。试了半天没试出来……是正好不支持这个化简呢,还是我没有找到?是我漏了什么重要假设吗?(我试验过的假设只有d属于Reals……)
    草红样 11-7
  • 5
    ArcTan[cos[170 D]/(cos[170 D] + 8/3 sin[40 D] sin[80 D])] 怎么得结果 ?
  • 2
    各位大佬,我的代码是这样的 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}] 结果得到非数值,这样应该如何解决?
  • 1
    输入: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]也没用。求助应该怎么解
    Repentanze 12-17
  • 2
    想请教一下,为什么不管是式子大于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
  • 1
    拓扑学中的一些拓扑空间,像莫比乌斯带,克莱因瓶这些,它们的图怎么用mathematica比较方便地画?或者说有什么软件能比较方便地画出这些图?
    liyajx 12-15
  • 7
    mathematica input 中怎样把平方打在数字或符号的右上角,就像我们平时作业书写一样?不是" ^2 ".
  • 7
    也算是对这帖(http://tieba.baidu.com/p/2964416898)的2,3,4条的解说吧。有什么意见或建议欢迎提出: http://note.youdao.com/share/?id=abd51087f44c0b6a41ff6022d549dc4
  • 58
    贴吧新弄了个贴条,贴条只能同时存在10个并且只能实际只能显示26字,于是我就在想着是不是该总结些Mathematica新手最常见的问题贴到那上面滚动播出(虽然很多人或许不会看,但有总比没有好嘛……)。试着弄了个草稿(二十戒),大家可以看看,哪些需要修改,哪些需要删节或补充: 二十戒 1 不要为了节约硬盘选择过老的版本。至少要用版本7 2  自带帮助是最好的教材,提问前先打开软件按下F1仔细查查 3 将光标停在不认识的函数前/中/后再
    Koko 3-25
  • 3
    在尝试求解把vk/vkdig这样的变上限积分携带自变量代入到微分方程中求解时,发现它的计算速度非常的慢。请问该怎么优化呢?
  • 4
    背景:我正在尝试求解一个不稳定性问题,它最终被我归结为含有一个未知参数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
    Repentanze 12-11
  • 4
    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
    汝之不惠 12-10
  • 2
    为什么会显示“无法求得一个明确的导数公式,考虑使用选项设置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[
  • 6
    我算了许多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
    xzcyr 12-7
  • 2
    求! 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> <
    xzcyr 12-7
  • 9
    如题,想要实现“矩阵直乘再乘积等于对应乘积的直乘”运算,可以通过定义 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),尝试模式匹配后发现不成功,甚至两项矩阵都未能实现想要的结果。请问有没有大佬能解答一下。代码如下:
    xzcyr 12-7
  • 2
    \[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
    xzcyr 12-7
  • 1
    理论上来说,第一类曲面积分对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 +
    xzcyr 12-7
  • 5
  • 4
    求助,各位大佬,程序中间有没有错误,现在的情况是没有办法出图 代码如下 中间要对特征值进行处理,最后画图,原本的问题是特征值排序的问题,我用了一位大佬解答的程序片段,现在的情况是没有办法出图,不清楚是哪里出了问题,结果是后面的图片 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\), \
    xzcyr 12-7
  • 1
    \[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](
    xzcyr 12-7
  • 3
    各位大佬,我在计算特征值时,出现错序 跳变问题,直接画图就能出现。两个都出现了,应该是同一个问题,劳烦各位大佬解答,代码不长,只是占的页数比较大,请见谅。 程序一 \[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
    xzcyr 12-7
  • 1
    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}] 这段代码有啥错误吗,我的输出(白色底)和原文(绿色底)不一样
    xzcyr 12-7
  • 4
    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] 求解微分方程的代码,但是运行后却直接把原式还给我了?不知道怎么回事?求大佬教教
    xzcyr 12-7
  • 1
    我用mathematica13.0计算一个很简单的极限并要求给出step by step solution, 发现答案没错但给的过程有明显错误(划线处),怎么回事? Limit[(2*Sin[p*Pi] - 2*p*Pi*Cos[p*Pi])/p^2, p -> Plus[Infinity]]
    xzcyr 12-7
  • 1
    画一个三维的箭头,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}] 好像没有三维箭头设置长宽比的选项,只能改大小,请问怎么让下面的图片中箭头形状变
    xzcyr 12-7
  • 3
    代码是 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[
    Repentanze 11-26
  • 27
    此帖是本吧初代吧主 @mm_酱 发在果壳小组的一篇教程。如今果壳小组已死,mm_酱近期又不太活跃,这里姑且重发一份: https://note.youdao.com/share/?id=058e6037396d925af1f4abe4d54a52 mm_酱你要是想整理下出个v2就吱一声。
    碳化物xj 11-29
  • 1
    把求根公式换成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 =
    xzcyr 9-7
  • 4
    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. 怎么搞成线性的啊
  • 4
    想请问一下吧友我想用mathematica的dsolve解微分方程组,我带入的参数(如Y(T))都是矩阵形式的,但是dsolve一直不输出,想请问下如何解决,是不是必须一条一条把方程列出来啊。
  • 26
    无论是内置的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,并且还会随着非对角元
  • 1
    \[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]
    liyajx 11-19
  • 4
    mma 同样作为 lisp-like 语言有 Throw/Catch 有 Dialog 甚至有 Goto 却没有 call/cc 的官方实现也是够奇怪的
    xzcyr 10-5
  • 23
    链接不发了。楼下简述更新内容
    DXDX450 11-13
  • 3
    最近在重看有限元相关的东西,结果在看到等参元这部分内容的时候发现了一个数学上的小问题,姑且拿出来水一帖。考虑下面这个参数方程: 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,其余
    xzcyr 8-3
  • 20
    其实本来想用的标题是《【交流】“这些问题你换软件也没用!”——受限于数学发展而无法或很难求解却经常有人想用Mathematica解的问题》。Mathematica作为一个软件,必然无法超越现有的数学理论,也就是说,对于尚未在理论上得到很好解决的问题,Mathematica是不太可能直接求解出来的——这话说出来大家都明白,但是,一般用户常常很难分辨哪些问题目前在理论上是很难或无法求解的,故开此交流帖,大家可以来谈谈自己所知道的相关问题。为避免
    氘化氢 11-30

  • 发贴红色标题
  • 显示红名
  • 签到六倍经验

赠送补签卡1张,获得[经验书购买权]

扫二维码下载贴吧客户端

下载贴吧APP
看高清直播、视频!

本吧信息 查看详情>>

小吧:小吧主共10

会员: Expr

目录: 程序设计及人工智能

友情贴吧