No BSD License
by Tony Roberts
08 Sep 1998
A solves a set of differential algebraic equations f(t,y,y')=0 where y'=dy/dt.
Watch this File
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.
This file inspired Jean Le Rand D'alambert Reduction Method (Update:22 06 07).
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?