一点都不乱，是这两个

   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

你要用的是第一个，它调用第二个

[此贴子已经被作者于2008-12-3 23:45:17编辑过]