PRO maxlik, img_name, psf_name, Niter, maxlik

rawim= double( readfits(img_name) )
help, rawim
print, total(rawim)

tit=string(img_name,psf_name,f='(A," maxlik deconvolved by ",A)')

function_2d,id,rawim, DATASET_NAME='rawim',WIDGET_TITLE=tit,/UNITY_ASPECT

psf=readfits(psf_name)
psf=psf / total(psf, /DOUBLE)


;Max_Likelihood,rawim,psf,maxlik,FT_PSF=psf_ft
for ii=0,Niter-1 do begin
  Max_Likelihood,rawim,psf,maxlik,FT_PSF=psf_ft
  if ((ii mod ceil(Niter/28.)) EQ 0) then $
    function_2d,id,maxlik, DATASET_NAME=string(ii)
endfor
function_2d,id,maxlik, DATASET_NAME=string(ii)
return
end