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

游龙游龙

( i& O+ |1 O" Z% P* ~# O麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
行列式求不出，画图找不到根
(*定义材料常数*)
2 {0 t: ?- P% c' Z) H1 cSubscript[\[Nu], TL] = 0.25;
" g) [! t: v  r" C+ O' x: |5 }Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
6 }+ i. U  \6 p6 DSubscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
7 V& r+ M% x1 I* f+ GSubscript[\[Nu], LT] =
  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
* `0 @8 U! q2 b/ o\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
Subscript[C, L] = Subscript[EE, L]/\[Mu];
Subscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
) n  _, v0 v9 m0 h) \) fSubscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
' X2 [0 W6 @: ]" ](*脱层几何参数*)
5 e0 W, N: A7 Q8 O\[Rho] = 1600;(*脱层密度*)
c = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
" L( e* |; c) a( }* \l = 1;(*为梁总的长度*)
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
Subscript[l, 2] = 0.3;(*为梁2区的长度*)
) u' S, c; ]3 a% q6 eSubscript[l, 3] = 0.4;(*为梁3区的长度*)
) p& n3 r, n7 e5 _; {- O6 e" s+ ~Subscript[l, 4] = 0.4;(*为梁4区的长度*)
9 L7 C% M4 e1 p$ a9 fnn1 = h/t;(*nn1为总铺层厚度*)
" i4 K  z2 O( j8 s  H) \# f# a5 unn2 = 2.0;(*为上子铺层层数*)
5 X  F$ t4 j5 t0 _! i6 E  Cnn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
( K" e% X  l5 z- A% f' TSubscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
9 n; K. K# L% A# a4 w& p  PSubscript[h, 1] = h;(*梁1区的厚度*)
Subscript[h, 2] = h;(*梁2区的厚度*)
& S( F* m+ u4 e/ `/ [/ _9 t(*求各单层刚度*)
7 s6 x3 d  U9 U4 w- l2 ragk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)


; w* B% @  a, t7 I0 {# F' ~' i( sck = Cos[agk];(*列表求各铺层角度余弦*)
4 O) j1 U: y! Q9 o. wsk = Sin[agk];(*列表求各铺层角度正弦*)
s11k = 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];(*求各区刚度*)
A11 = \!\(
  I; \- i4 D% \* |2 B\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
: O" J# X* B: IA21 = A11
A31 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
1 P( p" p% b6 c; B+ l# \+ w6 h% g     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
A41 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
. J! f& b0 Z5 A& r2 c3 o3 `' s" ~
/ }6 e# [3 d- O
1 l" \2 U6 j# \% ePrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
D11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
0 Y+ r) x0 m/ u  DD21 = D11
D31 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
- P8 M+ Y3 B0 s8 J    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
3 n' O( j! t4 R& f- S% M) G" [; z\*UnderoverscriptBox[\(\[Sum]\), \(k =
* c, D9 h+ N. M1 f% c    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
- w) A* ]2 Q2 b% W& E0 \# Y# h, W\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
! \/ o2 p  |& @0 u" G7 K9 Rk1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
+ v9 `$ i3 m6 v9 Z+ N# ea1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
\[Beta]1 = h/Subscript[l, 1]; \[Beta]2 = h/Subscript[l, 2]; \[Beta]3 \
9 g. T  J# t) P; ^% {$ _= h/Subscript[l, 3]; \[Beta]4 = h/Subscript[l, 4];
b1 = Subscript[l, 1]/Subscript[l, 3]; b2 = Subscript[l, \
1]/Subscript[l, 4]; b3 = Subscript[l, 2]/Subscript[l, 3]; b4 = \
4 L' C) f* q6 d5 O- |0 ASubscript[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);
" I" Z) U$ \' ^* M% S( o, M( S% hr11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
a2*\[Omega]^2, (4)^-1]; r13 = Power[
$ P* b3 X0 _# V/ o% {: {% n) k% L; ^: t a3*\[Omega]^2, (4)^-1]; r14 = Power[a4*\[Omega]^2, (4)^-1];
BB = {{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 M/ j6 M  W1 {: A; o  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
0 _" l" R: {' K% Q   Sin[r12], 0, 0},
9 `/ B/ a8 P7 O& j6 r  {0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
   Cos[r12], 0, 0},
- l4 e; ]* J5 u- w$ l  {-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
8 f4 S+ Y" w! ^   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
/ V" B5 E3 l  l+ P   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
  X' P% Y+ W, ^* ~" J9 m    DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
" f3 ]7 d! g8 l4 s- k/ r  {-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
   0, -(r11)^3*DD1*(\[Beta]1)^3*Cosh[r11],
) D. x" _& E+ e   0, (r13)^3*DD3*(\[Beta]3)^3, (r14)^3*DD4*(\[Beta]4)^3, -(r11)^3*
    DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
8 }) _* f2 i1 o8 w   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
   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*
: H3 w. ?+ D2 ^( \! h    Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
6 A! z5 J3 D! R   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
    Cos[r14], 0, 0,
   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
- l2 Q: E. _# e5 L: X$ O    Sin[r14]},
  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
    Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
   (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],
" J  S2 e* B: F, t9 H+ o4 J6 S   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
+ h8 n' t; H, a) Z   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
/ L$ `# S) j  z   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
) F2 ?: z3 I3 `   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
! {4 F& Q8 x! T, r   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
! J- |* `4 a$ M6 L3 H/ @6 H( c  z   0, -r13*b1},
- g! l; Y) U0 m2 z- I" f% C  {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,
: J* G: `" K+ O+ x$ `$ Z5 {   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
   0, 0, -r14*b2},
  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
    0, -Sin[r13]},
6 J3 }7 T. x) i0 J3 T0 N3 c+ Z  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
+ G% y8 a9 c8 r8 k   r13*b3*Sin[r13],
   0, 0, r12, -r13*b3*Cos[r13], 0},
1 X' k% x. k( c# u  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
9 K0 A' U- C* `3 w  `    0, -Sin[r14]},
! e' n% P; X" m3 Q" g  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
   r14*b4*Sin[r14],
# k& {( J* q, _* _. U   0, r12, 0, -r14*b4*Cos[r14]}}
aa1 = Det[BB]
Plot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]
. `% y1 H5 G/ n8 f+ j- w
! }8 B. M* _- T+ n  L

百年孤独

这是什么语言编程？

游龙游龙

百年孤独 发表于 2014-7-8 09:30
这是什么语言编程？

* ~. ]6 R2 M" U6 P# i! p+ h8 dmathematical编的