Maxima似乎没有提出这个包含误差函数的方程的解析解。这里的自变量是“p”,要求解的因变量是“x”。
see an illustration of equation follow link
(%i3) solveexplicit:true$ ratprint:false$ fpprintprec:6$
(%i4) eqn: (sqrt(%pi)*(25*2^(3/2)*p-25*sqrt(2))*erf(1/(25*2^(3/2)*x))*x+1)/(25*p) = 0.04;
(%i5) solve (eqn, x);
(%o5) []
(%i6) eqn, [p=2,x=0.00532014],numer;
(%o6) 0.04=0.04
任何帮助或指向正确的方向表示赞赏。
据我所知,Maxima无法解决包含erf
的方程式。您可以通过find_root
获得数值结果:
(%i5) find_root (eqn, x, 0.001, 0.999), p=2;
(%o5) 0.005320136894034347
至于符号解决方案,我稍微研究了等式。人们可以把它变成erf(something/x)*x = otherstuff
形式,或者等同于erf(y) = somethingelse*y
,其中y = something/x
和somethingelse = otherstuff/something
如果我没有弄错的话。我对这种形式的方程式一无所知,但也许你可以找到一些东西。
是的,解决只能做多项式。我使用系列扩展来获得小的x值,并且精度足够好。
(%i11) seriesE: 1$
termE: erf(x)$
for p: 1 unless p > 3 do
(termE: diff (termE, x)/p,
seriesE: seriesE + subst (x=0, termE)*x^p)$
seriesE;
(%o11) -(2*x^3)/(3*sqrt(%pi))+(2*x)/sqrt(%pi)+1
但是,“表达式比配置设置允许的时间更长!”