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

游龙游龙

2 z* @: B9 j" t5 u* Z- y1 ^) i麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
1 `# g& U! |4 v4 F% r% x9 I行列式求不出，画图找不到根
4 n& c( u# Y) z% G7 Y8 d) O(*定义材料常数*)
Subscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
8 Y$ [7 N* |1 T3 KSubscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
0 n6 x0 ]! g' u: U5 j: E+ vSubscript[\[Nu], LT] =
$ j# `: g' h: y0 E7 n2 }# f  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
\[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];
Subscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
: B) R4 a" C( {2 L6 K1 z(*脱层几何参数*)
% t( r5 [5 u0 S0 }! L# \\[Rho] = 1600;(*脱层密度*)
c = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
' N9 J, b4 [/ q1 A, K' ?3 rl = 1;(*为梁总的长度*)
3 s9 Y4 I* n. M* s- O! Z/ pSubscript[l, 1] = 0.3;(*为梁1区的长度*)
Subscript[l, 2] = 0.3;(*为梁2区的长度*)
) X' p$ n1 ?3 n- c! Y1 L: }Subscript[l, 3] = 0.4;(*为梁3区的长度*)
8 _: }) b) U; w! p8 RSubscript[l, 4] = 0.4;(*为梁4区的长度*)
+ T7 H% D& e7 }$ q" n, g: Onn1 = h/t;(*nn1为总铺层厚度*)
nn2 = 2.0;(*为上子铺层层数*)
nn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
) @/ n. V, x( k/ i+ s( ]/ oSubscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
: a( \2 ^+ ~- H% s7 q$ _1 TSubscript[h, 1] = h;(*梁1区的厚度*)
Subscript[h, 2] = h;(*梁2区的厚度*)
/ Y0 p. [) e* ~; D5 P7 a0 @(*求各单层刚度*)
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)
+ W/ \4 B1 ?8 r# @- W

ck = Cos[agk];(*列表求各铺层角度余弦*)
sk = Sin[agk];(*列表求各铺层角度正弦*)
5 S! T+ U* W% y8 B' d; D9 v7 ?8 Js11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
             sk^4*Subscript[C, T] +
& n  \5 W8 Z: v; q8 W* l  P3 Y   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
A11 = \!\(
$ q5 B3 x: K5 n\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
: K1 D6 h2 g: w0 SA31 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
" X6 f. Z; h6 Y# z% I4 K) x% E     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
/ z3 h. A' t+ w: R- n, qA41 = (\!\(
$ x& R/ S' C& N; @, _; W& T; Z. C\*UnderoverscriptBox[\(\[Sum]\), \(k =
     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
  G) ~4 {+ g3 f& r7 }" {
4 t- o: Q" N" y1 D* i" H
* ^/ A' q# p# I* ?. z; bPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
4 j* v) |- f. H8 vD11 = \!\(
# G6 {, E: n2 b\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
, C2 K9 H1 Q. I\*FractionBox[\(t^3\), \(12\)])\)\)
; }0 Y. ]3 c8 {+ \D21 = D11
. E  q2 R( r7 e+ D7 `; k  z8 ZD31 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
, c4 x3 {3 n' n# B- X    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
8 u: F( G$ j4 zD41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
! H. Y, c. b7 z: o0 ^    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
# u# b, w  P: `. D; Ik1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
3 o) Q0 o' A$ v8 |) ]4 a6 Z2 u) z# ba1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
4 ~7 o" n2 N" a; U DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
" m5 x0 |- Z; W; U\[Beta]1 = h/Subscript[l, 1]; \[Beta]2 = h/Subscript[l, 2]; \[Beta]3 \
= h/Subscript[l, 3]; \[Beta]4 = h/Subscript[l, 4];
, Z9 b$ b: \& E- S" |b1 = Subscript[l, 1]/Subscript[l, 3]; b2 = Subscript[l, \
* Y: I5 ?7 n( \' E$ n& K  \6 k4 O% F4 C1]/Subscript[l, 4]; b3 = Subscript[l, 2]/Subscript[l, 3]; b4 = \
4 t; k$ w8 z4 ~' g  WSubscript[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);
. p$ F% O. _3 B7 L9 O- |r11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
a2*\[Omega]^2, (4)^-1]; r13 = Power[
$ d. }6 Y. e' t6 [ 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},
  I/ y- x4 j- }& S  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
   Sin[r12], 0, 0},
' p3 J) a, u4 @# _  {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],
+ i8 R6 M! {" Z8 n& v( a   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
% C) W8 P5 t, K& }# v# t   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
    DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
. K9 ^/ e) G; z7 Q* w  {-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
0 x5 P2 A4 w6 I- y# H! C   0, -(r11)^3*DD1*(\[Beta]1)^3*Cosh[r11],
1 V/ \+ y' {0 B- u9 {2 C, w& L3 D   0, (r13)^3*DD3*(\[Beta]3)^3, (r14)^3*DD4*(\[Beta]4)^3, -(r11)^3*
# N  l2 c) I& k1 |    DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
' ]3 s0 I9 r. A$ H   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
- z8 e4 Z+ c6 s) k, I5 L1 o* w/ h7 N   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*
- s2 [' u; ?, }, ~9 A7 a    Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
1 m: u8 m5 S3 ^9 u! [( n   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,
* `1 \, o$ ]$ g5 O   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
$ @7 k4 i9 M& }' J9 y    Cos[r14], 0, 0,
# \' U) g5 i4 K, [/ m   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
8 M2 z. b+ o1 a! {$ e: _    Sin[r14]},
  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
9 P5 E6 w" k3 v    Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
% c" ]5 O+ n9 y* z   (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],
   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
+ a% o" ?& w3 E" T   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
4 |! r6 i7 x! u" u! q9 h  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
   Sin[r11], 0, 0, 0},
3 Q# V' p! z$ W1 e" K  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
3 h: c. q* E. M+ t: ~! c* m" u  X" l   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
" J8 u: W; W; ^+ I. Z& ?0 D# O4 a: ~   0, -r13*b1},
  {Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
   Sin[r11], 0, 0, 0},
! c. s% u7 q9 ?$ F  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
   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,
* @  ]% j3 v: B) d- w! Z2 v    0, -Sin[r13]},
; U1 j  u: P3 J0 w  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
, u3 [8 f- z8 F" |3 s% d% n4 f   r13*b3*Sin[r13],
4 z7 U7 a. C$ }/ q. @   0, 0, r12, -r13*b3*Cos[r13], 0},
- ]" X# X; Z& A, @4 D0 F0 Y  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
    0, -Sin[r14]},
  z. _2 @+ \8 N% [  {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]
Plot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

游龙游龙

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

; }% R; v% y% f& P5 j7 F" h) s& f- Z; lmathematical编的

百年孤独

这是什么语言编程？