In [34]:

    

import numpy as np
from sympy import *
from sympy.abc import x as xval


    

In [51]:

    

theta, r, a1, a2, b1, b2, lam1, lam2 = symbols('theta r a1 a2 b1 b2, lam1 lam2')
a3, a4, b3, b4 = symbols('a3 a4 b3 b4')

delta = 0.5354409455


    

In [52]:

    

p1 = lam1 * (a1 * sin(delta * theta) + b1 * cos(delta * theta))


    

In [53]:

    

p1


    

    
Out[53]:





lam1*(a1*sin(0.5354409455*theta) + b1*cos(0.5354409455*theta))

In [54]:

    

p2 = lam2 * (a2 * sin(delta * theta) + b2 * cos(delta * theta))


    

In [55]:

    

j = solve(( - diff(p2, theta)).evalf(subs={theta:0}), a1)[0]


    

    




---------------------------------------------------------------------------
IndexError                                Traceback (most recent call last)
<ipython-input-55-908396efcd82> in <module>()
----> 1 j = solve(( - diff(p2, theta)).evalf(subs={theta:0}), a1)[0]

IndexError: list index out of range

In [ ]:

    

diff(p1, theta)


    

In [ ]:

    




    

In [ ]:

    

xval


    

In [56]:

    

diff(p1, theta)


    

    
Out[56]:





lam1*(0.5354409455*a1*cos(0.5354409455*theta) - 0.5354409455*b1*sin(0.5354409455*theta))

In [59]:

    

p1 = lam1(delta * a1 * cos(delta * theta) - delta * sin(delta * theta))


    

In [ ]: