[SciPy-User] Interpolator coefficients for Radau method in solve_ivp

[Travis Czechorski]

Hello,

I am trying to figure out how to calculate the following coefficients in
the source code for the Radau integration method.

# Interpolator coefficients.P = np.array([ [13/3 + 7*S6/3, -23/3 - 22*S6/3,
10/3 + 5 * S6],
  [13/3 - 7*S6/3, -23/3 + 22*S6/3, 10/3 - 5 * S6],
 [1/3, -8/3, 10/3]])
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.python.org/pipermail/scipy-user/attachments/20200429/1f4efb35/attachment.html>