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

游龙游龙

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
行列式求不出，画图找不到根
(*定义材料常数*)
  R; D$ e: s5 VSubscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
$ ~  C+ Q; c! p% ESubscript[\[Nu], LT] =
3 u3 _' u0 g9 C  I  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
  j$ x2 B( H8 _8 K( o) `\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
; ]; c+ c& ?" M, G9 qSubscript[C, L] = Subscript[EE, L]/\[Mu];
Subscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
Subscript[C, T] = Subscript[EE, T]/\[Mu];
3 d: K) H7 g" ISubscript[C, LT] = Subscript[G, LT];
(*脱层几何参数*)
: g" s! _7 f2 v+ N4 _' L" _\[Rho] = 1600;(*脱层密度*)
9 P# {8 u# H( U4 o% v: w3 F* Bc = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
l = 1;(*为梁总的长度*)
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
4 a, m$ V+ ^4 {% ZSubscript[l, 2] = 0.3;(*为梁2区的长度*)
5 H7 S; N7 i( c! aSubscript[l, 3] = 0.4;(*为梁3区的长度*)
Subscript[l, 4] = 0.4;(*为梁4区的长度*)
9 [9 X4 W9 \) c6 \0 C, P" jnn1 = h/t;(*nn1为总铺层厚度*)
4 S+ }! P4 p% k% E8 tnn2 = 2.0;(*为上子铺层层数*)
7 @$ ?! @' r' U3 j4 L6 a# b4 S/ cnn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
1 J6 K6 n9 _5 iSubscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
Subscript[h, 2] = h;(*梁2区的厚度*)
(*求各单层刚度*)
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)


ck = Cos[agk];(*列表求各铺层角度余弦*)
5 n% }) |7 \8 V5 Vsk = Sin[agk];(*列表求各铺层角度正弦*)
s11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
! i, T5 H9 \# h  I6 ^; |$ ~$ K             sk^4*Subscript[C, T] +
   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
A11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
A31 = (\!\(
. U9 q6 ^/ B% C\*UnderoverscriptBox[\(\[Sum]\), \(k =
" ~7 \6 q/ i; j  V     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
A41 = (\!\(
( h- I) f3 E& F- T( x\*UnderoverscriptBox[\(\[Sum]\), \(k =
     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)

5 x# }/ w7 B  e3 O+ s
1 V5 ^+ l6 t- k: ?4 o# dPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
; q  s( X- B  _  U* iD11 = \!\(
: Q% D8 B! R& R4 o\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
% ~0 s1 F; w3 Z" v/ i4 _8 m\*FractionBox[\(t^3\), \(12\)])\)\)
D21 = D11
D31 = \!\(
. ]# p) m, d& w* W\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
' @( U. n7 }, d* c\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
: \2 z) O! H9 h2 _& V2 {D41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
3 a# J, i: O" Y% w+ zk1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
- S* t+ E" b. q! o( _a1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
% ?% D9 X5 o# L' [0 o4 V* Q DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
: K0 b4 J! S* \- d9 g& w 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 \
; O- l; n8 k  d  l0 e1 D. Z. Q= 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];
/ ?+ `' H- S+ j* j& T) e# L% |, HDD1 = 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[
' ^; R8 p3 u3 T" ?( ~% o a2*\[Omega]^2, (4)^-1]; r13 = Power[
/ a# V# W/ q* t- h4 s/ l- U a3*\[Omega]^2, (4)^-1]; r14 = Power[a4*\[Omega]^2, (4)^-1];
  y/ r% T2 j) J& ABB = {{1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0},
- z; F6 `. n' k  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
/ D; o9 P+ m$ n  L# v8 G7 e   Sin[r12], 0, 0},
0 F' [( r, z6 z0 K! M  {0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
   Cos[r12], 0, 0},
1 ~. U. a1 I# k4 I8 `; v" O) e  {-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
8 ]0 t( @" f0 s7 ]   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
5 C7 Y# z3 F/ L3 M   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
' s# N1 k6 T& m* g8 @) `9 p  E/ {' G    DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
# S/ B6 r6 ]# f* x. Y4 E  {-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
- |( H0 ]; y- S- }   0, -(r11)^3*DD1*(\[Beta]1)^3*Cosh[r11],
   0, (r13)^3*DD3*(\[Beta]3)^3, (r14)^3*DD4*(\[Beta]4)^3, -(r11)^3*
+ Q+ w  }$ t' j+ r5 ^) ?    DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
. b! g+ c- n/ X9 ^   0, -(r13)^3*DD3*(\[Beta]3)^3, -(r14)^3*DD4*(\[Beta]4)^3},
* Z: k3 ]8 u+ C+ t; p8 |  {0, -(r12)^2*DD2*(\[Beta]2)^2, (r13)^2*DD3*(\[Beta]3)^2*
) y: p0 U4 ]. O. g! s    Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
+ ?  Z1 b+ M/ y. R9 o   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,
4 V0 y# n+ n/ j  g1 h9 x   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
0 u' [$ Y. K9 y' F) v. A% g    Cos[r14], 0, 0,
   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
. R$ M% t: j8 @+ C( k, R    Sin[r14]},
  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
5 C  S# p9 ]- x& m1 t    Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
7 Q0 \( m" i0 i4 N8 ]5 D   (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],
   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
: `: b9 I& l7 ?7 i+ P  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
% X8 y- e0 I( k' m& C% p' @   Sin[r11], 0, 0, 0},
! g& e& N/ V+ f; J! }$ q, y  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
   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,
/ R% O# K3 w% X$ _/ W   Sin[r11], 0, 0, 0},
! k7 ], E/ @, O3 O8 l5 [  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
3 F( S. Z* m: ~7 l' w   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]},
! R& \9 M9 Y0 r# f* f4 m  {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,
4 @5 _1 E. o: N; |    0, -Sin[r14]},
; _2 S+ c2 B. V0 Y+ Y  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
2 U! V1 B- L/ m$ ~   r14*b4*Sin[r14],
2 Z% k# s' y- m- t   0, r12, 0, -r14*b4*Cos[r14]}}
8 q. V2 t& a- h5 }aa1 = Det[BB]
  u" I0 L9 c; r, s4 J4 K5 LPlot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

6 A, K/ d9 ~, X) b5 J, U' r% o

百年孤独

这是什么语言编程？
, g- i6 ?* ?- ?7 k  |

游龙游龙

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

mathematical编的