中心差分法

问：

我编写了一段代码，使用中心差分法在MATLAB中求解常微分方程。 脚本文件是

clc;
clear;
close all;
% Step size and paramteres
h = 0.0001;
tspan = 0:h:10;
initial_condition = 3;

% ode45 solver
[t,y] = ode45(@central_diff_f1,tspan,initial_condition);
plot(t,y);

% using the central difference
[t1,y1] = central_diff_f2(tspan,initial_condition,h);
hold on
plot(t1,y1)
legend('ODE 45','Central-Difference');

在这里，我尝试重叠解决方案，我的函数文件如下，

function [t1,y1] = central_diff_f2(tspan,initial_condition,h)

y1 = zeros(size(tspan));
y1(1) = initial_condition;
t1 = tspan;
% implementing central difference approximation
for i=2:1:length(tspan)
    y1(i) = y1(i-1) + h*(3*y1(i-1)-6);
end

现在我的问题是，

我们知道，在数学上使用中心差分，

如果我们重新排列，我们得到
y（t+h） ≈ y（t-h）+2h⋅y′（t） 我们代入方程 y' = 3y - 6，我们得到
y（t+h） ≈ y（t-h）+2h⋅（3y-6）

现在比较这个用代码编写的，
y1(i) = y1(i-1) + h*(3*y1(i-1)-6);

为什么我们把 t-h 作为 i-1，为什么 h 的系数取 1？（h的系数必须为2）