This question already has an answer here:
My questions is, that I see people using R´s lm() (linear regression model) with
Y ~ X^2 e.g. here: Simple non-linear regression problem
But I dont see how and why it works, hence lm() is clearly stated as a linear estimator (using OLS), why does it work for quadratic, clearly non-linear estimation, as well?