function [z] = anorm(alpha)
% function [z] = anorm(alpha)
%
% approximate area under the standardized normal density function
% taken from coarse resolutiontable by cubic interpolation
%
% input  :	alpha		- single sided area in %
%
% output :	z		- integration border to reach alpha
%
% version 0.1.0		last change 14.05.1998

% G.Krahmann, LODYC Paris, May 1998

% hand typed !!! table 
tab = [0.0	,	0.5000;...
       0.1	,	0.4602;...
       0.2	,	0.4207;...
       0.3	,	0.3821;...
       0.4	,	0.3446;...
       0.5	,	0.3085;...
       0.6	,	0.2743;...
       0.7	,	0.2420;...
       0.8	,	0.2119;...
       0.9	,	0.1841;...
       1.0	,	0.1587;...
       1.1	,	0.1357;...
       1.2	,	0.1151;...
       1.3	,	0.0968;...
       1.4	,	0.0808;...
       1.5	,	0.0668;...
       1.6	,	0.0548;...
       1.7	,	0.0446;...
       1.8	,	0.0359;...
       1.9	,	0.0287;...
       2.0	,	0.0228;...
       2.1	,	0.0179;...
       2.2	,	0.0139;...
       2.3	,	0.0107;...
       2.4	,	0.00820;...
       2.5	,	0.00621;...
       2.6	,	0.00466;...
       2.7	,	0.00347;...
       2.8	,	0.00256;...
       2.9	,	0.00187];

% exit, if out of borders
if ~isempty(find(alpha<0.187))
  warning('Alpha too small.')
end

% interpolate table
z = interp1(tab(:,2),tab(:,1),alpha/100,'cubic');
