Numerical methods for polymeric systems