How to speed up the symbolic calculation of the Jacobian?

Hello,

I am trying to solve a system of multivariate nonlinear equations using the Newton iteration method. The equations contain symbolic parameters. The goal is to have the variables expressed as functions of the symbolic parameters. The problem is that the Jacobian calculation using the jacobian(f,v) syntax is very slow due to the high non-linearity of the equations to be solved. I have also attempted to use jacobianest(fun,x0), but calculating the Jacobian numerically won't work because I have symbolic parameters in the equations. How could I speed up the symbolic Jacobian calculation in this case? Thanks in advance.