%                             ncchi2_7.m
%  Scope:   This MATLAB macro computes the probability density function of the 
%           non-central chi-square distribution with seven degrees of freedom.
%  Usage:   xoutput = ncchi2_7(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) 1997-00 by LL Consulting. All Rights Reserved.

function   xoutput = ncchi2_7(x,lambda)

sum = 1./(1.875*sqrt(pi));
temp = sum;

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

a = 1./(8.*sqrt(2));
xoutput = a.*sum.*(exp(-0.5.*(x+lambda))).*x.*x.*(sqrt(x));