%                             ncchi2_2.m
%  Scope:   This MATLAB macro computes the probability density function of the 
%           non-central chi-square distribution with two degrees of freedom.
%  Usage:   xoutput = ncchi2_2(x,lambda)
%  Description of parameters:
%           x       - input, scalar independent variable
%           lambda  - input, non-centrality parameter
%           xoutput - output, computed probability density function 
%  Remark:  This macro can be used for the computation of the integral of the
%           probability density function.
%  Last update:  04/12/00
%  Copyright (C) 1998 by LL Consulting. All Rights Reserved.

function  xoutput = ncchi2_2(x,lambda)

sum = 1;
temp = sum;

for k = 1:49
   termk = temp.*lambda.*x./(4.*k.*k);
   sum = sum + termk;
   temp = termk;
end

xoutput = 0.5.*sum.*exp(-0.5.*(x+lambda));