我正在尝试在 Matlab 中解决物理流问题（管道中的边界层流）。为此，我必须解决一个看起来像这样的 ODE 系统：
function dy = ODEFUN(t,y)
dy = dp + c;
end
问题是，该 dp 也是未知的，必须在下一步设置以满足某些物理条件。dp 没有明确的方程，但我只能在下一个时间步 t+1 时检查它是否正确。如果不正确，我需要回到上一个时间步并调整 dp.
似乎像 ODE15s 这样的标准 ODE 求解器无法做到这一点，或者至少我还没有找到一种方法来告诉他们。有没有可能实现这个算法？
感谢您的任何帮助！