[sisyphus] maxima && gnuplot && sbcl

[sisyphus] maxima && gnuplot && sbcl

[Rinat Bikov]

Здравствуйте, уважаемые!
На Maxima 5.18.1 пытаюсь выполнить следующую последовательность действий:
showtime:true$
C1:125000$
C2:15$
Q:2000$
P:1/4$
F(b):=integrate(%e^(-s^2/(2*Q^2)),s,minf,b)$
f(x):=((1-P)*C1*(1-F(x)) + P*C2*Q*(x/2 + 2*%pi/4) * F(x/2 - sqrt(2*%pi)/4));
plot2d(f,[x,-3000,-2990]);

После чего maxima зависает на целый час, sbcl грузит процессор на 100 процентов.
Потом выдаётся сообщение:

plot2d: expression evaluates to non-numeric value everywhere in plotting range.
plot2d: nothing to plot.
Evaluation took 3140.9800 seconds (3539.6040 elapsed) using 23500.930 MB.
(%o8)                                false

График соответственно не строится.
Это баг или со float(f(x)) должно работать? :).
Собственно, сейчас жду ещё один часок, чтобы узнать ответ :).

P.S. В математике столько ждать не нужно для построения графика
этой функции, там считается за секунды. Неужели максима настолько крива?..

-- 
С уважением, Ринат Биков.

Re: [sisyphus] maxima && gnuplot && sbcl

[Rinat Bikov]

19 мая 2009 г. 20:19 пользователь Rinat Bikov написал:
> Это баг или со float(f(x)) должно работать? :).
> Собственно, сейчас жду ещё один часок, чтобы узнать ответ :).
Неа, не вышел каменный цветок...
f1(x):=float(f(x));
plot2d(f1,[x,-3000,-2990]);
plot2d: expression evaluates to non-numeric value everywhere in plotting range.
plot2d: nothing to plot.
Evaluation took 6569.1550 seconds (7225.5440 elapsed) using 45428.142 MB.

В upstream стоит писать об этом или это норма? :)

-- 
С уважением, Ринат Биков.

Re: [sisyphus] maxima && gnuplot && sbcl

[Rinat Bikov]

19 мая 2009 г. 21:41 пользователь Rinat Bikov написал:
> plot2d(f1,[x,-3000,-2990]);
Мда, plot2d(f(x),[x,-3000,2990]); работает :).

-- 
С уважением, Ринат Биков.

Re: [sisyphus] maxima && gnuplot && sbcl

[Alexey Borisenkov]

Rinat Bikov пишет:
> 19 мая 2009 г. 21:41 пользователь Rinat Bikov написал:
>> plot2d(f1,[x,-3000,-2990]);
> Мда, plot2d(f(x),[x,-3000,2990]); работает :).

Я, конечно, не совсем по теме, но, кажется, вашу функцию F(x) несложно 
получить из имеющейся в наличии функции erf(x)? Или это слишком простой 
путь и надо именно через интегрирование? :)

Re: [sisyphus] maxima && gnuplot && sbcl

[Rinat Bikov]

19 мая 2009 г. 21:53 пользователь Alexey Borisenkov написал:
> Я, конечно, не совсем по теме, но, кажется, вашу функцию F(x) несложно
> получить из имеющейся в наличии функции erf(x)? Или это слишком простой путь
> и надо именно через интегрирование? :)

Угу, можно получить, вот она:
7500*(x/2+%pi/2)*(1000*sqrt(2)*sqrt(%pi)-1000*sqrt(2)*sqrt(%pi)*erf((sqrt(%pi)-sqrt(2)*x)/8000))+93750*
(-1000*sqrt(2)*sqrt(%pi)*erf((sqrt(2)*x)/4000)-1000*sqrt(2)*sqrt(%pi)+1)
Только вот задание - найти минимум и построить график именно для той
функции, которую я приводил.


-- 
С уважением, Ринат Биков.