MATLAB: Kron efficiency

kron efficiency ode solver

Hi all, I have a matlab program whose bottleneck according to the profiler is the kron function. The program basically is a ode solver, but the dydt matrix is built by making extensive use of the kron function, so that the bottleneck of the ode solver is defining the dydt matrix.
Do you have any suggestion on how this could be optimized?
Thanks!