Как воспользоваться результатом, который выдаёт функция scipy.integrate.solve_ivp?

import math
import numpy as np
import matplotlib.pyplot as plt
import sympy as sp
import scipy
 
delta=0.12e-3
n=3000
w1=2*math.pi*n
E=2.1e11
mu=0.3
ro=8e3
y1=np.array([1,0]) #для однородной системы
r_int=np.linspace(80e-3,640e-3,1200)
p2=40e6
y0=np.array([0,0]) #для неоднородной системы
r1=80e-3
r2=640e-3
h2=30e-3
def h(r):
    if r <= 160e-3 and r >= 80e-3:
        return 140e-3
    if r > 160e-3 and r < 600e-3:
        return 80e-3 - 13 * (r - 160e-3) / 88
    if r >= 600e-3 and r <= 640e-3:
        return 30e-3
def Fy_0(r,Y): #для неоднородного ДУ
    Fy1=-((1+mu)/r)*Y[0]+((1-mu**2)/(r*E*h(r)))*Y[1]
    Fy2=(E*h(r)/r)*Y[0]+((mu-1)/r)*Y[1]-ro*h(r)*(w1**2)*r
    Fy_0=np.array([Fy1,Fy2])
    return Fy_0
 
def Fy_1(r,Y): #для однородного ДУ
    Fy1=-((1+mu)/r)*Y[0]+((1-mu**2)/(r*E*h(r)))*Y[1]
    Fy2=(E*h(r)/r)*Y[0]+((mu-1)/r)*Y[1]
    Fy_ne0=np.array([Fy1,Fy2])
    return Fy_ne0
 
y_1=scipy.integrate.solve_ivp(Fy_1,(r1,r2),y1,first_step=0.5e-3,max_step=0.5e-3)
print("y_1:")
print(y_1)
y_0=scipy.integrate.solve_ivp(Fy_0,(r1,r2),y0,first_step=0.5e-3,max_step=0.5e-3)
print("y_0:")
print(y_0)
y21_r2= 9.276e+09
y20_r2= -1.891e+10
c1=(p2*h2-y20_r2)/y21_r2
print("c1= ",c1)
 
Выдаваемый результат:
 
y_1:
  message: The solver successfully reached the end of the integration interval.
  success: True
   status: 0
        t: [ 8.000e-02  8.050e-02 ...  6.400e-01  6.400e-01]
        y: [[ 1.000e+00  9.920e-01 ...  8.563e-01  8.563e-01]
            [ 0.000e+00  1.820e+08 ...  9.276e+09  9.276e+09]]
      sol: None
 t_events: None
 y_events: None
     nfev: 6727
     njev: 0
      nlu: 0
y_0:
  message: The solver successfully reached the end of the integration interval.
  success: True
   status: 0
        t: [ 8.000e-02  8.050e-02 ...  6.400e-01  6.400e-01]
        y: [[ 0.000e+00 -1.530e-06 ... -1.058e+00 -1.058e+00]
            [ 0.000e+00 -1.593e+07 ... -1.891e+10 -1.891e+10]]
      sol: None
 t_events: None
 y_events: None
     nfev: 6727
     njev: 0
      nlu: 0
c1=  2.038723587753342