Differential Algebraic Equation Solvers

A solves a set of differential algebraic equations f(t,y,y')=0 where y'=dy/dt.

Description

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.

• Acknowledgements: This submission has inspired the following:; Jean Le Rand D'Alambert Reduction Method (update:22-06-07)
• MATLAB release: MATLAB 5.2 (R10)

How does this compare to MATLAB's stiff ODE solvers?

[nidhi sahu]

 

[Brankica Vracar]

 

[puru premi]

 

[phyllis ferguson]

 

[hassan hijaizn]

 

[¶­(dong) Õ׿¡(zhaojun)]

Good job!
But if you can tell us the theory of the solver?

[jahan jahedpoor]

I need this file
Please send for me
Thanks for help me

[yang ying]

useful programme

[Vassili Pastushenko]

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

f(x,y,p)=c(x,y)*P

where
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?

[An Jin]

Thank you

[ss ss]

 

[dsfdfd sdd]

 

[alaa taha]

thank`s for you

[gabriel Hidalgo]

 

[saidi abdelaziz]

I need in my research the DAE model to apply in my power system test