06.03.2013 Views

Capitolul 1 Ecuatii diferentiale de ordinul ˆıntâi rezolvabile prin ...

Capitolul 1 Ecuatii diferentiale de ordinul ˆıntâi rezolvabile prin ...

Capitolul 1 Ecuatii diferentiale de ordinul ˆıntâi rezolvabile prin ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Calculul simbolic al solut¸iilor ecuat¸iilor <strong>de</strong> <strong>ordinul</strong> n 69<br />

<strong>de</strong>qns - ecuat¸ia diferent¸ialǎ <strong>de</strong> orice ordin pe care dorim sǎ o<br />

rezolvǎm sau lista <strong>de</strong> ecuat¸ii diferent¸iale <strong>de</strong> <strong>ordinul</strong> întâi<br />

(în cazul sistemelor)<br />

vars - variabila in<strong>de</strong>pen<strong>de</strong>ntǎ sau lista variabilelor in<strong>de</strong>pen<strong>de</strong>nte<br />

trange - domeniul <strong>de</strong> <strong>de</strong>finit¸ie al variabilei in<strong>de</strong>pen<strong>de</strong>nte<br />

inits - lista <strong>de</strong> condit¸ii init¸iale<br />

xrange - domeniul <strong>de</strong> variat¸ie al primei variabile <strong>de</strong>pen<strong>de</strong>nte<br />

yrange - domeniul <strong>de</strong> variat¸ie al celei <strong>de</strong>-a doua variabile<br />

<strong>de</strong>pen<strong>de</strong>nte<br />

options - diferite opt¸iuni: modul <strong>de</strong> afi¸sare al solut¸iei, metoda <strong>de</strong><br />

rezolvare, etc.<br />

> with(DEtools):DEplot(eq3,x(t),t=-4..1,[[x(0)=0,D(x)(0)=2,<br />

(D@@2)(x)(0)=4]],x=-0.6..1.8,stepsize=.05,title=‘Solutia<br />

Problemei Cauchy‘);<br />

Figura 11<br />

4. Ecuat¸ia diferent¸ialǎ liniarǎ <strong>de</strong> <strong>ordinul</strong> al treilea cu coeficient¸i variabili<br />

<strong>de</strong> tip Euler:<br />

t 2...<br />

x + 5t¨x + 4˙x = ln t, t > 0 (2.33)<br />

> eq4:=t^2*diff(x(t),t,t,t)+5*t*diff(x(t),t,t)+<br />

4*diff(x(t),t)=ln(t);<br />

eq4 := t 2 d3<br />

dt 3x (t) + 5 t d2<br />

dt 2x(t) + 4 d<br />

dt x (t) = ln (t)<br />

> dsolve({eq4,x(2)=2,D(x)(2)=1/2,(D@@2)(x)(2)=3});

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!