function j=junk9(par,x) % fminsearch('junk9',[2,1.1/2],[],x) f1=lognpdf(x,1.1,.4); f2=gampdf(x,par(1),par(2)); j=max(abs(f1-f2));