FUNCTION fit_function,p

common func_xy, x, y

;The function we are trying to fit is
;y=((d+(k0/kl)^(1/(b-4))*D)/po))^(b-2)
;
;p(0)=ko
;p(1)=kl
;p(2)=b
;p(3)=po

a=1./(p(2)-4)
b=p(2)-2
k=float(p(0)/p(1))

;print,p(1)
print,max(abs(y-(((x+(k^a)*.17)/p(3))^b*$
                  (5e-7/(x+(k^a)*.17))^2)))

return,total(abs(y-p(1)*(((x+(k^a)*.17)/p(3))^b*$
	          (5e-7/(x+(k^a)*.17))^2)))^2

end