realpart strangeness, was: besselarray bug, was: Bessel plotting problem
- Subject: realpart strangeness, was: besselarray bug, was: Bessel plotting problem
- From: Andrej Vodopivec
- Date: Fri, 13 Jun 2008 19:24:12 +0200
Some more examples:
(%i1) create_list(bessel_j(1/3, -x), x, 1, 10), numer$
(%i2) map(x->realpart(x)/imagpart(x), %)$
(%i3) identify(%);
(%o3) [1/sqrt(3),1/sqrt(3),1/sqrt(3),1/sqrt(3),1/sqrt(3),1/sqrt(3),1/sqrt(3),1/sqrt(3),1/sqrt(3),1/sqrt(3)]
(%i4) create_list(bessel_j(1/5, -x), x, 1, 10), numer$
(%i5) map(x->realpart(x)/imagpart(x), %)$
(%i6) identify(%);
(%o6) [sqrt(2*sqrt(5)+5)/sqrt(5),sqrt(2*sqrt(5)+5)/sqrt(5),sqrt(2*sqrt(5)+5)/sqrt(5),sqrt(2*sqrt(5)+5)/sqrt(5),sqrt(2*sqrt(5)+5)/sqrt(5),sqrt(2*sqrt(5)+5)/sqrt(5),sqrt(2*sqrt(5)+5)/sqrt(5),sqrt(2*sqrt(5)+5)/sqrt(5),
sqrt(2*sqrt(5)+5)/sqrt(5),sqrt(2*sqrt(5)+5)/sqrt(5)]
(%i7) create_list(bessel_j(1/6, -x), x, 1, 10), numer$
(%i8) map(x->realpart(x)/imagpart(x), %)$
(%i9) identify(%);
(%o9) [sqrt(3),sqrt(3),sqrt(3),sqrt(3),sqrt(3),sqrt(3),sqrt(3),sqrt(3),sqrt(3),sqrt(3)]
(%i10) create_list(bessel_j(1/8, -x), x, 1, 10), numer$
(%i11) map(x->realpart(x)/imagpart(x), %)$
(%i12) identify(%);
(%o12) [sqrt(2)+1,sqrt(2)+1,sqrt(2)+1,sqrt(2)+1,sqrt(2)+1,sqrt(2)+1,sqrt(2)+1,sqrt(2)+1,sqrt(2)+1,sqrt(2)+1]
The code for identify is here:
http://wxmaxima.sf.net/identify.zip
HTH,
Andrej