实施Radau-IIA方法

Question

我正在寻找一个Radau-IIA方法的实现，在MATLAB中使用3阶来解决DAE系统。这是我的系统（Rober问题）：

y1'=-0.04*y1+10^4*y2*y3
y2'= 0.04*y1-10^4*y2*y3-3*10^7*y2^2
0=y1+y2+y3-1

初始值为

y1(0)=1 and y2(0)=y3(0)=0.

Runge-Kutta画面看起来像这样

1/3  |  5/12  -1/12
1    |  3/4    1/4
_____|_______________
     |  3/4    1/4

（抱歉，我不知道是否有可能在StackOverflow中编写LaTeX或其他内容。）

我不知道该怎么做..