可以在R中生成具有Holling IV类功能响应的捕食者猎物isocline吗？

Question

使用R，我试图使用捕食者猎物等倾线以图形方式分析具有Holling的IV型功能反应的捕食者 - 猎物相互作用。

这是我正在使用的脚本：

library(deSolve)
library(phaseR)
model.LV <- function(t, y, parameters){
  with(as.list(parameters),{
    N<-y[1] 
    P<-y[2]
    dN <- r*N*(1-N/K)-(1.511*N/((N^2/9.2724)+N+1.483))*P*N
    dP <- -b*P+(1.511*N/((N^2/9.2724)+N+1.483))*a*P*N

    list(c(dN,dP))
  })
}

params.LV<-c(r=0.8,K=10,b=4.2,a=6)

data.LV<-as.data.frame(lsoda(c(N=3,P=3),seq(1,100,by=0.5), model.LV, params.LV))

plot(data.LV$N, data.LV$P, type="l", col="blue", main="Trajectory for Type-IV Functional Response",
     xlab="Prey density N", ylab="Predator density P", xlim=c(0,5), ylim=c(0,5))
nullclines(model.LV, x.lim=c(0,5),y.lim=c(0,5), parameters=params.LV, system="two.dim", colour=c("green","red"), add=TRUE)

等倾线应该有一条抛物线，两条垂直线与它相交，但图形完全错误。这个模拟可以在R中运行吗？我已经在MATLAB中看到过它，并且被告知R可能不足以运行此模拟。

谢谢。