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

游龙游龙

5 r/ {5 M% \* T/ M) O
麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
* P. ~1 n% ?$ Z5 K行列式求不出，画图找不到根
(*定义材料常数*)
$ ?$ v$ T7 E' J, [* qSubscript[\[Nu], TL] = 0.25;
0 T  N( V7 p6 s9 CSubscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
- ]& \, G* ]) D$ w# I4 sSubscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
1 s2 J9 [5 n% A5 R$ h' SSubscript[\[Nu], LT] =
/ g0 T0 T0 Q' Y/ R5 I# N0 m  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
+ E+ x% B1 q* {- n- ^\[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];
3 W2 F2 p, p+ gSubscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
5 ~  E$ Z$ V5 d* r: f3 r(*脱层几何参数*)
\[Rho] = 1600;(*脱层密度*)
/ C' Q4 c/ I4 ^5 [, W8 J7 t( pc = 0.06;(*阻尼*)
  g0 T+ v( C7 H! M' j6 G2 {% Ot = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
2 j1 K" y2 _4 s6 U$ p; o/ fl = 1;(*为梁总的长度*)
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
- v- @" x4 @& c+ v- {$ OSubscript[l, 2] = 0.3;(*为梁2区的长度*)
1 A% O: n( {# ]2 Z$ K. PSubscript[l, 3] = 0.4;(*为梁3区的长度*)
) a6 [  v. {- K$ y3 R- hSubscript[l, 4] = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
) G3 y4 v8 o. L& C) a- Jnn2 = 2.0;(*为上子铺层层数*)
nn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
( o9 l+ F. f; Z; P8 k+ F3 M) eSubscript[h, 1] = h;(*梁1区的厚度*)
* z7 M3 _; h- S4 rSubscript[h, 2] = h;(*梁2区的厚度*)
: Z# V/ S' H5 w' C(*求各单层刚度*)
  v) D9 \. ~! U+ a. P; ~$ U  v7 ]. ?8 [agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)
; d- B! |# I: H! W, i, n2 \
: m8 T$ U* E0 N
; ]' v( S$ i4 c0 vck = Cos[agk];(*列表求各铺层角度余弦*)
sk = Sin[agk];(*列表求各铺层角度正弦*)
8 v" V8 N6 H0 q, ls11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
             sk^4*Subscript[C, T] +
. n0 @- K. w  T0 f# P  ~5 h; Z9 T   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
A11 = \!\(
( Z2 _! y' y6 n2 V+ I" m\*UnderoverscriptBox[\(\[Sum]\), \(k =
: k$ F% E7 x  x" r    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
* b( b4 Z( n; d; f; VA31 = (\!\(
( H! \" T! n' Y7 Y9 ?7 w\*UnderoverscriptBox[\(\[Sum]\), \(k =
, s1 O) K+ S4 |/ m, `     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
A41 = (\!\(
/ I" D& h' k8 j( g! {\*UnderoverscriptBox[\(\[Sum]\), \(k =
* E% r- g) V- C0 w/ O! C  A     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
! H9 p& h( `" N4 }
! H. {+ D1 j; o$ F1 I* v" @
, y" y9 ^4 h) T$ t) q8 P6 R! J' xPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
D11 = \!\(
5 ^3 L. z* T7 j2 \  h6 p3 f\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
0 z0 P2 l, h6 [7 V% K! ?# t\*FractionBox[\(t^3\), \(12\)])\)\)
" g0 }; \5 A! P8 z4 F/ ~2 b8 x. uD21 = D11
D31 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
( {/ E* W, o4 T3 c* p; ^\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
5 x- T. K6 L. I# G\*UnderoverscriptBox[\(\[Sum]\), \(k =
' o" W3 c: y$ g! q1 E; s# |* V  @1 A) b    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
* {5 _" i, K4 v. C; h8 D\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
" K6 ~  @+ S1 Bk1 = 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)/(
1 _. ]; N  q+ i/ ]6 e+ Y' W7 ` DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
! ~/ B5 K) j; t* j DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
" E( G/ S, ?8 V; C4 j7 B( U9 q& x+ v" q\[Beta]1 = h/Subscript[l, 1]; \[Beta]2 = h/Subscript[l, 2]; \[Beta]3 \
= 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];
DD1 = D11/(A11*h^2); DD2 = D21/(A11*h^2); DD3 = D31/(A11*h^2);
1 c5 U; h! C9 R- d; Z5 QDD4 = D41/(A11*h^2);
r11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
8 C2 W5 `+ v" A! S2 x$ r" R  G a2*\[Omega]^2, (4)^-1]; r13 = Power[
9 X: _0 F, [  h3 ~ 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,
1 m7 |! ^7 F# g& k' \9 i2 e8 X/ T2 f   Cos[r12], 0, 0},
4 |% G! I/ y# b+ y  \/ u  K5 f/ 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],
" q0 m+ O2 ]2 U! o- S4 f" {4 Z   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
    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],
% V: R  T6 ]" A* C3 L   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*
' @0 T1 z( P! x( i' S9 ]    Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
% S0 U( s& |2 I  V    Cos[r14], 0, 0,
   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
. L+ x! \1 H9 E: {7 S    Sin[r14]},
  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
0 _5 G2 n) J/ G" D* f; k9 f    Sinh[r14], 0, -(r12)^3*DD2*(\[Beta]2)^3,
   (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],
; [; a  b! }9 Z   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
  {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,
   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,
% k3 q. }% \* T  w' l) t   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],
6 t4 P7 ^2 O( i" E* a5 z& _  Y   0, 0, -r14*b2},
  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
* v- i: Y8 [; p0 J' p# w; _5 a    0, -Sin[r13]},
  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
! B/ h3 ?- m# G   r13*b3*Sin[r13],
   0, 0, r12, -r13*b3*Cos[r13], 0},
. ]% E# v" H4 x4 R; n8 q  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
" j; a; l) M& @! y. N+ B    0, -Sin[r14]},
. y7 N) ]. j9 m) I2 M  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
) l! G% u0 ?& i! h4 P   r14*b4*Sin[r14],
5 B7 Z% G( Y- W   0, r12, 0, -r14*b4*Cos[r14]}}
/ x( [$ f% G, x7 f& g* @0 Kaa1 = Det[BB]
( @; E5 x/ O) a( }) {Plot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]
0 n% b, @- l( X  e; d, s2 g

百年孤独

这是什么语言编程？
) R# W5 i3 K1 p2 i3 O8 e

游龙游龙

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

mathematical编的