You should use solver. Set up a cell that contains your total RMS error for all the sub-problems, and then get solver to minimise that. You may need to apply weightings to the individual error terms.
Goal seek won't work, it can only change one cell to drive another cell towards a target value. You need to change multiple cells.
Another option that might work better is to solve each cubic individually with solver, driving that from a macro.