[b, Ib, e, s, Is] = lsfit(y,x,C)
Input : y dependant variate (column vector)
x regressor variates
C confidence level for the confidence intervals
(default0.95).
Output : b vector of point estimates,
Ib confidence intervals,
s estimated standard deviation of error
with confidence interval Is,
e residuals.
An intercept not included in x is automatically
added as an extra column.
See also LINREG and LSSELECT.