%% fedi.m is the Fermi-Dirac distribution function.
%%

function [fd,fdi]=fedi(cf,xv)

% cf = co-efficents for function
% xv = x-axis values
% fd = fermi Dirac values

fdi=1./(exp((xv-cf(3))./cf(4))+1);
fd=cf(1).*fdi+cf(2);
