23.08.2018 Views

algebralineal2

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

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

1.5 Cálculo con Python 25<br />

esto es,<br />

arctan √ 2<br />

2 = 0.615479708670387<br />

que al ser un método y no una función, no precisa importación.<br />

1 5 2 Operaciones avanzadas<br />

El módulo SymPy permite realizar cálculos más sofisticados con los números<br />

complejos: por ejemplo, reconoce la fórmula de Euler, y es capaz de simplificar<br />

operaciones trigonométricas:<br />

1 >>> from sympy import Symbol ,E,I,re<br />

2 >>> x= Symbol (’x’, real = True )<br />

3 >>> a=E **( I*x)<br />

4 >>> b=a. expand ( complex = True )<br />

5 >>> b<br />

6 I* sin (x) + cos (x)<br />

7 >>> c=a **3<br />

8 >>> d=c. expand ( complex = True ). expand ( trig = True )<br />

9 >>> d<br />

10 -3* cos (x) - I* sin (x) + 4*I* cos (x) **2* sin (x)<br />

11 + 4* cos (x) **3<br />

12 >>> re(d)<br />

13 -3* cos (x) + 4* cos (x) **3<br />

14 >>> f=c. expand ( complex = True )<br />

15 >>> re(f)<br />

16 cos (3* x)<br />

17 >>> (re(f)-re(d)). expand ( trig = True )<br />

18 0<br />

Nótese el uso de la etiqueta adicional real=True en la función Symbol, cuyo<br />

significado es evidente: asume que la variable simbólica x es real. El número e<br />

viene dado por E y la parte real es re. El método expand sirve para desarrollar<br />

expresiones, y admite etiquetas adicionales para usar desarrollos con complejos<br />

o trigonométricos (admite otras varias como logaritmos, potencias, etc.).<br />

Los cálculos prueban la fórmula:<br />

cos(3x) = −3 cos(x) + 4 cos 3 (x)<br />

¿Puede el lector encontrar una fórmula análoga para sen(3x)?<br />

El módulo SymPy también permite encontrar las raíces complejas de polinomios.<br />

La función solve nos proporciona las raíces de una ecuación del siguiente<br />

modo:<br />

1 >>> from sympy import Symbol , solve<br />

2 >>> x= Symbol (’x ’)<br />

3 >>> solve (x**4 -1 ,x)<br />

4 [1 , -1, -I, I]

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

Saved successfully!

Ooh no, something went wrong!