kolmcdf.m

来自「带延时的微分方程」· M 代码 · 共 29 行

function p = kolmcdf(x)
% KOLMCDF Kolmogorov cumulative distribution function (cdf).
%    P = KOLMCDF(X) computes the Kolmogorov cdf at the values in X.
% 
%    The size of P is the common size of X. A scalar input  
%    functions as a constant matrix of the same size as the other inputs.    
% Author: Sergiy Iglin
% e-mail: [e-mail]
% or: [e-mail]
% personal page: http://iglin.exponenta.ru

if nargin <  1, 
    error('Requires at least one input argument.');
end
p=zeros(size(x));
num=find(x>0);
xnum=x(num);
pnum=ones(size(num));
for k=1:1000,
   add=(-1)^k*exp(-2*k^2*xnum.^2);
   pnum=pnum+2*add;
   if norm(add,1)==0,
      break
   end   
end   
p(num)=p(num)+pnum;