Excellent! Many thanks.

I actually wish to fit a higher order polynomial, perhaps up to a cubic with cross terms, but the generalization to be made to your code should be straightforward.