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

游龙游龙

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
行列式求不出，画图找不到根
(*定义材料常数*)
0 ?: O0 i; y+ X/ E3 }* |9 _Subscript[\[Nu], TL] = 0.25;
. N3 G9 j( N& T$ Z% D1 SSubscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
; w4 H" Q0 C  N- ]Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
Subscript[\[Nu], LT] =
8 w2 K4 A6 U% Z: w8 n; Z  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
) f+ W" T4 l. ]1 U: o3 @# u8 p\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
# N# j4 }3 C: bSubscript[C, L] = Subscript[EE, L]/\[Mu];
& V/ [6 J2 P0 g- F( L3 GSubscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
Subscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
(*脱层几何参数*)
: m" Q8 R& t6 J/ s. Z5 }- X\[Rho] = 1600;(*脱层密度*)
" v$ K0 q% J6 J' Ac = 0.06;(*阻尼*)
- d! I6 Z. ^4 g/ e' v2 A4 a  k+ Yt = 0.01;(*t每层厚度*)
6 F& Q/ L3 e# q9 B( ~h = 0.05;(*h为梁总厚度*)
$ N$ X+ f. p# q0 Il = 1;(*为梁总的长度*)
* o$ `3 v& C4 hSubscript[l, 1] = 0.3;(*为梁1区的长度*)
3 S6 B2 X9 ?1 R9 hSubscript[l, 2] = 0.3;(*为梁2区的长度*)
7 G/ W4 C$ J4 X# k1 X) f- Q2 A  kSubscript[l, 3] = 0.4;(*为梁3区的长度*)
& S& w, `  i( [2 T  xSubscript[l, 4] = 0.4;(*为梁4区的长度*)
9 a6 _0 o) a3 K; p3 L2 S/ i; `5 hnn1 = h/t;(*nn1为总铺层厚度*)
nn2 = 2.0;(*为上子铺层层数*)
1 |- K$ z% e3 Hnn3 = h/t - nn2;(*下子铺层层数*)
8 e" S3 A' c# M# aSubscript[h, 3] = t*nn2;(*梁的上部子厚度*)
+ M0 @: O) e% F* E1 S% x4 B$ ^Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
Subscript[h, 2] = h;(*梁2区的厚度*)
! o  T9 ]% f5 A) S3 Q. L- l" c(*求各单层刚度*)
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)


ck = Cos[agk];(*列表求各铺层角度余弦*)
sk = 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 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
* i% a5 ^. P7 B    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
A31 = (\!\(
1 j( O6 \& t. _% x' r0 ?\*UnderoverscriptBox[\(\[Sum]\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
' u8 O1 d! y% d9 P* s9 ^" FA41 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
& V( v3 M0 ]" b' R& c5 I
* s2 Q& _; ~. x, H
( f  r( Y& t7 R4 _6 qPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
D11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
- M2 M0 T2 }9 Y' A\*FractionBox[\(t^3\), \(12\)])\)\)
D21 = D11
D31 = \!\(
: |9 x2 H' g3 F$ y' s+ @\*UnderoverscriptBox[\(\[Sum]\), \(k =
, M0 F! j4 a+ P2 O  i    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
/ Q( L4 b, m3 Z( s  h\*UnderoverscriptBox[\(\[Sum]\), \(k =
    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
4 h7 _" n3 f5 ?* u2 e9 D& N* `\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
# G& _  Y. ?, \+ L' }7 Ck1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
1 \  j1 A. ~9 e8 S/ j8 sa1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
1 I2 ?7 r6 |7 } DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
3 b; U; w8 o) e+ f4 q) @) } DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
2 s( a- n" `" H& P9 z\[Beta]1 = h/Subscript[l, 1]; \[Beta]2 = h/Subscript[l, 2]; \[Beta]3 \
= 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 = \
Subscript[l, 2]/Subscript[l, 4];
1 }: p! d) f: r% ODD1 = D11/(A11*h^2); DD2 = D21/(A11*h^2); DD3 = D31/(A11*h^2);
DD4 = D41/(A11*h^2);
6 Q3 Y6 f' M  H5 C* |# ir11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
& t3 F8 x. M5 ]) P9 b* d8 i5 ^0 L a2*\[Omega]^2, (4)^-1]; r13 = Power[
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},
+ z4 |+ N  C& a  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
% u3 v! |, i  A  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
; _3 p$ R. {* q0 w$ \* a   Sin[r12], 0, 0},
+ h9 z( X, C5 N+ h4 w  {0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
   Cos[r12], 0, 0},
  {-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
5 p- A- N3 X# [) B2 w   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
0 m2 ^1 v3 I1 V& r2 ?5 O2 Y3 |   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
5 u. V1 m/ S' J4 ]. T    DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
! @0 S5 ]8 F& y  {-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
   0, -(r11)^3*DD1*(\[Beta]1)^3*Cosh[r11],
   0, (r13)^3*DD3*(\[Beta]3)^3, (r14)^3*DD4*(\[Beta]4)^3, -(r11)^3*
, U  x& V1 e1 r( H$ x* v    DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
   0, -(r13)^3*DD3*(\[Beta]3)^3, -(r14)^3*DD4*(\[Beta]4)^3},
  b8 D1 d6 K( F0 A  {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*
    Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
   -(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*
7 w: {" g- j$ k, a) l    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,
8 U4 }' b- ~3 u5 S: `   (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],
# q# n+ M: B5 x# {   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
7 o" n, e4 c, P% G7 R0 [7 ?   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
5 e6 Z, b+ r, h" i. W% T  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
   Sin[r11], 0, 0, 0},
& _' J0 X2 }6 c/ v" c- B2 Q/ e: x5 _  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
, `# c0 Z1 W+ R3 G' p" E  V   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
- m9 `+ R/ @) O% Z0 C1 L   0, -r13*b1},
* q, Y' z, e4 g7 g2 T, c  {Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
3 c; o- N$ B$ G! l( u* U   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
* z* F& r+ O* x   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
6 @8 V, B/ ~: L+ f" V8 Z7 K, \5 `   0, 0, -r14*b2},
  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
    0, -Sin[r13]},
  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
   r13*b3*Sin[r13],
7 p  P+ b$ t2 v; }" K7 j   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,
" |# |0 E5 j' J* b5 M% g    0, -Sin[r14]},
# O; |1 y" H% v+ L! u7 @% h  {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]}}
7 N3 W8 N4 B, J. `/ Z* `4 eaa1 = Det[BB]
3 F' ]& @9 K* R7 h7 Y4 x* `0 \5 W: W  APlot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

; w3 x9 c) J. p- R

百年孤独

这是什么语言编程？
4 M/ g  k5 x3 N# W, C

游龙游龙

百年孤独 发表于 2014-7-8 09:30
& Y) L% F3 q. {3 q  G- @% B4 ~这是什么语言编程？

7 v* ~( K& D6 F( omathematical编的