I'm attempting to construct a formula . However there's something amiss ; hence I'm seeking some help on this from the group. This is the formula : comment{ Series Circuit . (c) copyright 2016 sciwise@ihug.co.nz Edward.Q.Montague } DiodeRLC { f = 1/real(p1) z = pixel x = real(z) y = imag(z) p = 0.039757 h = 0.071743 s = -1 R = 1 L =0.001 C = 10^6 C = 1/C K = 10^23 K = (1.3806485279)/K twopi = 6.283185307179586 T = 270 q = 10^19 q = (1.602176620898)/q Id = 10^12 Id = 1/Id a1 = q/(K*T) b1 = (K*T)/(q*C*L) b2 = (K*T)/(Id*q*L) b3 = R/L t1 = 0: v = twopi*t1 v = 100*sin(w) q = s*a1*x q = exp( q ) u = y w = s*a1*y*y + s*b2*q*y + s*b3*y + b2*q*w + b1*q + s*b1 a = x+p*u b = y+p*w q = s*a1*a q = exp( q ) c = x+h*(u+b) d = y + h*(w + s*a1*b*b + s*b2*q*b + s*b3*b + b2*q*v + b1*q + s*b1) x = c y = d t1 = t1 + f z = x+flip(y) |y|