We formulate a Least Squares function.

Input := 

h[a_,b_] = Sum[(f[data[[i]][[1]]] - data[[i]][[2]])^2,
		{i,1, Length[data]}];

We can look at several contour plots to get an idea of where minimum values of least square function could be.