Dear reader, hallo Morgan I took a short look to you frm's about distributions. well that's to say only to the frm of the normal distribution. I hope you don't mind I changed it a bit. gaussian_distm (XAxis) {;standard distribution take real(p2)=1, the other variables 0 if(real(p4)==0) bailout=100 else bailout=real(p4) endif sigma = real(p1) o = imag(p1) c1 = 2*sqr(sigma) z = o, c = pixel/sqrt(pi*c1) : z = c*exp(-sqr(z-o)/c1) |z|<bailout } and not gaussian_distm2 (XAxis) {;standard distribution take real(p2)=1, the other variables 0.Real (p2) must never be 0 if(real(p4)==0) bailout=100 else bailout=real(p4) endif sigma =p2 mu =p1 z = 0, c =1/(sqr(2*pi)*sigma),scalingfactor=50, c1=scalingfactor*pixel: z = c*exp(-.5*((z-mu)/sigma)^2)*c1 |z|<bailout } notes: 1.The 50 is only there because then you get when using p1=(0,0) and sigma=(1,0), that's to say the standard normal distribution, directly something on the screen. If you don't put it in then you have only a part of the basin of attraction. 2.I don't know why you omit the fasctor .5 in the exonent of e. 3. I don't see any reason not to allow complex numbers for mean and variance. 4. I don't know ehy first squaring sigma and then take the squareroot of sigma Cheers, Jos Hendriks