A numerical study of non-gaussianity in the curvaton scenario