We study the curvaton scenario using gauge-invariant second order perturbation theory and solving the governing equations numerically. Focusing on large scales we calculate the non-linearity parameter f_nl in the two-fluid curvaton model and compare our results with previous analytical studies employing the sudden decay approximation. We find good agreement of the two approaches for large curvaton energy densities at curvaton decay, Omega_dec, but significant differences of up to 10 percent for small Omega_dec.