[y0new, yp0new] = decic(odefun, t0, y0, fixed_y0, yp0, fixed_yp0)
[y0new, yp0new, resnrm] = decic(odefun, t0, y0, fixed_y0, yp0, fixed_yp0, options)
decic adjusts the free components of y0 and yp0 so that the residual odefun(t0,y0,yp0) is small. Components marked by fixed_y0 and fixed_yp0 are kept fixed.
| Item | Details |
|---|---|
| Problem form | Fully implicit residual F(t,y,yp) = 0. |
| Fixed components | fixed_y0 and fixed_yp0 mark values that must not change. |
| Outputs | Consistent initial values y0mod and yp0mod. |
| Used with | ode15i or an ode object with a fully implicit equation type. |
f = @(t,y,yp) yp + y;
[y0, yp0] = decic(f, 0, 1, 1, 0, 0);
[t, y] = ode15i(f, [0 1], y0, yp0)| Version | Description |
|---|---|
| 2.0.0 | initial version |