I have a function that I have written, attached file, and would like to
pass theta, thetad, phi and phid through it as symbolic variables.
Therefore hoping for an output in terms of these variables that I can then
integrate. Is this possible??
--
You received this message because you are subscribed to the Google Groups
"sympy" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
To post to this group, send email to [email protected].
Visit this group at https://groups.google.com/group/sympy.
To view this discussion on the web visit
https://groups.google.com/d/msgid/sympy/1004b533-73fe-4fc6-b051-a91b2a3c3fe9%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
def C_vm(L, iradius, oradius, density, basemass, eemass, Base, theta, thetad, phi, phid, n):
import math
import numpy as np
C_vm = np.matrix([[math.pi*density*(iradius - oradius)**2*((1/2)*math.sin(theta[0, 0])*L[1, 1] - 1/2*math.cos(theta[0, 0])*L[1, 0])*L[1, 0], math.pi*density*(iradius - oradius)**2*((1/2)*math.sin(theta[0, 0])*math.sin(theta[1, 0])*L[2, 0] + (1/2)*math.sin(theta[0, 0])*math.cos(theta[1, 0])*L[2, 1] + (1/2)*math.sin(theta[1, 0])*math.cos(theta[0, 0])*L[2, 1] - 1/2*math.cos(theta[0, 0])*math.cos(theta[1, 0])*L[2, 0])*L[2, 0], math.pi*density*(iradius - oradius)**2*((-math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*L[3, 1] + (math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*L[3, 0] + (1/2)*(math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*L[3, 0] + (1/2)*(math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*L[3, 1])*L[3, 0], math.pi*density*(iradius - oradius)**2*((3/2)*((-math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*L[4, 1] + (3/2)*(-(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) + (math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*L[4, 0] + (1/2)*(-(math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) + (math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*L[4, 1] + (1/2)*((math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (-math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*L[4, 0])*L[4, 0], eemass*((1/2)*(-(-(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) + (-math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) - math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.sin(theta[4, 0]) - math.cos(theta[0, 0])*math.cos(theta[4, 0]))*L[5, 0] + (1/2)*(((math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (-math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) - math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.cos(theta[4, 0]) + math.sin(theta[4, 0])*math.cos(theta[0, 0]))*L[5, 1] + (1/2)*(-(-(math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.sin(theta[4, 0]) + math.sin(theta[0, 0])*math.cos(theta[4, 0]))*L[5, 1] + (1/2)*((-(math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.cos(theta[4, 0]) + math.sin(theta[0, 0])*math.sin(theta[4, 0]))*L[5, 0] - (-(-math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - (math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.sin(theta[4, 0])*L[5, 1] + (-(-math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - (math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.cos(theta[4, 0])*L[5, 0] - ((1/2)*(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - 1/2*(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.sin(theta[4, 0])*L[5, 1] + ((1/2)*(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - 1/2*(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.cos(theta[4, 0])*L[5, 0])], [math.pi*density*(iradius - oradius)**2*(-1/2*math.sin(theta[0, 0])*L[1, 0] - 1/2*math.cos(theta[0, 0])*L[1, 1])*L[1, 0], math.pi*density*(iradius - oradius)**2*((1/2)*math.sin(theta[0, 0])*math.sin(theta[1, 0])*L[2, 1] - 1/2*math.sin(theta[0, 0])*math.cos(theta[1, 0])*L[2, 0] - 1/2*math.sin(theta[1, 0])*math.cos(theta[0, 0])*L[2, 0] - 1/2*math.cos(theta[0, 0])*math.cos(theta[1, 0])*L[2, 1])*L[2, 0], math.pi*density*(iradius - oradius)**2*((1/2)*(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*L[3, 0] + (1/2)*(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*L[3, 1] + (math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*L[3, 1] + (-math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*L[3, 0])*L[3, 0], math.pi*density*(iradius - oradius)**2*((1/2)*(-(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) + (math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*L[4, 1] + (1/2)*((math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (-math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) - math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*L[4, 0] + (3/2)*(-(-math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) + math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) + (-math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*L[4, 0] + (3/2)*((math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (-math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*L[4, 1])*L[4, 0], eemass*((1/2)*(-(-(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.sin(theta[4, 0]) - math.cos(theta[0, 0])*math.cos(theta[4, 0]))*L[5, 1] + (1/2)*((-(math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.sin(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (math.sin(theta[0, 0])*math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[0, 0])*math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.cos(theta[4, 0]) - math.sin(theta[4, 0])*math.cos(theta[0, 0]))*L[5, 0] + (1/2)*(-(-(-math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) + math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) + (math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.sin(theta[4, 0]) - math.sin(theta[0, 0])*math.cos(theta[4, 0]))*L[5, 0] + (1/2)*(((-math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) + math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.cos(theta[4, 0]) + math.sin(theta[0, 0])*math.sin(theta[4, 0]))*L[5, 1] - ((1/2)*(-math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) + math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - 1/2*(-math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.sin(theta[4, 0])*L[5, 1] + ((1/2)*(-math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) + math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - 1/2*(-math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.cos(theta[4, 0])*L[5, 0] - (-(math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - (-math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.sin(theta[4, 0])*L[5, 1] + (-(math.sin(theta[1, 0])*math.sin(theta[2, 0])*math.cos(theta[0, 0]) - math.cos(theta[0, 0])*math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - (-math.sin(theta[1, 0])*math.cos(theta[0, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[0, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.cos(theta[4, 0])*L[5, 0])], [0, math.pi*density*(iradius - oradius)**2*(-1/2*math.sin(theta[1, 0])*L[2, 1] + (1/2)*math.cos(theta[1, 0])*L[2, 0])*L[2, 0], math.pi*density*(iradius - oradius)**2*((-math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.cos(theta[1, 0])*math.cos(theta[2, 0]))*L[3, 0] + (-math.sin(theta[1, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[1, 0]))*L[3, 1])*L[3, 0], math.pi*density*(iradius - oradius)**2*((3/2)*(-(-math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) + (-math.sin(theta[1, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*L[4, 1] + (3/2)*((-math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*L[4, 0])*L[4, 0], eemass*(-(-(-math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (-math.sin(theta[1, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.sin(theta[4, 0])*L[5, 1] + (-(-math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) - (-math.sin(theta[1, 0])*math.cos(theta[2, 0]) - math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.cos(theta[4, 0])*L[5, 0] - (-1/2*(-math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) + (1/2)*(math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.sin(theta[4, 0])*L[5, 1] + (-1/2*(-math.sin(theta[1, 0])*math.sin(theta[2, 0]) + math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.cos(theta[3, 0]) + (1/2)*(math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.sin(theta[3, 0]))*math.cos(theta[4, 0])*L[5, 0] + (-1/2*(math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) + (1/2)*(math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.cos(theta[4, 0])*L[5, 1] - ((1/2)*(math.sin(theta[1, 0])*math.sin(theta[2, 0]) - math.cos(theta[1, 0])*math.cos(theta[2, 0]))*math.sin(theta[3, 0]) - 1/2*(math.sin(theta[1, 0])*math.cos(theta[2, 0]) + math.sin(theta[2, 0])*math.cos(theta[1, 0]))*math.cos(theta[3, 0]))*math.sin(theta[4, 0])*L[5, 0])]])
return(C_vm)
