y = 0.3*exp(-2.5*x) + 1 + randn(size(x))/50;
ft = fittype( 'a*exp(-b*x)+c', 'independent', 'x' );
opts = fitoptions( 'Method', 'NonlinearLeastSquares','start',[1 1 1]);
[fitresult, gof] = fit(x,y, ft, opts);
Now, look at what fit returns.
Can we extract what you want now? Of course.
How would you get things like confidence intervals on the parameters? When you don't know how to interact with an object, try this:
Do you see anything that might be useful? How about confint?