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

游龙游龙

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
行列式求不出，画图找不到根
(*定义材料常数*)
9 E+ Z: i9 S/ {. l% aSubscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
) X- S+ w6 _. r  @" @Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
Subscript[\[Nu], LT] =
  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
2 ]& I9 F% Y- l8 ?& l4 p\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
Subscript[C, L] = Subscript[EE, L]/\[Mu];
5 w' _7 N9 o+ F: V) k* tSubscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
Subscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
(*脱层几何参数*)
\[Rho] = 1600;(*脱层密度*)
c = 0.06;(*阻尼*)
7 B" x1 b; ?+ A' Jt = 0.01;(*t每层厚度*)
! ~" p% H! A- H8 [, G( @h = 0.05;(*h为梁总厚度*)
* a, i6 J6 I! u! y% i# v0 Yl = 1;(*为梁总的长度*)
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
Subscript[l, 2] = 0.3;(*为梁2区的长度*)
/ P- U6 k6 M, Y. }: b0 WSubscript[l, 3] = 0.4;(*为梁3区的长度*)
Subscript[l, 4] = 0.4;(*为梁4区的长度*)
# ~9 R, E% P# t  N/ B% |nn1 = h/t;(*nn1为总铺层厚度*)
nn2 = 2.0;(*为上子铺层层数*)
! y6 |9 J. N; h. L" fnn3 = h/t - nn2;(*下子铺层层数*)
* q+ J& J9 x, p4 WSubscript[h, 3] = t*nn2;(*梁的上部子厚度*)
Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
- f5 I" e% |3 b( OSubscript[h, 2] = h;(*梁2区的厚度*)
(*求各单层刚度*)
( m. S- J$ `! z3 _agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)


ck = Cos[agk];(*列表求各铺层角度余弦*)
" D, y( H- ]2 u, ?" Wsk = Sin[agk];(*列表求各铺层角度正弦*)
s11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
0 R. J, o0 k# O             sk^4*Subscript[C, T] +
   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
A11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
8 M# V% h9 Q8 bA21 = A11
$ I7 T/ Z. d2 [4 W2 }  YA31 = (\!\(
6 s* A; X) E& E2 t1 k- e\*UnderoverscriptBox[\(\[Sum]\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
0 v/ }- \! _$ v  A0 e* LA41 = (\!\(
: R) S% }8 h: O# f: |\*UnderoverscriptBox[\(\[Sum]\), \(k =
     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)

# w0 P2 f. |  o
7 `. N% h) }* y. U. Q) {Print["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
D11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
* [' w% Q% u7 u. D3 {5 ]! e$ ED21 = D11
D31 = \!\(
6 _1 ~7 a1 D0 H6 I% Z" w/ U\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
% z( G- R, H# A. i5 X$ o) XD41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
0 A  S& u; e! E9 t\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
k1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
* D) v% N$ D) ?3 [7 G3]/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)/(
DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
, t# e8 k/ S$ o* c2 |  R3 J: Q7 k, 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];
- ]& Q3 U( T1 [( Zb1 = Subscript[l, 1]/Subscript[l, 3]; b2 = Subscript[l, \
+ ~/ b" Z9 T0 d+ m1 k8 T5 M) k1]/Subscript[l, 4]; b3 = Subscript[l, 2]/Subscript[l, 3]; b4 = \
7 H* {, y1 [$ s/ ^% A2 ySubscript[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);
5 p" g* s4 n" o4 R0 W1 i) ^r11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
% s9 i9 s. V% H a2*\[Omega]^2, (4)^-1]; r13 = Power[
$ S# @- ^; L  S& C" ` a3*\[Omega]^2, (4)^-1]; r14 = Power[a4*\[Omega]^2, (4)^-1];
! _/ m8 z/ o8 o& ]9 XBB = {{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},
  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
3 L1 A5 S6 f" E) F   Sin[r12], 0, 0},
  {0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
5 |$ |" K% H+ b. |$ c* L! r   Cos[r12], 0, 0},
9 b" Y; q2 Q  O' o0 S6 a5 I  {-(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,
   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
, G; w$ V4 K8 e! m6 y. D* W& Y+ }   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
; g- o3 o# u  j( J; t$ X4 Q    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*
. s  Q) J7 U: h" O- B, Z    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},
& W& j4 w, p- _; S7 H# b/ c* W  {0, -(r12)^2*DD2*(\[Beta]2)^2, (r13)^2*DD3*(\[Beta]3)^2*
" A- M) {  t5 h+ M6 i) e% i2 w    Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
8 e7 L& X; V( _- {9 T- S& v   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,
" |8 b8 E4 L; g3 G" }4 _   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
3 @6 y2 m5 b8 ?, I    Cos[r14], 0, 0,
   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
' A# z; B( D. D( _9 C( ~; G" |    Sin[r14]},
; f) `8 p( |1 h: m  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
1 w3 n! a3 [7 @- ~& ?) V    Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
" ?, z$ k% W, Z2 k0 y0 ?! @; a   (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],
+ Y- }) O, }4 A- X3 p# U   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
5 X6 j+ ]. N; c* p, N  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
; t1 ^2 {  W4 r: O1 k: N   Sin[r11], 0, 0, 0},
: ?! S" T. `# V6 m8 J  {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,
4 W! i. }7 A: L' M" F7 B# A- I8 a   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
& u! @/ E2 t- h6 y8 m( D# B4 T   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
0 N% i" G0 V  m, E( ~+ }2 s   0, 0, -r14*b2},
  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
; B8 n' E# r! x5 e    0, -Sin[r13]},
" K1 w8 `- }% M  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
5 D' u$ o/ Q1 t: g8 ]0 f   r13*b3*Sin[r13],
   0, 0, r12, -r13*b3*Cos[r13], 0},
! {  N1 U9 i& u! ^1 w7 R0 b  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
, A* ~# ?9 p! }5 p    0, -Sin[r14]},
. v" r/ v0 X" w5 b  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
; H; e4 y' b! b3 P( G   r14*b4*Sin[r14],
7 G9 D9 l  r& F1 d   0, r12, 0, -r14*b4*Cos[r14]}}
aa1 = Det[BB]
0 I. W7 o0 H# I) E0 I( o! ]Plot[aa1, {\[Omega], -10, 10}]
" o6 f' j$ L, d9 z& ?2 xss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

百年孤独

这是什么语言编程？
/ S/ L0 ]. O* a& g2 h) B

游龙游龙

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

mathematical编的