File Exchange

image thumbnail

Differential Algebraic Equation Solvers

version 1.0 (7.24 KB) by

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

3.41667
12 Ratings

9 Downloads

Updated

No License

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.

Comments and Ratings (17)

Very nice.

saidi abdelaziz

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

gabriel Hidalgo

alaa taha

thank`s for you

dsfdfd sdd

ss ss

An Jin

Thank you

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?

yang ying

useful programme

jahan jahedpoor

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

¶­(dong) Õ׿¡(zhaojun)

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

hassan hijaizn

phyllis ferguson

puru premi

Brankica Vracar

nidhi sahu

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

MATLAB Release
MATLAB 5.2 (R10)

Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.

» Watch video