Rhinoscript | Boy surface

http://mathworld.wolfram.com/BoySurface.html

Option Explicit
'2010년 8월 27일 금요일
'Script written by Sanghoon Yoon @ +plastic
'Script copyrighted by +plastic

main
Sub main()
 
  Dim PP, N, ST, i, u,v, X,Y,Z, A,B,C, arrPoints 
    PP=Rhino.Pi
    N=50
    ST=PP/(N-1)
 
  ReDim arrPoints(N*N-1)

    i=0
  For u=0 To PP Step ST
    For v=0 To PP Step ST
   
      A= Cos(u)*Sin(v) : B=Sin(u)*Sin(v) : C=Cos(v)

      X=0.5*(2*A^2-B^2-C^2+2*B*C*(B^2-c^2)+C*A*(A^2-C^2)+A*B*(B^2-A^2))
      Y=Sqr(3)/2*(B^2-C^2+(C*A*(C^2-A^2)+A*B*(B^2-A^2)))
      Z=(A+B+C)*((A+B+C)^3+4*(B-A)*(C-B)*(A-C))/5

      arrPoints(i) =Array(X,Y,Z)
      i=i+1
    Next
  Next

  Rhino.AddSrfPtGrid Array(N,N), arrPoints

End Sub