     e.push norm(f)
 
     d = (sdot(e,i) for i in axes)
-    r = norm(d)
 
     if opts.traditionalLayout
       [y,z,x] = d


  s = axesCovariance.map Math.sqrt
  axes = transpose(axes)

  c = if degrees then 180/Math.PI else 1

  scales =
    upper: 1

         s[2]*c1]

    d = (sdot(e,i) for i in axes)

    if opts.traditionalLayout
      [y,z,x] = d
