A fast algorithm for S-estimation in robust heteroscedastic regression