[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.