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

游龙游龙

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
) L1 L9 ^7 ]8 M! C3 P3 _( p( \: z" G行列式求不出，画图找不到根
(*定义材料常数*)
Subscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
  M6 V$ D9 ?3 |' ^+ G7 X) Y$ xSubscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
3 M( u, j8 m; P/ x+ c1 rSubscript[\[Nu], LT] =
  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
Subscript[C, L] = Subscript[EE, L]/\[Mu];
& E$ k/ e5 a8 g# F) sSubscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
Subscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
(*脱层几何参数*)
\[Rho] = 1600;(*脱层密度*)
) \& X4 z! @  c0 S& oc = 0.06;(*阻尼*)
) Y2 ^4 D& q4 r, w% k. Z9 pt = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
l = 1;(*为梁总的长度*)
9 a$ k9 r8 k! w3 i) F! z& ~/ PSubscript[l, 1] = 0.3;(*为梁1区的长度*)
! F( x/ |* `  r" cSubscript[l, 2] = 0.3;(*为梁2区的长度*)
Subscript[l, 3] = 0.4;(*为梁3区的长度*)
! M/ X- `- j! ~" Z$ [6 ySubscript[l, 4] = 0.4;(*为梁4区的长度*)
1 o" A( a! }8 N) ]1 B' Y5 n3 lnn1 = h/t;(*nn1为总铺层厚度*)
nn2 = 2.0;(*为上子铺层层数*)
8 b  R& U) D0 E- Qnn3 = h/t - nn2;(*下子铺层层数*)
% P: C5 ~  M& p; oSubscript[h, 3] = t*nn2;(*梁的上部子厚度*)
& t+ |* R& s/ RSubscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
/ U7 b: s: Y+ M4 q/ N7 Y. q9 K0 nSubscript[h, 2] = h;(*梁2区的厚度*)
(*求各单层刚度*)
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)
' g) d) @# O5 S6 M' W
8 S- W) u9 v* l
ck = Cos[agk];(*列表求各铺层角度余弦*)
sk = Sin[agk];(*列表求各铺层角度正弦*)
s11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
' T/ l" h  B: [# {5 M7 \             sk^4*Subscript[C, T] +
   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
! E# ]' H2 J( j/ s- _% [A11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
  k1 a9 X. H. N    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
A31 = (\!\(
( L: ~2 \: Y1 Y/ S1 P( L4 B\*UnderoverscriptBox[\(\[Sum]\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
" }6 w5 J0 A7 i- OA41 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
% e, {( g5 [6 ?/ D, K, f% h     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
) t5 X" _, _& G; s

Print["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
D11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
1 z" @; {6 Z6 ]0 u% t3 d    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
9 C7 O4 C: i. |0 q\*FractionBox[\(t^3\), \(12\)])\)\)
5 m/ R, [! H* v0 m7 S, w, qD21 = D11
8 r0 F  n% B* Z' {2 {4 ID31 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
0 r' R+ a% z' ]" ]    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
" _; r1 s/ y& s. e: N; G4 ]\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
1 [5 q' z5 L" @: d9 i! Wk1 = 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)/(
% V; `% Q; ~5 q$ S  c DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
( B8 d4 V: b$ `, j9 \7 S 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 \
= h/Subscript[l, 3]; \[Beta]4 = h/Subscript[l, 4];
& R* U0 b) _( d- F7 B' ^b1 = Subscript[l, 1]/Subscript[l, 3]; b2 = Subscript[l, \
2 V9 b& J6 P) v3 I9 F8 |$ g1]/Subscript[l, 4]; b3 = Subscript[l, 2]/Subscript[l, 3]; b4 = \
1 w) m  u2 i1 A6 L; z5 Z' ]) USubscript[l, 2]/Subscript[l, 4];
DD1 = D11/(A11*h^2); DD2 = D21/(A11*h^2); DD3 = D31/(A11*h^2);
& j  m0 l5 q3 ?% e% g6 D3 d; h3 L) |& cDD4 = D41/(A11*h^2);
  S/ j' l% i  Rr11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
4 Z3 @$ m8 u+ W# I7 h a2*\[Omega]^2, (4)^-1]; r13 = Power[
& b. `! J* r; f9 |% l' G 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},
  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
   Sin[r12], 0, 0},
  {0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
   Cos[r12], 0, 0},
3 R) P  n( A( ], z! V) x  {-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
9 b$ T5 W) F2 w4 Y* a1 v7 x) H+ l( s   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
+ X8 j9 N# u2 U4 N- L! J   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
( Y  j2 G  S9 K   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
9 W( E2 h5 w( p& Z0 g3 f    DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
8 \- M' w( v( z, V" @. R- Y+ }  {-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
   0, -(r11)^3*DD1*(\[Beta]1)^3*Cosh[r11],
6 E2 g4 {; ?8 I6 g   0, (r13)^3*DD3*(\[Beta]3)^3, (r14)^3*DD4*(\[Beta]4)^3, -(r11)^3*
    DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
! P, m$ E0 ~2 @, H4 B   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
   0, -(r13)^3*DD3*(\[Beta]3)^3, -(r14)^3*DD4*(\[Beta]4)^3},
2 g  K( M* T  ^9 s1 Q  {0, -(r12)^2*DD2*(\[Beta]2)^2, (r13)^2*DD3*(\[Beta]3)^2*
7 Q: O3 c  D- V/ \0 v* H4 y    Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
  H3 F9 ^0 Y3 w# @2 j   0, 0, (r13)^2*DD3*(\[Beta]3)^2*Sinh[r13], (r14)^2*DD4*(\[Beta]4)^2*
  F' N0 Q7 l8 A. ^$ F6 o# m! j# E/ W* a    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*
' L& v( i/ D9 ~/ A4 V  N1 E    Sin[r14]},
" |* S8 N3 ]" @& F9 u  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
0 Q' u! n0 b$ k, G  u+ V, Q" r5 e& }& q    Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
- D4 p) ~/ Z# P: H) o/ \   (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]},
  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
+ X) o- ?/ J$ D) Q  J! R1 O- w   Sin[r11], 0, 0, 0},
! t# r$ h- S6 S- c5 s1 A" Q  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
: Q( W6 @. J/ z, W/ W1 G; K   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
9 Q" A  H2 p- a9 l   0, -r13*b1},
  {Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
5 z9 p. ^$ E5 {& {( d   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],
   0, 0, r12, -r13*b3*Cos[r13], 0},
& m+ \* O+ p2 V3 Y; j. p1 T) V  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
    0, -Sin[r14]},
  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
; t# m4 U& i% A' k: J) H' r% ]   r14*b4*Sin[r14],
   0, r12, 0, -r14*b4*Cos[r14]}}
7 E6 U: ]5 X7 H7 n6 ~* y6 C+ uaa1 = Det[BB]
Plot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

: w# _0 z& ]5 H9 u' d- t/ u' z$ @0 B

百年孤独

这是什么语言编程？

游龙游龙

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

mathematical编的