fit with least squares and store outliers

Hello! I have a vector called "time" and i want to find all the outliers after fitting with a least squares model. The outliers will be all the numbers of my vector that have a distance greater than 3 standard deviations. I want the outliers stored in a new vector. For example, in the code below an obvious outlier is 24th number of my vector, as I can see from the scatter plot too. I really know nothing about curve fitting so detailed answers will be appreciated. Thanks on advance.

time = [5.5869 12.5278 4.2830 6.8060 6.5090 10.7813 23.9145 12.2830 8.0845 6.3399 8.2545 8.5110 13.8448 10.5343 7.5985 4.8459 5.4595 9.5206 14.2282 40.1542 10.4116 12.4295 13.0836 152.6972 21.6085 24.9892 9.9827 7.4979 16.3898 15.8891 11.7220 15.2996 10.6531 6.4329 11.3927 7.7482 5.4688 5.1611 5.9166 7.1803 13.6313 8.0200 7.0225 9.3118 4.9429];

xAxis = 1: length(time);

yAxis = time;

scatter(xAxis, yAxis)

refline; %Superimpose a least-squares line on the scatter plot.