A boundary integral method for two-dimensional (non)-Newtonian drops in slow viscous flow