麻烦大家帮忙看下mathematical编的程序，不知道哪里出错了

游龙游龙

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
行列式求不出，画图找不到根
9 _/ ]' i0 O' h  f( P7 b(*定义材料常数*)
Subscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
- U7 T4 m: P. `6 w* pSubscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
Subscript[\[Nu], LT] =
! Q* Z3 j, x+ g  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
5 h- y3 i0 u$ X4 a, p/ nSubscript[C, L] = Subscript[EE, L]/\[Mu];
7 k# f- L+ N- e, D- w! @! X$ J% P& CSubscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
: o8 O9 F, A& ~; ]Subscript[C, T] = Subscript[EE, T]/\[Mu];
5 I  _/ c9 H+ y1 g: p: dSubscript[C, LT] = Subscript[G, LT];
(*脱层几何参数*)
  H+ c1 w, H/ {$ L1 j5 i0 P\[Rho] = 1600;(*脱层密度*)
c = 0.06;(*阻尼*)
) k0 B! W) Q" B8 X# C! mt = 0.01;(*t每层厚度*)
! t) U+ N2 f7 i& Z( r2 Fh = 0.05;(*h为梁总厚度*)
l = 1;(*为梁总的长度*)
% Q8 Z+ Y# Q$ ^Subscript[l, 1] = 0.3;(*为梁1区的长度*)
1 J( v( [( c5 z1 Z1 @Subscript[l, 2] = 0.3;(*为梁2区的长度*)
) v9 Z$ n1 B1 ~, ]) ESubscript[l, 3] = 0.4;(*为梁3区的长度*)
5 e( d7 Z+ c. O& G2 l6 j* PSubscript[l, 4] = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
nn2 = 2.0;(*为上子铺层层数*)
nn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
$ O; m" `3 p4 |, t& VSubscript[h, 2] = h;(*梁2区的厚度*)
(*求各单层刚度*)
4 \1 n) S# U2 Y" ?/ qagk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)


. g. i4 K6 p& j* h2 s6 [3 Uck = Cos[agk];(*列表求各铺层角度余弦*)
5 J& _; _. z+ R0 ask = Sin[agk];(*列表求各铺层角度正弦*)
+ M4 P1 Z& M9 k! Y9 }& Ls11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
             sk^4*Subscript[C, T] +
   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
4 K2 t6 o! `3 f) _( ?- f, mA11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
  e/ h" a( G) X% y7 P5 b    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
; u6 S4 Q5 {- }# C- @  {A31 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
8 M3 n6 U3 P# |- b5 ^     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
4 m: Q/ A! k8 h5 CA41 = (\!\(
9 L% q# T& D3 ^( p3 y' n4 X\*UnderoverscriptBox[\(\[Sum]\), \(k =
0 A$ S1 i8 m1 a) e     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
0 U& F0 h; |  z9 T& `

' M& M. d- Z  H$ k$ k8 J  u+ IPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
' t% C2 e1 H8 Y5 t+ y: Z$ gD11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
( @) o% Q3 n/ ^; Q) ~\*FractionBox[\(t^3\), \(12\)])\)\)
D21 = D11
D31 = \!\(
3 U- \$ G! i6 i1 y8 Z1 i! M1 }\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
2 g& c0 M# E0 |; B! {; [7 E' A- B7 p\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
- l8 g- x( S# h& C' N9 @* V    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
k1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
3 v. h4 \$ b4 X5 _/ E# t0 y; w; Va1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
& y9 D! t& J. W$ O2 ?# u DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
- ~1 x; E" D4 u; J' {: m" n DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
: L( O6 Q: D$ [2 k! ]$ H) u" j5 O\[Beta]1 = h/Subscript[l, 1]; \[Beta]2 = h/Subscript[l, 2]; \[Beta]3 \
= h/Subscript[l, 3]; \[Beta]4 = h/Subscript[l, 4];
. A3 c+ L6 Q5 d7 `+ q. Ab1 = Subscript[l, 1]/Subscript[l, 3]; b2 = Subscript[l, \
1]/Subscript[l, 4]; b3 = Subscript[l, 2]/Subscript[l, 3]; b4 = \
Subscript[l, 2]/Subscript[l, 4];
DD1 = D11/(A11*h^2); DD2 = D21/(A11*h^2); DD3 = D31/(A11*h^2);
DD4 = D41/(A11*h^2);
r11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
; ^- s" `* I$ l. D9 n+ O a2*\[Omega]^2, (4)^-1]; r13 = Power[
a3*\[Omega]^2, (4)^-1]; r14 = Power[a4*\[Omega]^2, (4)^-1];
, a! `9 U2 K$ bBB = {{1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
8 G  ^' b$ S6 l+ K; v3 u, I6 h, _  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
   Sin[r12], 0, 0},
  {0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
   Cos[r12], 0, 0},
$ h* j4 v8 R% v2 h/ t3 Y  {-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
- f$ x# f2 J3 _+ m8 m8 N   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
1 Q( B: ]/ L& w   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
1 ]/ ^* G5 S! ?7 ?% V. E1 W   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
; B( c# w1 ]  a3 c5 s4 m& [    DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
  {-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
4 U$ t, l1 K$ P" p) ]  f   0, -(r11)^3*DD1*(\[Beta]1)^3*Cosh[r11],
   0, (r13)^3*DD3*(\[Beta]3)^3, (r14)^3*DD4*(\[Beta]4)^3, -(r11)^3*
( T2 k# z$ q( F+ U2 G9 Y4 d    DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
2 L: @/ P1 p+ I7 O% G0 f. p   0, -(r13)^3*DD3*(\[Beta]3)^3, -(r14)^3*DD4*(\[Beta]4)^3},
  {0, -(r12)^2*DD2*(\[Beta]2)^2, (r13)^2*DD3*(\[Beta]3)^2*
    Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
   0, 0, (r13)^2*DD3*(\[Beta]3)^2*Sinh[r13], (r14)^2*DD4*(\[Beta]4)^2*
" a- B8 y8 K3 E  f* Q# {. W5 h" `; p! G: j    Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
4 w' E. F( t; x6 q4 {5 G   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
    Cos[r14], 0, 0,
" @: T2 M' b$ @4 L   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
    Sin[r14]},
  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
) N  T" z: |6 [, ]    Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
' v% L/ \3 m! I/ {   (r13)^3*DD3*(\[Beta]3)^3*Cosh[r13], (r14)^3*DD4*(\[Beta]4)^3*
    Cosh[r14], 0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sin[r13],
, S, k( B3 l2 x+ h% y9 O   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
& }# Z/ _7 p# \   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
% q9 E9 T$ J: T$ ^   Sin[r11], 0, 0, 0},
2 C" r0 d6 M( d( c/ H  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
) ?' S# z2 P" N) v  k& `. Z( ]8 T   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
% P! p7 G7 i. B+ q. L+ B   0, -r13*b1},
  {Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
& M- |( m5 F3 ~% L7 |2 ^   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
! y. K% ?9 n9 e& v   0, 0, -r14*b2},
$ l" m0 t2 C4 N$ S  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
2 g$ l4 R- h9 E, r0 @2 Z    0, -Sin[r13]},
  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
   r13*b3*Sin[r13],
   0, 0, r12, -r13*b3*Cos[r13], 0},
  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
- p- j5 Z! ?  _0 t8 l9 K8 C    0, -Sin[r14]},
) y/ B2 e0 A, f  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
   r14*b4*Sin[r14],
   0, r12, 0, -r14*b4*Cos[r14]}}
aa1 = Det[BB]
0 i/ a. u  b8 K/ X. ~Plot[aa1, {\[Omega], -10, 10}]
. I1 _) }$ Y" z$ G( s$ wss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

8 v/ A1 j/ x: I; B' |& x/ D7 T

百年孤独

这是什么语言编程？

游龙游龙

百年孤独 发表于 2014-7-8 09:30
) Q) t7 |- v0 w1 t这是什么语言编程？

" t$ w1 \2 `/ }mathematical编的