Согласно этому описанию улучшенный метод Эйлера - это просто метод Рунге-Кутты второго порядка ... в любом случае, я не вижу причин, по которым вы не могли бы использовать улучшенный метод Эйлера с тремя уравнениями , Просто примените ту же процедуру к каждой переменной в отдельности, например
x(i+1) = x(i) + 0.5 * dt * (
a*(y(i)-x(i)) +
a*(
y(i)-x(i) +
dt*a*(y(i)-x(i))
)
)
и аналогично для y
и z
.