R中使用ode23的非数字参数错误

Question

为了理解和学习ode23s方法在R中的工作原理，我使用了一本来自数值分析书的简单系统并实现了该系统。我遇到错误

t + y [1] ^ 2中的错误：二进制运算符的非数字参数

这是我的代码

#x'=t+x^2-y
#y'=t^2-x+y^2
#x(0)=3, y(0)=2

model  <- function(t,y) as.matrix(c(t+y[1]^2-y[2], t^2-y[1]+y[2]^2))
sol <- ode23s(model, 0, 20, c(3,2), jac = NULL, hmax = 0)

在我看来t不是数字，并且我也用作.numeric（t）但也没有用。您能解释一下如何解决此错误吗？