[sisyphus] maxima && gnuplot && sbcl

Rinat Bikov bikoz.r на gmail.com
Вт Май 19 19:19:34 MSD 2009


Здравствуйте, уважаемые!
На 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. В математике столько ждать не нужно для построения графика
этой функции, там считается за секунды. Неужели максима настолько крива?..

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


Подробная информация о списке рассылки Sisyphus