* [sisyphus] maxima && gnuplot && sbcl
@ 2009-05-19 15:19 Rinat Bikov
2009-05-19 16:41 ` Rinat Bikov
0 siblings, 1 reply; 5+ messages in thread
From: Rinat Bikov @ 2009-05-19 15:19 UTC (permalink / raw)
To: ALT Linux Sisyphus discussions
Здравствуйте, уважаемые!
На 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. В математике столько ждать не нужно для построения графика
этой функции, там считается за секунды. Неужели максима настолько крива?..
--
С уважением, Ринат Биков.
^ permalink raw reply [flat|nested] 5+ messages in thread
* Re: [sisyphus] maxima && gnuplot && sbcl
2009-05-19 15:19 [sisyphus] maxima && gnuplot && sbcl Rinat Bikov
@ 2009-05-19 16:41 ` Rinat Bikov
2009-05-19 16:48 ` Rinat Bikov
0 siblings, 1 reply; 5+ messages in thread
From: Rinat Bikov @ 2009-05-19 16:41 UTC (permalink / raw)
To: ALT Linux Sisyphus discussions
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 стоит писать об этом или это норма? :)
--
С уважением, Ринат Биков.
^ permalink raw reply [flat|nested] 5+ messages in thread
* Re: [sisyphus] maxima && gnuplot && sbcl
2009-05-19 16:41 ` Rinat Bikov
@ 2009-05-19 16:48 ` Rinat Bikov
2009-05-19 16:53 ` Alexey Borisenkov
0 siblings, 1 reply; 5+ messages in thread
From: Rinat Bikov @ 2009-05-19 16:48 UTC (permalink / raw)
To: ALT Linux Sisyphus discussions
19 мая 2009 г. 21:41 пользователь Rinat Bikov написал:
> plot2d(f1,[x,-3000,-2990]);
Мда, plot2d(f(x),[x,-3000,2990]); работает :).
--
С уважением, Ринат Биков.
^ permalink raw reply [flat|nested] 5+ messages in thread
* Re: [sisyphus] maxima && gnuplot && sbcl
2009-05-19 16:48 ` Rinat Bikov
@ 2009-05-19 16:53 ` Alexey Borisenkov
2009-05-19 16:58 ` Rinat Bikov
0 siblings, 1 reply; 5+ messages in thread
From: Alexey Borisenkov @ 2009-05-19 16:53 UTC (permalink / raw)
To: ALT Linux Sisyphus discussions
Rinat Bikov пишет:
> 19 мая 2009 г. 21:41 пользователь Rinat Bikov написал:
>> plot2d(f1,[x,-3000,-2990]);
> Мда, plot2d(f(x),[x,-3000,2990]); работает :).
Я, конечно, не совсем по теме, но, кажется, вашу функцию F(x) несложно
получить из имеющейся в наличии функции erf(x)? Или это слишком простой
путь и надо именно через интегрирование? :)
^ permalink raw reply [flat|nested] 5+ messages in thread
* Re: [sisyphus] maxima && gnuplot && sbcl
2009-05-19 16:53 ` Alexey Borisenkov
@ 2009-05-19 16:58 ` Rinat Bikov
0 siblings, 0 replies; 5+ messages in thread
From: Rinat Bikov @ 2009-05-19 16:58 UTC (permalink / raw)
To: ALT Linux Sisyphus discussions
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)
Только вот задание - найти минимум и построить график именно для той
функции, которую я приводил.
--
С уважением, Ринат Биков.
^ permalink raw reply [flat|nested] 5+ messages in thread
end of thread, other threads:[~2009-05-19 16:58 UTC | newest]
Thread overview: 5+ messages (download: mbox.gz / follow: Atom feed)
-- links below jump to the message on this page --
2009-05-19 15:19 [sisyphus] maxima && gnuplot && sbcl Rinat Bikov
2009-05-19 16:41 ` Rinat Bikov
2009-05-19 16:48 ` Rinat Bikov
2009-05-19 16:53 ` Alexey Borisenkov
2009-05-19 16:58 ` Rinat Bikov
ALT Linux Sisyphus discussions
This inbox may be cloned and mirrored by anyone:
git clone --mirror http://lore.altlinux.org/sisyphus/0 sisyphus/git/0.git
# If you have public-inbox 1.1+ installed, you may
# initialize and index your mirror using the following commands:
public-inbox-init -V2 sisyphus sisyphus/ http://lore.altlinux.org/sisyphus \
sisyphus@altlinux.ru sisyphus@altlinux.org sisyphus@lists.altlinux.org sisyphus@lists.altlinux.ru sisyphus@lists.altlinux.com sisyphus@linuxteam.iplabs.ru sisyphus@list.linux-os.ru
public-inbox-index sisyphus
Example config snippet for mirrors.
Newsgroup available over NNTP:
nntp://lore.altlinux.org/org.altlinux.lists.sisyphus
AGPL code for this site: git clone https://public-inbox.org/public-inbox.git