以下是引用cc_ruc在-10-27 22:53:00的发言：
函数很乱，好像不能直接用
一点都不乱，是这两个
function p = bvnl( dh, dk, r )
%
% A function for computing bivariate normal probabilities.
% bvnl calculates the probability that x < dh and y < dk.
% parameters
% dh 1st upper integration limit
% dk 2nd upper integration limit
% r correlation coefficient
function p = bvnu( dh, dk, r )
%
% A function for computing bivariate normal probabilities.
% bvnu calculates the probability that x > dh and y > dk.
% parameters
% dh 1st lower integration limit
% dk 2nd lower integration limit
% r correlation coefficient
你要用的是第一个，它调用第二个
