如上面的截图所示,我想要计算出能使F(x)=1/2的x值,代码如下:
- with(student):with(stats):
- e := 1:a := 1: z := 1: th_0 := 1:
- p := y->statevalf[pdf, normald](y) ;
- c := r->statevalf[cdf, normald](r);
- Inter := x -> c(sqrt(a+e)*(th_0-x*e/(a+e)-z*a/(a+e))) ;
- g := (xt, t) ->sqrt(e)*p(sqrt(e)*(xt-t)) ;
- f := (t, x) ->piecewise(th_0 < t, sqrt(a+e)*p(sqrt(a+e)*(t-x*e/(a+e)-z*a/(a+e)))/Inter(x), t <= th_0, 0) ;
- h := (x, xt)->int(f(t, x)*g(xt, t), t = -infinity .. infinity) ;
- F := x->int(h(x, xt), xt = x .. infinity) ;
- solve(F(x)-1/2);
可是为什么maple会报错呢?报错内容是:
Error, (in solve) cannot solve expressions with int(int(statevalf[pdf, normald](xt-t)*2^(1/2)*statevalf[pdf, normald](2^(1/2)*(t-(1/2)*x-1/2))/statevalf[cdf, normald](2^(1/2)*(1/2-(1/2)*x)), t = 1 .. infinity), xt = x .. infinity) for x
我刚刚接触maple而论文时间要求比较紧张,请大家帮帮忙,这是哪里出了问题,又该怎么改正,求指教!谢谢


雷达卡



京公网安备 11010802022788号







