Hi,
I try to solve system of symbolic equatios:
#------------------------------------------------------------------------------------------------------------------------------------
var ("A_0, B_0, C_0, D_0, r, h, E, q, v, e11, k11, m11, n11, Q")
v=0.3
#E=2*10^11
q=10*10^3
r_1=.1
r_2=.2
#h=0.004
D=(E*(h^3))/(12*(1-v^2))
W_0=(q*r^4)/(64*D)
f_0=A_0*r^2*ln(r) + B_0*r^2 + C_0*ln(r) + D_0
W=W_0+f_0
W1=(e11*r^4+k11*r^(-4)+m11*r^6+n11*r^(-2))*cos(4*Q)
W_diff_1 = diff(W,r,1)
W_diff_2 = diff(W,r,2)
W1_diff_1 = diff(W1,r,1)
W1_diff_2 = diff(W1,r,2)
M=(-D)*(W_diff_2 + v*1/r*W_diff_1)
MM=D*(W1_diff_2 + v*(-12*e11*r^2 - 20*k11*r^(-6) - 10*m11*r^4 -
18*n11*r^-4)*cos(4*Q))*cos(4*Q)
M_diff_1 = diff(M,r,1)
eq1=(M(r=r_2)==0)
eq2=(M(r=r_1)==0)
eq3=(W(r=r_1)==0)
eq4=(W_diff_1(r=r_1)==0)
eq5=(M(r=r_2)+MM(r=r_2)==-cos(4*Q)*r_2*M_diff_1(r=r_2)-M(r=r_2))
eq6=(M(r=r_1)+MM(r=r_1)==0)
eq7=(W1(r=r_2)==0)
eq8=(W1(r=r_2)==-cos(4*Q)*r*W_diff_1(r=r_2))
solve([eq1, eq2, eq3, eq4, eq5, eq6, eq7, eq8], A_0, B_0, C_0, D_0,
e11, k11, m11, n11)
#-----------------------------------------------------------------------------------------------------------------------------------------
and get an error:
raise ValueError, "Unable to solve %s for %s"%(f, args)
ValueError: Unable to solve "my_system_of_symbolic_equations"
Can anybody help me to do it?
I dont't understand what is my mistake...
--~--~---------~--~----~------------~-------~--~----~
To post to this group, send email to [email protected]
To unsubscribe from this group, send email to
[email protected]
For more options, visit this group at
http://groups.google.com/group/sage-support
URLs: http://www.sagemath.org
-~----------~----~----~----~------~----~------~--~---
