四阶龙格库塔方法（ODE45）自变量的步进如何设置

MILAN3

示例程序如下：

f=@(t,y)y+t^2;                   %通过函数名柄定义函数 常微分方程 dy/dt=y+t^2

tspan=[1 4];                        %自变量T的取值范围

y0=-2;                                 %初始值，即T=1时，Y=-2

[t,y]=ode45(f,tspan,y0);      %采用四阶龙格库塔方法求解微分方程

plot(t,y) %作图

计算后发现，程序将T的取值范围分了40份，即T的步进为3/40＝0.075，算法进行了40次迭代

请教高手，如何设定T的步进，如设置T的步进为3/50=0.06，让算法进行50次迭代

inndoor

http://www.mathworks.cn/support/solutions/en/data/1-1TJ3GZ/

定步长积分的，两个欧拉三个RK。