We find the best fitting (in least squares sense) parameters a and b.

Input := 

sol = FindMinimum[h[a,b],{a,40},{b,.02}]

Output =

{0.00601426, {a -> 53.9126, b -> 0.0401633}}