[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]])

I read that it uses a cubic polynomial to interpolate, can someone provide
me with the process to calculate these coefficients?

Thanks,

On Wed, Apr 29, 2020 at 5:10 PM Travis Czechorski <tjczec01 at gmail.com>
wrote:

> 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/dbcbd426/attachment.html>