20.06.2013 Views

Acceso al documento en PDF - Biblioteca Nacional de Maestros

Acceso al documento en PDF - Biblioteca Nacional de Maestros

Acceso al documento en PDF - Biblioteca Nacional de Maestros

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.

Com<strong>de</strong>x/Infocom Arg<strong>en</strong>tina'97 Roberto E. C<strong>al</strong>igaris, Graciela A. Mansilla y MartaC. C<strong>al</strong>igaris<br />

En t<strong>al</strong>es casos se <strong>de</strong>be solicitar <strong>al</strong> programa que resuelva la ecuación difer<strong>en</strong>ci<strong>al</strong> <strong>en</strong><br />

forma numérica y para ello basta con escribir una N <strong>de</strong>lante <strong>de</strong> DSolve y un<br />

pequeño cambió <strong>en</strong> las instrucciones como se indica a continuación:<br />

NDSolve[{i V 1<br />

[t]+(R/L)i l<br />

[t]+l/(L c3 i[t]==<br />

fern[t]/L,i[O]-0,i~[01==03,i,~t,0,4.31<br />

y Mathematica nos indica que ha efectuado las cu<strong>en</strong>tas <strong>en</strong> forma numérica con el<br />

m<strong>en</strong>saje:<br />

{{i -> Inter~olatingFunction[{O., 4.}, O]}}.<br />

El interv<strong>al</strong>o <strong>de</strong> tiempo pedido <strong>al</strong> programa fue (0,4) y para cada caso es el<br />

resultado o <strong>de</strong> un análisis cu<strong>al</strong>itativo <strong>de</strong>l problema, o, lo que es más común, por el<br />

inf<strong>al</strong>ible sistema <strong>de</strong> prueba y error. Si el lector quiere dibujar la respuesta <strong>en</strong>tre los<br />

v<strong>al</strong>ores <strong>de</strong> t 0 seg y 4 seg es libre <strong>de</strong> hacerlo. Se lo advertimos. Nosotros<br />

mostramos el sigui<strong>en</strong>te dibujo:<br />

Plot[Ev<strong>al</strong>uate[i[t]/.%],{t,0,0.002),PlOtStYle-><br />

{Text["i (t) ll,, {0.001,6 lo"{-633131;<br />

- Eì<br />

7.10<br />

-6<br />

6.10<br />

-6<br />

5.10<br />

i(t)<br />

-6J<br />

4.10<br />

I / -<br />

O.OQO5 0.001 0.0015 0.002<br />

Observe el lector los v<strong>al</strong>ores <strong>de</strong> i y <strong>de</strong> t.<br />

El otro ejemplo discutido sigui<strong>en</strong>do con los v<strong>al</strong>ores ya conocidos es:<br />

R-0.560; L=lO; c= 0.1;<br />

los comandos que se utilizan y la curva que correspon<strong>de</strong> a la int<strong>en</strong>sidad <strong>de</strong><br />

corri<strong>en</strong>te <strong>en</strong> el circuito es:<br />

NDSolve[(i "[t]+(R/L)i l<br />

[t]+l/(L c) i[t]==<br />

fern[t]/L,i[0]==0,i 1<br />

[0]==03,i,{t,0,5031<br />

81

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

Saved successfully!

Ooh no, something went wrong!