How to fit symetrical shifted functions from bvp4c with anonymous functions ?

Hi First, if you understood my question, congratulations as I'm not sure I understand it myself...

Here is the problem. I've got 1D profiles which look like gaussian I need to fit with output of bvp4c. Problem is that my profiles are almost symetrical but shifted along the x coordinate and that bvp4c only accept x = [0 ...].

I wanted to use kind of "multi-lines anonymous functions" and nlinfit to first fit the profile and then find the shift as a coefficient, but it seems you can't do multi-lines in anonymous functions...

Any idea?

Thanks a lot again if you can understand my problem and then fix it :-)