[dfeinresult,t]=dfein(t);

function [dfeinresult,t]=dfein(t);
dfeinresult=[];
persistent aln ; 

;
% common :: ;
global dfeinx_1; if isempty(dfeinx_1), dfeinx_1=0; end;
global dfeinx_2; if isempty(dfeinx_2), dfeinx_2=0; end;
global dfeinx_3; if isempty(dfeinx_3), dfeinx_3=0; end;
%% common /dfeinx/ x , a , fkm;
%% common /dfeinx/ dfeinx_1 , dfeinx_2 , dfeinx_3;
if isempty(aln), aln=0; end;
aln =(dfeinx_3-t).*dfeinx_1 - dfeinx_2.*log(t);
dfeinresult = exp(aln);
csnil=dbstack(1); csnil=csnil(1).name(1)~='@';
if csnil&&~isempty(inputname(1)), assignin('caller','FUntemp',t); evalin('caller',[inputname(1),'=FUntemp;']); end
end %function dfein