A numerical scheme for stochastic PDEs with Gevrey regularity