An automated tanh-function method for finding solitary wave solutions to non-linear evolution equations