Bogacki Shampine的定制Simulink离散时间积分器模块

Question

我正在尝试使用Bogacki Shampine规则在Simulink中创建自己的离散时间积分器。规则的一般公式（当它只是时间的函数时）为：

y(n+1) = y(n) + (t/9)*(2*s1+3*s2+4s3) 

其中：

s1 = x(n)
s2 = x(n+h/2)
s3 = x(n+3h/4)

也等于：

y(n) = y(n-1) + (t/9)*(2*s1+3*s2+4s3) ; 

其中：

s1 = x(n-1)
s2 = x(n-h/2)
s3 = x(n-h/4)

然后，我将结果与使用ode3的简单积分器块（Bogacki Shampine）进行了比较。结果彼此接近，但并不过分。 另外，我不确定我是否以正确的方式创建了此集成器。由于Bogacki Shampine是3阶。我以为我应该使用3个单位延迟，但是2个就足够了。

如何改善这一点或创建另一个以获得更准确的结果？