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

游龙游龙

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
' `- G1 f  c5 s' z行列式求不出，画图找不到根
' c9 q4 w0 r% a. U7 B' g& r, ](*定义材料常数*)
4 ~* T+ z7 @  Z2 {& y6 G7 TSubscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
$ _% ~7 i; |& p7 [, X5 hSubscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
4 @2 z; c! n3 s- @; O: T4 [5 c4 nSubscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
Subscript[\[Nu], LT] =
1 r  L; x1 p8 _$ U8 G. v  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
9 ?# g. G. i3 SSubscript[C, L] = Subscript[EE, L]/\[Mu];
& V' E4 W$ G- W1 f' OSubscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
Subscript[C, T] = Subscript[EE, T]/\[Mu];
0 y! X9 O. I7 w& E1 ySubscript[C, LT] = Subscript[G, LT];
3 x3 ~1 [+ o$ v$ p6 ?" R$ O4 ^1 s(*脱层几何参数*)
\[Rho] = 1600;(*脱层密度*)
c = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
( Q1 A: q+ \. j: b% nl = 1;(*为梁总的长度*)
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
0 y7 t: ]# ~" U# e  }' W& b4 ZSubscript[l, 2] = 0.3;(*为梁2区的长度*)
4 _; h# H/ z8 c1 {& H+ YSubscript[l, 3] = 0.4;(*为梁3区的长度*)
Subscript[l, 4] = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
( Y/ t. A" r+ j/ Znn2 = 2.0;(*为上子铺层层数*)
& k& u. p# T+ r1 U; D( m0 q0 |nn3 = h/t - nn2;(*下子铺层层数*)
* b% D3 z4 b! C8 }: U1 {; w; ~1 eSubscript[h, 3] = t*nn2;(*梁的上部子厚度*)
4 l$ K+ w# V/ D; K8 V$ a  USubscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
Subscript[h, 2] = h;(*梁2区的厚度*)
2 c: z2 a- Y. y( I0 L0 R3 T(*求各单层刚度*)
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)
  d0 b# i- f& M2 _9 e' o
! |8 M$ J  d) a$ {6 a
ck = Cos[agk];(*列表求各铺层角度余弦*)
% z/ X/ d2 H) G8 jsk = Sin[agk];(*列表求各铺层角度正弦*)
s11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
             sk^4*Subscript[C, T] +
/ |- P7 Y# c1 u- O3 u" z   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
A11 = \!\(
; \( T6 ^5 T' q& Z2 w\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
+ O% `3 {2 R, D5 fA21 = A11
0 Y% z! h, F# Q' WA31 = (\!\(
  c9 c. e; t- C' d3 N\*UnderoverscriptBox[\(\[Sum]\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
A41 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
8 V! X$ s* P' @. a  N& Y; e+ q
( m7 i1 v5 r3 O/ ~
5 R7 m8 a2 U2 Z: g$ mPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
8 J! H, D! u$ j; _4 |& Z. XD11 = \!\(
( |$ T; v% V' [2 j0 A+ \\*UnderoverscriptBox[\(\[Sum]\), \(k =
9 V2 s! u" y! |0 I& Q    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
" X/ S) R* K: x2 V" X+ kD21 = D11
2 m1 u" C2 u* ]' k" G9 L2 LD31 = \!\(
  N- M  I4 ^6 x; b\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
; k. x( Z( T0 m" I\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
8 K" F: l' h& p; K3 w$ x\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
: C9 D  @* ~) \0 E4 Y3 a! l7 `k1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
a1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
8 X$ M) P4 h) }2 j- j2 Y 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 \
* i1 n9 }2 A2 O5 ~$ A= h/Subscript[l, 3]; \[Beta]4 = h/Subscript[l, 4];
  l$ m0 W. s1 B0 c" Mb1 = Subscript[l, 1]/Subscript[l, 3]; b2 = Subscript[l, \
( Q5 [$ W4 n& r5 `  k& I/ \6 h1]/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);
% B" ~2 }% K) g+ H5 L) Y4 |DD4 = D41/(A11*h^2);
r11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
5 ~" M: p4 e2 @$ I a2*\[Omega]^2, (4)^-1]; r13 = Power[
4 U' }& ^, Y7 s. a1 _1 s* `- h( 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},
2 e; a. }5 `6 g' t5 a7 M" z$ s& A+ q  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
& t# c0 M, W8 e9 |   Sin[r12], 0, 0},
- ~* L" W4 Y3 I* i( v$ {, y" u( m3 K3 @  {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],
   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
$ f- w6 u3 M2 l) T5 x/ i# S& n   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
- |7 U6 P4 C) \% q/ T; y3 j9 @0 `    DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
  {-(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*
    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},
  {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,
; E5 m/ n) I  O3 i0 L9 X$ o   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
    Cos[r14], 0, 0,
! l8 ?2 A, i, ^/ G6 p  X% 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, W( Y2 h3 S3 g' U    Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
   (r13)^3*DD3*(\[Beta]3)^3*Cosh[r13], (r14)^3*DD4*(\[Beta]4)^3*
: _! c2 u) t7 t  J: K' _    Cosh[r14], 0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sin[r13],
/ f  x  S. t* H3 g0 ~1 a# c, G   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
$ ?9 u  L2 O+ c/ f; i. t   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
: ~7 @3 B6 s$ v. r' e; p5 [1 t  b" {   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
: z4 K; [1 e# P2 ?& f; l  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
' o- M7 U% F) `$ s; ^& O   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
   0, -r13*b1},
  {Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
   Sin[r11], 0, 0, 0},
9 R  P; j, t" W0 a+ b% g  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
0 k1 L; m/ ], ?' f   0, 0, -r14*b2},
  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
    0, -Sin[r13]},
* P3 Q. F& y6 ?1 q) J  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
$ w: Q9 b% p4 `" T   r13*b3*Sin[r13],
" f( b0 k7 ?. x   0, 0, r12, -r13*b3*Cos[r13], 0},
! r5 j. t) g3 k1 v) L1 L! ~. K  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
+ a! @% c8 _0 F    0, -Sin[r14]},
  {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]}}
3 J# j: |6 A) ~$ Aaa1 = Det[BB]
1 ^" t* z8 F4 g2 W, ePlot[aa1, {\[Omega], -10, 10}]
4 j9 U4 u9 m2 w% |9 d* |5 }, @" ~ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

百年孤独

这是什么语言编程？
% |5 ]% r3 N8 z$ b6 U+ Y% x# b

游龙游龙

百年孤独 发表于 2014-7-8 09:30
1 h  H1 A2 Z7 U0 V, [8 N: I8 d. D" E这是什么语言编程？

2 n- P) Q9 `8 K3 q: f: Lmathematical编的