how do i use the bootstrap method with a custom fit function?

i am trying to globally fit 17 double Gaussian derivative functions and then use the bootstrap method to analyze the error in the parameter fits. looking at the examples for [bootstat,bootsam] = bootstrp(...), i can't figure out how to implement my function for the @corr that works for the example given online.

the function i'm using is somewhat complicated: y0+A*((f/g)-B)*(exp(-((B-f/g)^2)/(2*(Bhyp^2+(S^2)*(f/g)^2))))/(Bhyp^2+(S^2)*(f/g)^2)+A*((f/g)-B)*(exp(-((B-f/g)^2)/(2*(Bhyp1^2+(S1^2)*(f/g)^2))))/(Bhyp1^2+(S1^2)*(f/g)^2)

but i should be able to label it as something and use it in the bootstrap function correct?