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: 

 

non  eta  edozein funtzio diferentziagarri diren (ez dute linealak izan beharrik), eta 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 

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 bat Newtonen bigarren mugimendu-legea da. Horren arabera, F indarraren pean dagoen gorputzaren desplazamenduaren eta t denboraren arteko erlazioa ekuazio diferentzial baten bidez adierazten da:

 

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  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:

 

Oro har, n ordenako ekuazio diferentzial arrunt inplizitu bat honelakoa da:

 

Beste sailkapen batzuk ere badaude:

 Autonomo
 x aldagaiarekiko dependentziarik ez duen ekuazio diferentziala
 Lineal
 Ekuazio diferentziala lineala dela esaten da baldin eta F idatz badaiteke y-ren deribatuen konbinazio lineal gisa: 
 
non a i (x) eta r (x) funtzioak x aldagaiarekiko jarraituak diren.
 r(x) funtzioaren arabera ondoko bi klasifikazio garrantzitsu egin daitezke:

 Homogeneo
 r(x) = 0 denean, eta ondorioz, soluzio "automatikoa" daukanean y = 0 soluzio tribiala. 
 Ez-homogeneo
 r(x) ≠ 0 denean. 

 Ez-lineal
 Konbinaziolineal baten moduan idatzi ezin den  ekuazio diferentziala.

 EDAen sistema 
Elkar nahasian dauden hainbat ekuazio diferentzialek ekuazio sistema bat osatzen du. Baldin eta  y bektore bat bada eta bere elementuak funtzioak badira, hau da,  y bektorearen osagai bakoitza funtzio bat bada,: y(x) = [y1(x), y2(x),..., ym(x)], eta  F funtzio bektorial bat bada, hau da,  F bektoreko elementu bakoitza y-ren eta haren deribatuen funtzio bat bada, orduan, n ordenako eta m dimentsioko ekuazio diferentzial arrunten sistema esplizitua da:

 

Zutabe-bektore gisa:

 

Funtzio bektorialeko  f1(), f2(),..., fm() funtzioek ez dute nahitaez linealak izan beharrik. Horren pareko sistema inplizitua da:

 

non  0= (0, 0, ..., 0) zero bektorea den. Matrize gisa

 

 itxurako sistema baterako, iturri batzuek  eskatzen dute  EDA [sistema] inplizitua izendatzeko  matrize jacobiarra  ez-singularra izatea: ez-singularra izateko baldintza betetzen duen EDA inplizituen sistema transforma daiteke EDA esplizituen sistema batean. Iturri berberetan, jacobiar singularra duten  EDA inplizituen sistemei ekuazio aljebraiko diferentzial (EAD) deritze. Bereizketa hori ez da terminologia hutsa; EADek funtsean ezaugarri desberdinak dituzte, eta, oro har, EADen sistemak ebaztea zailagoa da EDAen sistemak (ez-singularrak) ebaztea baino. Irizpide honi jarraituz, ordena altuagoko ekuazio sistemetan, hurrengo deribatuak, matrize Hessiarra eta hurrengoak,  ez-singularrak izatea eska daiteke, Dena den, 1 ordena baino altuagoko edozein EDA alda daiteke eta berridatz daiteke (gainera, ohikoena horixe egitea da) lehen ordenako EDAen sistema gisa. Ondorioz, jacobiarraren singularitatearen irizpidea nahikoa da ordena edozein dela ere.

EDO sistema baten portaera fase-aldaketen bidez ikus daiteke.

Soluzioak 
Ekuazio diferentzial bat emanda

 

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 

 

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

 

Luzapenik ez duen soluzio bati soluzio maximal deritzo.  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 . 

EDA linealaren testuinguruan, soluzio partikular terminoa EDAren edozein soluziori ere badagokio (hasierako baldintzak bete ez arren), zeina soluzio homogeneoari gehitzen zaion (EDA homogeneoaren soluzio orokor bat); Eta honek, jatorrizko EDAren soluzio orokor bat eratzen du. Terminologia hori erabiltzen da artikulu honetako asmatze-metodoaren atalean, eta maiz erabiltzen da koefiziente zehaztugabeen metodoa eta parametroen aldaketa eztabaidatzen direnean.

Lehenengo ordenako EDA baten soluzioa

Har dezagun lehen ordenako ekuazio diferentzial arrunta:

Bere soluzio orokorrak K konstante arbitrario baten menpekotasuna izango du, eta soluzio hori honako itxurako soluzioa izango da:

K konstantearen edozein baliorentzat ekuazio diferentzialaren baldintza bete egingo du soluzioak. Gainera, ekuazio diferentzialak hasierako baldintzaren bat bete behar badu, esate baterako 

baldintza bada, orduan beti aurkitu ahal izango dugu K konstante arbitrarioarentzat  K0 balioa zeinarentzat   funtzioak, ekuazio diferentziala betetzeaz gain, hasierako baldintza ere bete egingo duen. Beti ere suposatzen da  puntua soluzioaren esistentzia eta bakartasuna betetzen den tarteren batean dagoela.

Iraupen mugatuko soluzioak 
EDA autonomo ez-linealen kasuan, posible da, zenbait baldintzatan, iraupen finituko soluzioak garatzea.  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:

 

Iraupen mugatuko soluzioa onartzen du:

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. mendeko 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 diferentzialen 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:Izan bedi, era berean R laukizuzena x-y planoan definitua: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 eta tarte horretan hasierako balioaren problemak soluzioa dauka. Hau da, soluzioa dauka eta gainera bakarra da. F lineala izateko murrizketarik 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, tarte ireki maximo bakarra dago (ziur aski infinitua):

 

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

Baldin eta , orduan bi aukera daude:

 leherketa denbora mugatuan: 
 Definizio-eremutik irtetea:  

non  den F definituta dagoen tarte irekia, eta  bere muga.

Kontuan izan soluzioaren gehienezko domeinua

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

 Adibidez,

 

Horrek esan nahi du  dela eta hori  denez, lokalki Lipschitz jarraitua da, beraz, Picard-Lindelöf teorema betetzen du.

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

 

eta horren domeinu maximoa da:

 

Horrek argi erakusten du tarte maximoa hasierako baldintzen araberakoa izan daitekeela. y-ren domeinutzat har liteke  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 , izan ere

 

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,

 

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

 

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

 

bektore-notazioarekin era trinkoagoan adieraz daiteke:

 

non

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  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,  eta   integrazioko alegiazko aldagaiak dira (batutze-indizeen analogo jarraituak), eta  notazioak adierazten du  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 , 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:

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)

Ekuazio diferentzialak
Kalkulu diferentziala