heat map with ODEs solutions

Hi every body, I'm trying to map the solution universe of a set of differential equation which I've optimized using lsqcurvefit to an objective function and i want to plot the simulation results obtained when sweeping unknown parameters as a heat map accordingly to their optimal fit to my experimental data. Thanks in advance.