I have a set of data points [z,y] , and I want to fit a function with two parameters to this data. I have done this as follows:
f=@(p,q) p(1)*((p(2) + q.^1).^-1);
p=lsqcurvefit( f, [1.4,1], z,y, [0.5,0.5],[3,3]);
But I want the area under the curve in the interval [0 1] to be a total of 1. Is it possible to build this in as a condition?