%
%	Inversion by Non-linear iterative algorithm; 
%
%	Usage:
%
%		[ad,vd]=nlia_original(scat,i_type)	
%
%	where scat is the FULLY CORRECTED data (after removing zscat and multiplying by dcal)
%       and i_type is instrument type,vd is volume distribution.  
%%
%	*** Needs files m32x32a.mat and m32x32b.mat in \matlab directory.***
%
function [ad]=nlia_(scat,m2, H)
ad = ones(size(scat));
hi=input('ready?')

for iter=1:50
   calc=m2*ad;
   for j=1:32
      ad(j)=ad(j)*sum(H(:,j).*scat./calc);
   end
end


