好的，谢谢您提醒！
Public Function MCoption(S, X, T, r, sigma, n, a)
'S：现价
'X：执行价
'T：到期日期间的时长
'r：连续复利的无风险利率
'sigma：股票价格的波动率
'epsilon：正态分布中的随机抽样
'n：模拟n次
'a：看涨看跌（看涨为1，看跌为0）
Dim i As Integer
Dim Se, P, sum As Variant

sum(0) = 0
For i = 1 To n Step 1
Se(i) = S * Exp((r - sigma * sigma / 2) * T + sigma * RndNorm(0, 1) * Sqr(T))
P(i) = Exp(-r * T) * Application.Max((X - Se(i)) * (-1) ^ a, 0)
sum(i) = sum(i - 1) + P(i)
Next i

MCoption = sum(n) / n

End Function

Public Function RndNorm(Mean As Double, Std As Double)
    Dim V1 As Double, V2 As Double, r As Double
    Do
        V1 = 2 * Rnd - 1
        V2 = 2 * Rnd - 1
        r = V1 ^ 2 + V2 ^ 2
    Loop Until r < 1
   
    RndNorm = Mean + V2 * Sqr(-2 * Log(r) / r) * Std
End Function

关于那个无法调用函数的问题，刚刚已经解决了，重新打开了一次又能用了，可能是打开启用宏之后需要重新打开文件？

Public Function MCoption(S, X, T, r, sigma, n, a)

Dim i As Integer
Dim sum, Se, P As Double

sum = 0

For i = 1 To n Step 1
Se = S * Exp((r - sigma * sigma / 2) * T + sigma * RndNorm(0, 1) * Sqr(T))
P = Exp(-r * T) * Application.Max((X - Se) * (-1) ^ a, 0)
sum = sum + P

Next i

MCoption = sum / n

End Function

Public Function RndNorm(Mean As Double, Std As Double)
    Dim V1 As Double, V2 As Double, r As Double
    Do
        V1 = 2 * Rnd - 1
        V2 = 2 * Rnd - 1
        r = V1 ^ 2 + V2 ^ 2
    Loop Until r < 1
   
    RndNorm = Mean + V2 * Sqr(-2 * Log(r) / r) * Std
End Function

你试试这个，我不太明白你为什么要将sum，Se和P设成array，设成标量就可以。而且我记得VBA里面array的size是要明确声明的，即便一开始不知道size用了variant，之后要用的时候还是要用ReDim来声明size。可能这就是一直报错的原因吧。

感谢大大！！！算出结果了！！

就是我之前说的“我记得VBA里面array的size是要明确声明的，即便一开始不知道size用了variant，之后要用的时候还是要用ReDim来声明size。可能这就是一直报错的原因吧。” 你的二叉树程序里面的vvec 有ReDim

其实这里你也不用redim，因为你已经知道了大小n因此直接定义Dim vvec As double (1 to n) 就行。我之前用Variant主要有两种情况：
一种是你要往array里装的变量类型未知，比如你既可能要装string也可能要装double
另一种是你在最初定义的时候还无法确定array的大小，比如你的array大小取决于后面计算出的某些量，这时候先定义成variant，不声明size然后等知道了确切的大小再redim size。

还得请教您：当n取到10,000还可以运行，渠道100,000就不行了，请问有解决的方法么？

Public Function MCoption(S, X, T, r, sigma, n, a)

Dim i As Long
Dim sum, Se, P As Double

sum = 0

For i = 1 To n Step 1
Se = S * Exp((r - sigma * sigma / 2) * T + sigma * RndNorm(0, 1) * Sqr(T))
P = Exp(-r * T) * Application.Max((X - Se) * a, 0)
sum = sum + P

Next i

MCoption = sum / n

End Function

把i定义成Long就行了

十分感谢大大！

三人行，必有我师。
学问，约等于debug。