When you say "non-linear" do you really mean nonlinear in x? Or do you mean just a regular polynomial, which can be, say a cubic, which can be solved with regular linear least squares like polyfit() gives you? You can have a fit that's not a straight line and it's still linear least squares because it's linear in the coefficients, not the x terms. I attach a demo for you that uses polyfit().
