Oh, thanks. [If we could parenthesize constants, I guess it would be that 1p1j1000 is 1p(1j1000) rather than (1p1)j1000. And, that is documented at http://www.jsoftware.com/help/dictionary/dcons.htm - though I wonder what the advantage is.]
As for the closed form representation issue, I'm not sure whether that should be significant or not. Still, it seems to be the way things are right now. Thanks again, -- Raul On Thu, Jul 20, 2017 at 11:18 AM, HenryRich <[email protected]> wrote: > i: is not restricted to Gaussian integers: > > i: 1p1 j. 1000 > > _3.14159 _3.13531 ... > > > 1p1j1000 is not what you think it is (not what I thought it was either): > > 1p1j1000 > > 1.16714j2.91674 > > 1p1 ^ 1j1000 > > 1.16714j2.91674 > > > We'd better get Martin's Constants page finished! > > > (1&o.@*:) is not in the domain of d., but I wouldn't call that a'd. issue' - > there is no closed-form integral for this function. > > > Henry Rich > > > > > > On 7/20/2017 2:16 AM, Raul Miller wrote: >> >> load'plot' >> plot ((2&o.@*:t.i.100)&p.d._1;(1&o.@*:t.i.100)&p.d._1) 1p1*i:1j1000 >> >> This has two domain error workarounds which might interest people. >> >> First, i:1p1j1000 does not work, but 1p1*i:1j1000 does. In other >> words, i: has a gaussian integer domain (perhaps because that is >> simpler to document) even though it might make sense to allow the real >> part to be non-intger when the complex part is non-zero. >> >> Second, expressions such as 1&o.@*:d._1 do not work, but >> (1&o.@*:t.i.100)&p.d._1 does. In other words taylor series polynomials >> can sometimes be a workaround for d. domain issues. >> >> (Finally, note that plot ;/|:+.(^@j.@*:t.i.100)&p.d._1] 1p1*i:1j1000 >> gives a similar plot, but that transform is quirky for me to explain.) >> >> Thanks, >> > > > > --- > This email has been checked for viruses by AVG. > http://www.avg.com > ---------------------------------------------------------------------- > For information about J forums see http://www.jsoftware.com/forums.htm ---------------------------------------------------------------------- For information about J forums see http://www.jsoftware.com/forums.htm
