标签: maxima

我无法找到如何做到这一点:

x=a+b
y=c+d
z=x*y

我希望输出

z=ac+ad+bc+bd

不

z=xy

我试图告诉maxima有关hermite多项式的递归关系:

我的第一个表达是这样的:

phi[0]:exp(-1/2*x^2);

phi[1]:sqrt(2)*x*phi[0];

wxplot2d([phi[0],phi[1]], [x,-5,5]);

到目前为止一切顺利,但我现在要通过以下方式定义所有其他内容:

phi[n]:sqrt(2/n)*x*phi[n-1] - sqrt((n-1)/n)*phi[n-2];

这只是炸弹(堆栈溢出).我真的想这么说

wxplot2d(phi [10],[x,-5,5])会给我一个明智的图片吗？

如何告诉Maxima解决以下问题？("解决"部分):

我做了:

load(distrib);
fpprec: 100;
bftorat:true;
solve(2*bfloat(cdf_normal(x,0,1))-1=0.99999999999999999968130594071b0, [x]);
%,numer

拿到:

(%o1) "/usr/share/maxima/5.32.1/share/distrib/distrib.mac"
(%o2) 100
(%o3) true
`rat' replaced -1.99999999999999999968130594071B0 by -199999999999999999968130594071/100000000000000000000000000000 = -1.99999999999999999968130594071B0
`rat' replaced 5.0B-1 by 1/2 = 5.0B-1
`rat' replaced 5.0B-1 by 1/2 = 5.0B-1
`rat' replaced 7.071067811865475244008443621048490392848359376884740365883398689953662392310535194251937671638207864B-1 by 118807941462947422469655519336079782367473013592460/168019802134529020067676914738440478110633605571601 = 7.071067811865475244008443621048490392848359376884740365883398689953662392310535194251937671638207864B-1
`rat' replaced 7.071067811865475244008443621048490392848359376884740365883398689953662392310535194251937671638207864B-1 by 118807941462947422469655519336079782367473013592460/168019802134529020067676914738440478110633605571601 = 7.071067811865475244008443621048490392848359376884740365883398689953662392310535194251937671638207864B-1
(%o4) [x=
(168019802134529020067676914738440478110633605571601*inverse_erf(99999999999999999968130594071/100000000000000000000000000000))/118807941462947422469655519336079782367473013592460]
inverse_erf: inverse_erf(1.0) is undefined.
 -- an error. To debug this try: debugmode(true);

进一步尝试(看看理性替换是否影响inverse_erf):

inverse_erf(9.9999999999999999968130594071b?1);
gamma_incomplete: continued fractions failed for gamma_incomplete(5.0b-1, 1.675965338889773975600843228854238162008399514002414970690458801529039878850010279559673197036592113916729947593696740535214189646774061729913734402353264492788885098143556404059170138591463120333687838496039284224858192931635551067412157341539627014907074717352945374476804912353312948754404927014555821149803440423871160460111635311071245528519256957555845418916034380535359079516879576795825857468710891474077746896341697834315575814989209244705740463652472196503944998297956825510866851943203353716451062616549067258800559231646552924469724160521456041856694702333938138297284123098699530288993519920353577729741393726951293645734447176179175560311787410907660921783058989196733914852345086206761158383783714360001646880454976428470630991611033582649957934844195398077660932657131767136966243415075424193909691302431604307524134764326959890911322928113456129784004060990585972427799869290459688878023964671879763390452043333273689436077956597199441415992202082578463153853017929328667898523224007b0).
 -- an error. To debug …

在maxima中,有没有办法为子表达式应用变量替换？例如,替换x+ywith的实例z.

subst 适用于琐碎的案件,但不仅限于此.

(%i92) subst(x + y = foo, x + y);
(%o93) foo
(%i94) subst(x + y = foo, x + y + z);
(%o95) z + y + x

我有一个最大值方程，我想找到这个方程中变量的数值解。我想知道是否有人知道 maxima 中的哪个函数与 mathematica 中的 FindRoot 最相似？

似乎有几种方法可以做到这一点，但我希望根据 mathematica 代码检查结果，因此我对类似于 FindRoot 的东西感兴趣。

本

我有一个方程组，其形式为：

eq1 : q = a*x+b*y;
eq2 : r = c*x+d*y;

我希望它获得系数矩阵，在这个简单的情况下是：

T : matrix([a,b], [c,d]);

这样我就可以将方程组表示为：

eq3 : Q = T . X;

T因为我想在进一步的矩阵运算中使用。（原始方程是我从问题中得到的，我不想手动将它们转换成矩阵形式，因为这就是这种程序的用途。但我就是无法弄清楚）

我想要一个 Maxima 函数来检测已经定义的函数。这可以被称为is_function，因此is_function(sin)和is_function(exp)将为真，但is_function(%pi)和is_function(foo)将为假（假设我还没有定义foo）。这已经存在了吗？

手册说

一些 Lisp 函数隐藏在 Maxima 包中，即以下 ... functionp ...

我不确定这到底意味着什么，但functionp(sin)返回false，即使我定义了 ，functionp(foo)也会返回。另外，如果我像以前一样定义，则失败并显示 message 。falsefoo(x) := x;fundef(sin)dispfun(sin)fundef: no such function: sinfundef(foo)foo(x)

我有几个方程式,用于指定椭圆元素之间的关系(次要和主要半轴,偏心率,中心和中心距离,焦点参数).我想通过所有其他元素的组合找到表达这些元素的其他公式.solve()在某些组合上完全正常,但在另一组合中则不然.这些案件之间有什么区别我能如何解决？

    sys:[b=a*sqrt(1-e^2),p=a*(1-e^2),c=a*(1+e),d=a*(1-e)];
    solve(sys,[a,b,e,p]);
    solve(sys,[a,b,e,c]);
    solve(sys,[a,b,e,d]);
    solve(sys,[a,b,p,c]);
    solve(sys,[a,b,p,d]);
    solve(sys,[a,b,c,d]);
    /* solve(sys,[a,e,p,c]); -- gives error */
    /* solve(sys,[a,e,p,d]); -- gives error */
    /* solve(sys,[a,e,c,d]); -- gives error */
    solve(sys,[a,p,c,d]);
    solve(sys,[b,e,p,c]);
    solve(sys,[b,e,p,d]);
    solve(sys,[b,e,c,d]);
    solve(sys,[b,p,c,d]);
    solve(sys,[e,p,c,d]);

这是"好"和"坏"案例的输出:

(%i5) solve(sys,[a,b,c,d]);
                                        2
                  p           sqrt(1 - e ) p          p          p
(%o5)   [[a = - ------, b = - --------------, c = - -----, d = -----]]
                 2                 2                e - 1      e + 1
                e  - 1            e  - 1

(%i6) solve(sys,[a,e,p,c]);

algsys: tried and failed …

当我使用Eclipse或Visual Studio或Xcode时,项目的结构会自动创建,并且可以立即开始.

现在我正在尝试在Allegro CL express版本上创建一个项目,如果我收集所有lisp文件并运行它,它会显示包错误,我认为这是因为文件没有按特定顺序执行,因为必须在某个文件中创建包必须先执行.

我无法理解如何解决这个问题.关于在lisp中创建大型项目以及如何逐步处理其大小和复杂性的信息很少或很少.

我想将Maxima移植到Allegro CL上.我可以在这里获得一些帮助吗？如果可以解释多个文件和包中的代码分解以及基本上如何加载整个系统.

我知道lisp的基础知识,但我不了解lisp项目的结构.