Hey @Antonio_Aguiar! You can reduce the cutoff_dim, but that isn’t ideal 
. Do you need the fock backend? The gaussian backend might be better here — it’s exact and less computationally intensive than an equivalent fock calculation. You would just need to do the following:
eng = sf.Engine('gaussian')
Although this isn’t a problem with your example, the gaussian backend cannot represent non-gaussian states 
. You can find out more information on the various backends available here: Circuits — Strawberry Fields 0.23.0 documentation