The ODE solution methods covered our class were based on approximating the integral in the definition of the average slope using Newton-Cotes methods. This resulted in specific Runge-Kutta methods. Would it be possible, to derive Runge-Kutta methods based on Gaussian quadrature? Can you help me explain, conceptually, why or why not?
