Ekuazio diferentzial arrunt

Matematikan, **ekuazio diferentzial arrunta** (EDA) ekuazio diferentzial bat da. Ekuazio diferentzial horretan ezezaguna aldagai independente bakarreko funtzio bat edota aldagai independente bakarreko hainbat funtzio dira, eta ekuazioan bertan ageri da bai funtzio ezezaguna[k] eta bai funtzio ezezagun horren edo horien deribatuak.  *Arrunta* terminoa erabiltzen da ekuazio diferentzial partzial terminoarengandik bereizteko, izan ere,  ekuazio diferentzial partzialetan funtzio ezezagunean hainbat aldagai independente ager daitezke, hau da, aldagai independente bat baino gehiago izan ditzake eta aldagai horietako bakoitzarekiko deribatu partzialez hitz egin beharko genuke.

## Ekuazio diferentzialak 
Ekuazio diferentzial arrunten artean, ekuazio diferentzial linealek garrantzi nabarmena dute zenbait arrazoirengatik. Fisikan eta matematika aplikatuan dauden funtzio gehienak, bai  oinarrizkoak eta bai bereziak,  ekuazio diferentzial linealen ebazpenak dira (ikus funtzio holonomikoa). Fenomeno fisikoak ekuazio ez-linealen bidez modelizatzen direnean, gehienetan ekuazio diferentzial linealen bidez hurbiltzen dira, ebazpena errazago lortzeko. Esplizituki ebatz daitezkeen EDA ez-lineal gutxi horiek ebazteko, ekuazioa EDA lineal baliokide bihurtzen da normalean (ikus, adibidez, Riccati ekuazioa). 

Ekuazio diferentzial lineala ekuazio diferentzial bat da, funtzio ezezagunean eta haren deribatuetan polinomio lineal batek definitzen duena, hau da, honako itxura duen ekuazioa: 

: a_0(x)y +a_1(x)y' + a_2(x)y'' +\cdots +a_n(x)y^{(n)}+r(x)=0,

non a_0(x), a_1(x),\cdots +a_n(x) eta r(x) edozein funtzio diferentziagarri diren (ez dute linealak izan beharrik), eta y, y', +\cdots y^{(n)}funtzio ezezaguna eta funtzio ezezagunaren '''' aldagaiarekiko deribatuak diren.

Hainbat EDA esplizituki ebatz daitezke funtzio eta integral ezagunen bidez. Hori ezinezkoa denean, soluzioaren Taylor seriea konputatzeko ekuazioa erabilgarria izan daiteke. Problema aplikatuetarako, ekuazio diferentzial arruntetarako zenbakizko metodoek soluzioaren hurbilketa bat eman dezakete, esate baterako Runge-Kutta metodoek.

## Testuingurua 
[alt=parabolic projectile motion showing velocity vector|eskuinera|thumb|250x250px|[Kanoi batetik jaurtitako jaurtigai baten ibilbideak Newtonen bigarren legetik ondorioztatzen den ekuazio diferentzial arrunt batek zehaztutako kurba bati jarraitzen dio.]]
Ekuazio diferentzial arruntak (EDA) matematikaren eta gizarte- eta natura-zientzien testuinguru askotan sortzen dira. Aldaketen adierazpide matematikoa deribatuen eta diferentzialen bidez egiten da eta ekuazio batek erlaziona ditzake hainbat diferentzial, deribatu eta funtzio. Horrela, ekuazio diferentzial batek deskribatzen du dinamikoki aldatzen den fenomeno bat, bere bilakaera edota aldaketa. Askotan, kantitateak definitzen dira beste kantitate baten aldaketa-tasa moduan  (adibidez, desplazamenduaren deribatuak denborarekiko) edo kantitateen gradiente moduan, eta horrela ekuazio diferentzialak sortzen dira.

Ekuazio diferentzialak agertzen diren alorrak dira, matematika-eremuetan, geometria eta mekanika analitikoa, eta zientzia-alor ezberdinetan aipa genitzake, besteak beste, fisika orokorrean,  astronomia (zeruko mekanika), meteorologia (modelizazio meteorologikoa), kimika (erreakzio-tasak), biologia (gaixotasun infekziosoak, aldakuntza genetikoa), ekologia eta populazioaren modelizazioa (populazioaren lehia), ekonomia (erreserben joerak, interes-tasak eta merkatuko prezioen aldaketak).

Matematikari askok ekuazio diferentzialak aztertu izan dituzte eta gaiari aurrerapen nabarmenak eskaini dizkiote, aipatzeko modukoak dira, Newton, Leibniz, Bernoulli familia, Riccati, Clairaut, d'Alembert eta Euler.

Adibide sinple ba*t* Newtonen bigarren mugimendu-legea da. Horren arabera, F indarraren pean dagoen gorputzaren desplazamenduaren eta t denboraren arteko erlazioa ekuazio diferentzial baten bidez adierazten da:

: m \frac{\mathrm{d}^2 x(t)}{\mathrm{d}t^2} = F(x(t))\,

*Ekuazioak m masa konstanteko partikula baten mugimendua mugatzen du.* Oro har, denboraren t unean partikularen x(t) posizioaren funtzio bat da F. Funtzio ezezaguna, hau da, x(t) funtzioa,  ekuazio diferentzialaren bi aldeetan ageri da, eta F(x(t)) notazioan adierazten da. 
 

## Definizioak 
Aurrerantzean, y mendeko aldagai bat izango da, eta *x* aldagai independentea, horrela, y = f(x) izango da x aldagaiaren mende dagoen funtzio ezezaguna. Diferentziazioaren notazioa aldatu egiten da autoreen *eta* lan eremuen arabera. Testuinguru horretan, Leibnizen notazioa  erabilgarriagoa da diferentziazioan eta integrazioan baina Lagrangeren notazioa  erabilgarriagoa da edozein ordenatako deribatuak adierazteko. Bestalde, Newtonen notazioa (\dot y, \ddot y, \overset{...}{y}) fisikan asko erabiltzen da denborarekiko ordena txikiko deribatuak adierazteko.

### Definizio orokorra 
Izan bedi F funtzioa x-ren, y-ren *eta* y-ren deribatuen funtzio bat. Horrela,  '''*n* ordenako ekuazio diferentzial arrunt'* *esplizitu'' deritzo  honako ekuazioari:

: F\left (x,y,y',\ldots, y^{(n-1)} \right )=y^{(n)}

Oro har, *n* ordenako ekuazio diferentzial arrunt *inplizitu* bat honelakoa da:. Ondorioz, jacobiarraren singularitatearen irizpidea nahikoa da ordena edozein dela ere.

EDO sistema baten portaera fase-aldaketen bidez ikus daiteke.

### Soluzioak 
Ekuazio diferentzial bat emanda

: F\left(x, y, y', \ldots, y^{(n)} \right) = 0

u funtzioa,  *non I tarte bat den,*  F-ren soluzio edo kurba integral bat izango da baldin eta u n aldiz diferentziagarria bada I tartean eta 

: F(x,u,u',\ \ldots,\ u^{(n)})=0 \quad x \in I.

Bi soluzio emanda,   eta , u izango da v-ren *luzapena* baldin eta  eta

: u(x) = v(x) \quad x \in J.\,

Luzapenik ez duen soluzio bati *soluzio maximal* deritzo. \mathbf{R} osoan definitutako soluzio bati *soluzio global* deritzo.

n ordenako ekuazio baten soluzio orokor batek  n integrazio-konstante independente arbitrario edukiko ditu.  *Soluzio* *partikular* bat soluzio orokorretik eratortzen da, integrazio-konstanteei balio partikularrak ezarriz;  balio partikular hauek aukeratzen dira ekuazioari ezarritako  "hasierako baldintzak edo mugako baldintzak" bete daitezen. Soluzioa singularra izango da baldin eta soluzioa ezin bada lortu  soluzio orokorreko konstante arbitrarioei balio definituak ezarriz .  Horrek esan nahi du, beren dinamikatik abiatuta, sistemak zero balioa lortuko duela azken unean, eta zeroan egongo dela gero betirako. Iraupen finituko soluzio horiek ezin dira funtzio analitikoak izan lerro erreal osoan, eta, amaierako denboran Lipschitz-funtzioak ez direnez, ez dute beteko soluzio bakarra izatearen propietatea,  Lipschitzen ekuazio diferentzialen soluzioek duten berezitasuna.

Adibidez, ekuazio honek:

: y'= -\text{sgn}(y)\sqrt,\,\,y(0)=1

Iraupen mugatuko soluzioa onartzen du:

: y(x)=\frac{1}{4}\left(1-\frac{x}{2}+\left|1-\frac{x}{2}\right|\right)^2
:

### Hasierako baldintzak 
Orokorrean, ekuazio diferentzialak bete beharreko hasierako edo mugako baldintzarik zehazten ez bada, ekuazioak ez du soluzio [partikular] bakarra izango. Ekuazio diferentzial lineal baten ordena n bada, hasierako edo mugako n baldintza behar dira ekuazio diferentziala eta hasierako edo mugako  baldintza guztiak betetzen dituen soluzio bakarra egon dadin. 

### Ebazpen metodoak 
Mota askotako Ekuazio Diferentzial Arruntak daude, mota bakoitzak ebazpena lortzeko bide ezberdina du. Ordenaren arabera bereiztea komeni da, izan ere lehen ordenakoak ebaztea errazagoa izan ohi da.

#### Soluzio analitikoak 
Badaude hainbat metodo Ekuazio Diferentzial Arrunt linealen soluzio analitikoak lortzea ahalbidetzen dutenak. Bereziki, ekuazio linealeko koefizienteak konstanteak edota periodikoak direnean, soluzioa eratzea ez da zaila izaten. Koefizienteak ez direnean konstanteak edo periodikoak, baina Taylor Serieen edo Laurenten serieen bidez adieraz daitezkeenean, hainbat murrizpen betetzen direnean Frobeniusen metodoa aplika daiteke. Beste aukera bat izan daiteke n ordenako ekuazio diferentzial arrunta n ekuazio diferentzial lineal dituen sistema batean bihurtzea.

Ezkuazio diferentzial arrunta ez bada lineala ez dago metodo orokorrik.

#### Zenbakizko  metodoak 
Besterik ezin denean, zenbakizko metodoak erabiltzen dira Ekuazio Diferentzial Arrunten soluzioen hurbilketak lortzeko. Metodo hauek ere asko dira eta problemaren arabera metodo egokiena aukeratzea komeni da.

## Teoriak 
### Soluzio singularrak 
Ekuazio diferentzial arrunten eta partzialen soluzio singularren teoria Leibnizen garaitik ikertu zen, baina XIX. mendearen erdialdetik aurrera bakarrik jaso du arreta berezia. Gaiari buruzko lan baliotsua, baina ez oso ezaguna, Houtainena da (1854). Darboux (1873) teoriaren buru izan zen, eta soluzio horien interpretazio geometrikoaren inguruko ikerketa lerroa ireki zuen. Hainbatek ekin zion bide horri, aipagarrienak Casorati eta Cayley. Azken horri zor zaio (1872) lehen ordenako ekuazio diferentzialen soluzio singularren teoria, 1900 inguruan onartutakoa.

### Koadraturetara murriztea 
Ekuazio diferentzialekin jarduteko lehen saiakerek ekuazio diferentziala koadraturetara murriztean oinarritu ziren. XVIII. me*n*deko algebraistek n graduko ekuazio orokorra ebazteko metodo bat aurkitzeko itxaropena izan zuten bezala, analisten itxaropena zen edozein ekuazio diferentzial integratzeko metodo orokor bat aurkitzea. Gauss-ek (1799), ordea, ekuazio diferentzial konplexuek zenbaki konplexuak behar dituztela erakutsi zuen. Hala, analistek funtzioen azterketa ordezkatu zuten, eremu berri eta emankorra sortuz. Cauchy izan zen ikuspegi horren garrantzia azpimarratzen lehena.  Hortik aurrera, benetako kontua ez zen ea soluzio bat posible ote zen funtzio ezagunen edo integralen bidez, baizik eta emandako ekuazio diferentzial bat nahikoa ote zen aldagai askearen edo askeen funtzio bat definitzeko, eta, hala bada, zein diren propietate bereizgarriak.

### Fuchsiar teoria 
Fuchsen bi memoriek beste ikuspegi berri bat sustatu zuten,  gerora Thomék eta Frobeniusek landutakoa. Collet-ek asko eragin zuen 1869tik aurrera. Sistema ez-lineal bat integratzeko bere metodoa Bertrand-i jakinarazi zitzaion 1868an. Clebschek (1873) integral abeliarren teorian bezala ekin zion teoriari. Azken hori eraldaketa arrazional baten pean aldatu gabe dagoen oinarrizko kurbaren propietateen arabera sailka daitekeen legez, ekuazio di*f*erentzialen bidez definitutako funtzio transzendenteak ere antzera sailkatzea proposatu zuen Clebsch-ek, hau da, bakoitzari dagokion f = 0 gainazalen propietate inbarianteen arabera, bat-batiko transformazio arrazionalen pean.

### Lieren teoria 
1870etik aurrera, Sophus Lieren lanak oinarri sendoagoak eman zizkion ekuazio diferentzialen teoriari. Frogatu zuen antzinako matematikarien integrazio-teoriek, Lieren taldeak erabiliz, jatorri berekoak direla, eta transformazio infinitesimal berak onartzen dituzten ekuazio diferentzial arruntek integratzeko zailtasun konparagarriak dituztela. Kontaktuen transformazioen gaia ere nabarmendu zuen.

Lieren taldeen teoria ekuazio diferentzialetan ziurtatu egin da: (1) ekuazio diferentzialak ebazteko ezagutzen diren ad hoc metodo asko bateratzen ditu, eta (2) soluzioak aurkitzeko modu berri eta indartsuak ematen ditu. Teoriak ekuazio diferentzial arrunt nahiz partzialetarako aplikazioak ditu. 

Soluzio orokorraren ikuspegi batek ekuazio diferentzialen simetriaren propietatea erabiltzen du, soluzioen transformazio infinitesimal jarraitua soluzioetara (Lie teoria). Talde teoria jarraitua, Lie algebrak eta geometria diferentziala baliagarri dira  ekuazio diferentzial linealen eta ez-linealen (partzialen) egitura ulertzeko, ekuazio integragarriak sortzeko, Lax pareak, errekurtsio-eragileak eta Bäcklund-en transformatua aurkitzeko, eta, azkenik, Ekuazio Diferentzialen soluzio analitiko zehatzak aurkitzeko. 

Simetria-metodoak matematikan, fisikan, ingeniaritzan eta beste hainbat alorretan sortzen diren ekuazio diferentzialei aplikatzen zaizkie.

### Sturm–Liouvilleren teoria 
Sturm–Liouvillen teoria bigarren ordenako ekuazio diferentzial lineal mota berezi baten teoria da. Bere soluzioak autobalioetan eta  dagozkien autofuntzioetan oinarritzen dira, bigarren mailako ekuazio lineal homogeneoen bidez definitutako eragile linealen autofuntzioetan hain zuzen ere. Problemak Sturm-Liouville Problemak (SLP) izenarekin identifikatzen dira  XIX. mendearen erdialedean inguruan J.C.F. Sturmek. eta J. Liouvillek,  aztertu zituztenetik. SLPek autobalio-kopuru infinitua dute, eta dagozkien autofuntzioek multzo osotua eta ortogonala osatzen dute, eta horrek hedaketa ortogonalak egiteko aukera ematen du. Ideia hori funtsezkoa da matematika aplikatuan, fisikan eta ingeniaritzan. SLPak ere erabilgarriak dira ekuazio diferentzial partzial batzuk aztertzeko.

## Soluzioen existentzia eta bakartasuna 
Zenbait teorema daude esaten dutenak hasierako balioaren problemek, tartean EDAek, soluzioa dutela eta soluzioa bakarra dela, bai lokalki eta bai globalki. Hauek dira bi teorema nagusiak:

Oinarrizko moduan, bi teoremek soluzio lokalak baino ez dituzte bermatzen, baina hedatu egin daitezke eta soluzio globalak bermatzeko egokitu, adibidez, Grönwall-en desberdintasunaren baldintzak betetzen badira.

Bestalde, bakartasun-teoremak ez zaizkie aplikatzen EAD sistemei, beren zati aljebraikotik (ez-lineala) abiatuta hainbat soluzio izan baititzakete.

### Existentzia lokalaren eta bakartasunaren teorema sinplifikatuta 
Teorema era sinplifikatuan honela azal daiteke. Izan bedi hasierako balioaren problema bat ekuazio diferentzial baten bidez eta hasierako balio baten bidez adierazia: y' = F(x,y)\,,\quad y_0 = y(x_0)*Izan bedi, era berean R laukizuzena x-y planoan definitua:*R = [x_0-a,x_0+a] \times [y_0-b,y_0+b]non a eta b errealak diren (sinbolikoki:  ),  ikurrak  biderkadura cartesiarra adierazten duen eta kortxeteek tarte itxia adierazten duten.

*F* eta ∂*F*/∂*y jarraituak baldin badira* R laukizuzenean, orduan, existitzen da  balio bat zeinak definitzen duen tarte *bat* I = [x_0-h,x_0+h] \subset [x_0-a,x_0+a]eta tarte horretan hasierako balioaren problemak soluzioa dauka. Hau da, soluzioa dauka eta gainera bakarra da. F lineala izateko murrizk*eta*rik ez dagoenez, F(x, y) forma hartzen duten ekuazio ez-linealei ere aplika dakieke, bai eta ekuazio-sistemei ere.

### Bakartasun globala eta soluzioaren domeinu maximoa 
Picard-Lindelöfen teoremaren hipotesiak betetzen direnean, existentzia lokala eta bakartasuna izaera global batera heda daitezke. Zehatzago esanda:

Hasierako (x0, y0) baldintza bakoitzerako, I_{\text{max}}tarte ireki maximo bakarra dago (ziur aski infinitua):

: I_{\max} = (x_-,x_+), x_\pm \in \R \cup \{\pm \infty\}, x_0 \in I_{\max}

eta hasierako baldintza hori betetzen duen edozein soluzio izango da I_\max domeinuan hasierako baldintzak betetzen dituen soluzioaren  murrizketa (funtzio murriztua).

Baldin eta x_\pm \neq \pm\infty, orduan bi aukera daude:

* leherketa denbora mugatuan:\limsup_{x \to x_\pm} \|y(x)\| \to \infty 
* Definizio-eremutik irtetea: \lim_{x \to x_\pm} y(x)\ \in \partial \bar{\Omega} 

non \Omega den F definituta dagoen tarte irekia, eta \partial \bar\Omega bere muga.

Kontuan izan soluzioaren gehienezko domeinua

* tarte bat da beti (bakartasun baldintza bete dadin)
* \R baino txikiagoa izan daiteke 
* (x0, y0) aukeraren araberakoa izan daiteke.

; Adibidez,

: y' = y^2

Horrek esan nahi du F(x,y) = y^2 dela eta hori C^1 denez, lokalki Lipschitz jarraitua da, beraz, Picard-Lindelöf teorema betetzen du.

Ingurune sinple honetan ere, soluzioaren domeinu maximoak ezin du \R osoa izan, soluzioa hau baita:

: y(x) = \frac{y_0}{(x_0-x)y_0+1}

eta horren domeinu maximoa da:

: \begin{cases}\R & y_0 = 0 \\[4pt] \left (-\infty, x_0+\frac{1}{y_0} \right ) & y_0 > 0 \\[4pt] \left (x_0+\frac{1}{y_0},+\infty \right ) & y_0 

Horrek argi erakusten du tarte maximoa hasierako baldintzen araberakoa izan daitekeela. y-ren domeinutzat har liteke \R \setminus (x_0+ 1/y_0), baina horrek tartea ez den domeinua suposatzen du; ondorioz, hasierako baldintzaren aldea kontrako aldetik deskonektatuta legoke  *eta*, beraz, ezin du berak era bakarrean zehaztu.

Domeinu maximoa ez da \R, izan ere

: \lim_{x \to x_\pm} \|y(x)\| \to \infty,

hain zuzen ere aipatutako teoreman agertutako bi aukera posibleetako bat.

## Ordena murriztea 
Ekuazio diferentzialak errazago ebatz daitezke ekuazioaren ordena murriztu badaiteke.

### Lehen ordenako sisteman bihurtzea 
*n* ordenako edozein ekuazio diferentzial esplizitu,

: \begin{array}{rcl}
  y_1'&=&y_2\\
  y_2'&=&y_3\\
  &\vdots&\\
  y_{n-1}'&=&y_n\\
  y_n'&=&F(x,y_1,\ldots,y_n).
\end{array}

lehe*n* mailako n ekuazio diferentzialen sistema gisa idatz daiteke, *i* = 1, 2,..., *n* balioentzat funtzio ezezagunen familia berri bat definituz

: y_i = y^{(i-1)}.\!

Lehen ordenako ekuazio diferentzial akoplatuen sistema n-dimentsionala hau da:

: \mathbf{y}'=\mathbf{F}(x,\mathbf{y})

bektore-notazioarekin era trinkoagoan adieraz daiteke:

: \mathbf{y}'=\mathbf{F}(x,\mathbf{y})

non

: \mathbf{y}=(y_1,\ldots,y_n),\quad \mathbf{F}(x,y_1,\ldots,y_n)=(y_2,\ldots,y_n,F(x,y_1,\ldots,y_n)).

## Soluzio zehatzen laburpena 
Ekuazio diferentzial batzuek soluzio zehatzak eta mugatuak dituzte. Honako tauletan hainbatt klase garrantzitsu daude.

Hurrengo taulan, , Q(x), P(y), Q(y), eta M(x,y), N(x,y) funtzioak x-ren eta y-ren menpeko edozein funtzio integragarridira, bestalde,   eta c zenbaki erreal konstanteak dira eta C_1, C_2, \ldots edozein balio arbitrario izan ditzaketen konstanteak dira (orokorrean konplexuak izan daitezke). Ekuazio diferentzialak beren forma baliokideetan erakusten dira,  integrazioaren bidez soluziora iristeko moduan.

Integralen soluzioetan, \lambda eta \varepsilon  integrazioko alegiazko aldagaiak dira (batutze-indizeen analogo jarraituak), eta  notazioak adierazten du F(\lambda) funtzioaren integrala egin λ aldagaiarekiko, eta integratu ondoren  λ = x aldaketa egin, konstanterik gehitu gabe.

### Ekuazio bereizgarriak 
### Lehen ordenako ekuazio orokorrak 
### Bigarren ordenako ekuazio orokorrak 
### n. ordenarainoko ekuazio linealak 
## Asmatze-metodoa 
EDA ebazteko beste metodo guztiek huts egiten dutenean, edo Ekuazio Diferentzialaren soluzioa nolakoa izan daitekeen susmoa dugun kasuetan, batzuetan posible izan daiteke ED bat ebaztea soluzioa asmatuz eta ekuazioan ordezkatuz, betetzen ote duen ikusteko.  Metodo hori erabiltzeko, ekuazio diferentzialaren soluzio bat asmatuko dugu, eta, ondoren, soluzioa ekuazio diferentzialean sartuko dugu, ekuazioa betetzen duen balioztatzeko.  Asmatutako soluzioak ekuazioa bete egiten badu, EDaren soluzio partikular bat izango dugu esku artean; betetzen ez badu, berriz hasiko gara eta beste suposizio bat lortzen saiatuko gara.  Adibidez, pentsa genezake ED baten ebazpenak y = Ae^{\alpha t}, izan ere, oso ohikoa den soluzioa da, eta fisikoki portaera sinusoidala du.

Homogeneoa ez den lehen ordenako EDAren kasuan, lehenik eta behin, Ekuazio Diferentzialaren zati homogeneoari (ekuazio bereizgarria ere deitzen zaio) soluzio bat aurkitu behar diogu, eta, ondoren, asmatze metodoaren bidez,  ekuazio ez-homogeneo osoari soluzio bat aurkitu behar diogu. Azkenik, bi soluzioak batzen ditugu EDAaren soluzio osoa lortzeko, hau da:

\text{soluzio osoa} = \text{soluzio homogeneoa} + \text{soluzio partikularra}

## EDAk ebazteko softwarea 
* Maxima, kode irekiko aljebra konputazionaleko sistema bat.
* COPASI, EDAk integratzeko eta aztertzeko software-pakete libreaa (Lizentzia Artistikoa 2.0).
* MATLAB, aplikazio informatiko teknikoa (MATrix LABoratory)
* GNU Octave, goi-mailako lengoaia, batez ere zenbakizko kalkuluetarako erabiltzen dena.
* Scilab, zenbakizko kalkulurako kode irekiko aplikazioa.
* Maple, kalkulu sinbolikoetarako aplikazio patentatua.
* Mathematica, kalkulu sinbolikoetarako aplikazio patentatua.
* SymPy, EDAk sinbolikoki ebatzi dizakeen Python paketea
* Julia (programazio-lengoaia), batez ere zenbakizko kalkuluetarako erabiltzen den goi-mailako lengoaia.
* SageMath, kode irekiko aplikazio bat, Python lengoaiaren  antzeko sintaxia erabiltzen du eta matematikaren zenbait adar hartzen ditu gaitasun sorta zabal batekin.
* SciPy, Python pakete bat, EDA integrazio-modulu bat barneratzen duena.
* Chebfun, MATLABen idatzitako kode irekiko pakete bat.
* GNU R, kode irekiko ingurune konputazionala, batez ere estatistikaren eremuan erabilia baina EDAk ebazteko paketeak dituena.

## Erreferentziak 
## Ikus, gainera 
* Ekuazio diferentzial
* Runge-Kutta metodo

## Kanpo estekak 
*  (bideoa) Ekuazio diferentzialak: koefiziente konstantedun ED lineal homogeneoak eta ez-homogeneoak (1)
*  (bideoa) Ekuazio diferentzialak: koefiziente konstantedun ED lineal ez-homogeneoak (2)
*  (bideoa) Ekuazio diferentzialak: koefiziente konstantedun ED lineal ez-homogeneoak (3)
*  (bideoa) Ekuazio diferentzialak: koefiziente konstantedun ED lineal ez-homogeneoak (4)
*  (bideoa) Ekuazio diferentzialak: koefiziente konstantedun ED lineal ez-homogeneoak (5)