This is very simple. Use matrices. I'll use your model as is... .

We don't need the epsilon since we are only fitting a curve.

Using your model the design matrix consists of three columns.

The first column is through , the second column is through and the last column is all 1's.

The vector Y consists of through .

The solution is the vector

where