[feinresult,t]=fein(t);

function [feinresult,t]=fein(t);
feinresult=[];
persistent aln ; 

;
% common :: ;
global feinx_1; if isempty(feinx_1), feinx_1=0; end;
global feinx_2; if isempty(feinx_2), feinx_2=0; end;
global feinx_3; if isempty(feinx_3), feinx_3=0; end;
%% common /feinx / x , a , fkm;
%% common /feinx / feinx_1 , feinx_2 , feinx_3;
if isempty(aln), aln=0; end;
aln =(feinx_3-t).*feinx_1 - feinx_2.*log(t);
feinresult = 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 fein