Как получить доступ к определенной точке данных из связанных решений ODE, используя условные

Question

Я решаю 2 связанных ODE, dP / dt и dM / dt в R, используя код ниже. Теперь я хочу использовать условное, так что для некоторого заданного значения P условное сохраняет соответствующее значение t и M в 2 векторах 'x' и 'y'. Как я могу получить доступ к этим конкретным значениям T и M из "Out".

state <- c(P=0.01,M=0) #initial cond

save_1 <- function(t,state,parmeters){

with(as.list(c(state,parameters)),{

  dP <- (-4)*((P)-(4*(mu)^3/(27)))*(M+P*t^3)/(t^2-2*a*t*M)

  dM <- (c*t^2)*(3*P+4*b-(16*(mu)^3*l/(27*e)))


  list(c(dP,dM))
})

} parameters<-NULL

# callingODE

out <- ode(y=state,func = save_1,times=times,parms=NULL)

}