% ‘½�€Ž®�i“ñŽŸŽ®�j‹ßŽ—

t = [1; 2; 3; 4; 5];
f = [4; 9; 2; 3; -3]; % a = 2; b = 4; c = -1;

A = [ ones(size(t)), t, t.*t ];
coefs = (A'*A)\(A'*f);

fprintf("%f %f %f\n",  coefs(1), coefs(2), coefs(3));
