dae2.m solves a set of differential algebraic equations (DAEs)
f(t,y,y')=0 where y'=dy/dt
with a 2nd order method starting from y0 at time t0 and finishing at time tfin where tspan=[t0 t1 ... tfin].
The method will also work well for stiff sets of ODEs.
See pendrun.m, penddae.m & pendg.m for a pendulum example.
See dae4.m and dae4o.m for higher order accurate versions.
I need in my research the DAE model to apply in my power system test
thank`s for you
Let p=dy/dx. Consider a single equation f(x,y,p)=0 (i.e. y=scalar) . Its general algebraic form (with respect to p) is
c(x,y)= row-vector of coefficients (1 by N>2, N being polynomial order), possibly functions of x,y.
P=p.^(N-1:-1:0)' = column of derivative powers
The equation c*P=0 solved for p has N-1 (more than one) roots.
The question: which of these roots will be selected and what are selection criteria?
I need this file
Please send for me
Thanks for help me
But if you can tell us the theory of the solver?
How does this compare to MATLAB's stiff ODE solvers?
Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.