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