The SUSY-QCD β function to three loops