高手帮忙看看这个VBA齿轮程序哪有问题（SolidWorks2007）

posinan

   经典图书

此例子来自张晋西编写的书。
用样条曲线拟合渐近线
Private Sub Userform_Initialize()
Me.Label1 = "齿 数"
Me.Label2 = "模 数"
Me.Label3 = "压力角"
Me.CommandButton1.Caption = "确 定"
Me.CommandButton2.Caption = "取 消"
'窗体上文本框赋初值
Me.TextBox1 = 21 '齿数
Me.TextBox2 = 12 '模数
Me.TextBox3 = 20 '压力角
End Sub
Private Sub CommandButton1_Click()
Dim points(9) As Double
Dim CZ As Double, CM As Double, CA As Double, CRa As Double
CZ = Me.TextBox1 '齿数
CM = Me.TextBox2 / 1000 '模数 ，/1000单位变为米
CA = Me.TextBox3 * 3.141 / 180 '压力角
'子程序计算出齿轮廓线的坐标points和顶圆半径cra
Call 齿轮廓线(CZ, CM, CA, points(), CRa)
Dim swapp As SldWorks.SldWorks
Dim swmodel As SldWorks.ModelDoc2
Dim nptdata(26) As Double
Dim vptdata As Variant
Dim swSkechSeg(1) As SldWorks.SketchSegment
Set swapp = Application.SldWorks
Set swmodel = swapp.ActiveDoc
Set swSkechSeg(0) = swmodel.CreateCircleByRadius2(0, 0, 0, CRa)
SwModle.InsertSketch2 True
nptdata(0) = -points(8): nptdata(1) = points(9): nptdata(2) = 0#
nptdata(3) = -points(6): nptdata(4) = points(7): nptdata(5) = 0#
nptdata(6) = -points(4): nptdata(7) = points(5): nptdata(8) = 0#
nptdata(9) = -points(2): nptdata(10) = points(3): nptdata(11) = 0#
nptdata(12) = points(0): nptdata(13) = points(1): nptdata(14) = 0#
nptdata(15) = points(2): nptdata(16) = points(3): nptdata(17) = 0#
nptdata(18) = points(4): nptdata(19) = points(5): nptdata(20) = 0#
nptdata(21) = points(6): nptdata(22) = points(7): nptdata(23) = 0#
nptdata(24) = points(8): nptdata(25) = points(9): nptdata(26) = 0# '样条曲线上的九个点
vptdata = nptdata
Set swSkechSeg(1) = SwModle.CreateSpline(vptdata) '创建齿轮廓线样条曲线
Dim bRet As Boolean
'绘制齿轮顶圆曲线
bRet = swmodel.CreatArcByCenter(0, 0, 0, points(8), points(9), 0, -points(8), points(9), 0)
swmodel.InsertSketch2 True
SwModle.ViewZoomtofit2 '整屏显示图形
End Sub
Sub 齿轮廓线(CZ As Double, CM As Double, CA As Double, points() As Double, CRa As Double)
Dim CR As Double, CRb As Double, CRf As Double, CSb As Double, Th(3)
CR = CM * CZ / 2 '齿轮分度圆半径
CRf = (CR - 1.25 * CM) '齿根圆半径
CRb = CR * Cos(CA) '齿轮基圆半径
CRa = CR + CM '齿轮顶圆半径
'齿轮基圆齿厚
CSb = Cos(CA) * (3.14 * CM / 2 + CM * CZ * (Tan(CA) - (CA)))
Th(1) = (3.14 * CM * Cos(CA) - CSb) / (2 * CRb)
Th(0) = Th(1) / 3
Th(2) = Th(1) + Tan(CA) - CA
'ACos---反余弦，自定义函数
Th(3) = Th(1) + Tan(Acos(CRb / CRa)) - Acos(CRb / CRa)
'第0点
points(0) = 0: points(1) = CRf
'第1点
points(2) = CRf * Sin(Th(0)): points(3) = CRf * Cos(Th(0))
'第2点
points(4) = CRb * Sin(Th(1)): points(5) = CRb * Cos(Th(1))
'第3点
points(6) = CR * Sin(Th(2)): points(7) = CR * Cos(Th(2))
'第4点
points(8) = CRa * Sin(Th(3)): points(9) = CRa * Cos(Th(3)) ‘各个圆与渐近线的焦点
'当基圆小于根圆，调整第1、第2点坐标，得到近似值
If CRb < CRf Then
'第1点
points(2) = points(6) * 0.2: points(3) = points(1) + 0.25 * CM * 0.03
'第2点
points(4) = points(6) * 0.7: points(5) = points(1) + 0.25 * CM * 0.8
End If
End Sub
Function Acos(X As Double) As Double '反余弦
Dim pi As Double
pi = 4# * Atn(1#) '45度=pi/4
If Abs(X) > 1# Then
MsgBox "cosX>1 ,Acos(X)函数出错", 1 + 16, "警告": Exit Function
Else
If Abs(X) = 1# Then
Acos = (1# - X) * pi / 2#
Else
Acos = pi / 2 - Atn(X / Sqr(-X * X + 1))
End If
End If
End Function
Private Sub CommandButton2_Click()
End
End Sub
Set swSkechSeg(0) = swmodel.CreateCircleByRadius2(0, 0, 0, CRa)这一行总是提示出现运行错误91，对象变量或with变量未设定。