[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