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

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

36 CAPITOLUL 1<br />

> dsolve(diff(x(t),t)=(t/(1+t))*(1-x(t)));<br />

<br />

et x (t) =<br />

1+t<br />

+ C1<br />

<br />

(e −t + e −t t).<br />

Dacǎ dorim ca solut¸iile sǎ fie afi¸sate sub formǎ parametricǎ, atunci se<br />

folose¸ste argumentul opt¸ional ‘parametric‘ ¸si obt¸inem:<br />

> dsolve(diff(x(t),t)=(t/(1+t))*(1-x(t)),x(t),‘parametric‘);<br />

x (t) = 1 − e−t<br />

C1 − e−t t<br />

C1 .<br />

Se mai poate utiliza ca argument opt¸ional ”metoda <strong>de</strong> rezolvare a<br />

ecuat¸iei”. Dacǎ dorim sǎ se rezolve ecuat¸ia diferent¸ialǎ ca o ecuat¸ie<br />

liniarǎ, atunci se folose¸ste argumentul opt¸ional [linear] ¸si obt¸inem:<br />

> dsolve(diff(x(t),t)=(t/(1+t))*(1-x(t)),x(t),[linear]);<br />

<br />

et x (t) =<br />

1+t<br />

+ C1<br />

<br />

(e −t + e −t t),<br />

iar dacǎ dorim sǎ se rezolve ecuat¸ia diferent¸ialǎ ca fiind o ecuat¸ie<br />

cu variabile separate, atunci folosim argumentul opt¸ional [separable] ¸si<br />

obt¸inem:<br />

> dsolve(diff(x(t),t)=(t/(1+t))*(1-x(t)),x(t),[separable]);<br />

x (t) = ( C1 et−1−t)e−t . C1<br />

Nespecificând metoda <strong>de</strong> rezolvare Maple va alege una dintre ele.<br />

Deoarece în secvent¸ele <strong>de</strong> mai sus nu s-a dat nici o condit¸ie init¸ialǎ, Maple a<br />

afi¸sat rǎspunsul cu ajutorul unei constante necunoscute. Dacǎ specificǎm ¸si<br />

condit¸ia init¸ialǎ atunci calculatorul va rezolva o problemǎ cu condit¸ii init¸iale<br />

(Problemǎ Cauchy) ¸si va afi¸sa solut¸ia acesteia.<br />

Pentru ecuat¸ia diferent¸ialǎ (1.62) vom consi<strong>de</strong>ra douǎ Probleme Cauchy<br />

<strong>de</strong>oarece domeniul <strong>de</strong> <strong>de</strong>finit¸ie al membrului drept este reuniunea<br />

(−∞, −1) × IR 1 ∪ (−1, +∞) × IR 1 .<br />

Dacǎ consi<strong>de</strong>rǎm t > −1 ¸si condit¸ia init¸ialǎ x(2) = 4, atunci se obt¸ine<br />

solut¸ia:<br />

> dsolve({diff(x(t),t)=(t/(1+t))*(1-x(t)),x(2)=4},x(t));<br />

x (t) =<br />

<br />

et 1+t − 1/3 e−2e2−4 e−2 <br />

(e −t + e −t t),<br />

iar dacǎ consi<strong>de</strong>rǎm t < −1 ¸si condit¸ia init¸ialǎ x(−2) = 0, atunci se<br />

obt¸ine solut¸ia:<br />

> dsolve({diff(x(t),t)=(t/(1+t))*(1-x(t)),x(-2)=0},x(t));

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

Saved successfully!

Ooh no, something went wrong!