you can use UISTACK on the handles returned by HERRORBAR. Something along these lines might work for you:
h = errorbar(...)
With these handles you can also change the properties of the plot symbols and the error bars.
Evgheny - If you seriously need to fit such a surface in spherical coordinates, I might be able to help you out with a separate, un-posted toolbox that is capable of fitting a surface in spherical coordinates. You would need to contact me directly.
Thanks, function does do its job great!
Does anyone know alternative of this function for closed surface (sphere-like).
If I try this function with spherical coordinates, it works fine, but there are problems with the seam (line of joint)