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