I'm not sure what is happening but I would guess that at some point the ^(1/2) gets turned into ^(0), and then your standard deviation goes from .06... to 1. I.e., it seems like maybe the preparser doesn't catch these nested loadings.
-M. Hampton On Jun 15, 1:54 pm, Pogon <vic...@saase.net> wrote: > Hi, > I'm not sure if it's a bug or it's me doing something wrong. > > I have two files: > > test1.sage containing nothing but > print numpy.random.normal(0,(2*0.0061*0.33)^(1/2),1) > > and > > test2.sage containing > load "test1.sage" > > I import numpy > sage: import numpy > > Now > sage: load "test1.sage" > returns values always smaller than 1 > thats the right distribution, the same i get when using the notebook- > interface > > but > sage: load "test2.sage" > very often returns values bigger than 1, > thats a whole different distribution > > My system is ubuntu-9.04-amd64 on Pentium Dual Core > sage-4.0.1 from 2009-06-06 --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to sage-devel-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
