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

游龙游龙

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
行列式求不出，画图找不到根
(*定义材料常数*)
Subscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
, q9 a: @  }+ N9 h, Z# BSubscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
% v- u) Q7 G& s, m0 A/ ^Subscript[\[Nu], LT] =
  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
9 W6 C0 {( f" b( i6 K1 s% W\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
; D9 \9 V" a# o* [" O" m" L( `3 [Subscript[C, L] = Subscript[EE, L]/\[Mu];
Subscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
- \. Q$ }4 k. n. d& cSubscript[C, T] = Subscript[EE, T]/\[Mu];
5 Z, X# x/ i! q" ?Subscript[C, LT] = Subscript[G, LT];
(*脱层几何参数*)
3 N. C% F3 ?0 K' s; H\[Rho] = 1600;(*脱层密度*)
  Q& e8 E. Z  cc = 0.06;(*阻尼*)
& p. k! B4 l# G8 h8 xt = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
# X! ]. `, A2 p( a" Z4 p. bl = 1;(*为梁总的长度*)
1 w) m, o* W; U  \+ u* e: {Subscript[l, 1] = 0.3;(*为梁1区的长度*)
& T% x6 w0 ]- BSubscript[l, 2] = 0.3;(*为梁2区的长度*)
6 K4 F1 _5 [# n3 USubscript[l, 3] = 0.4;(*为梁3区的长度*)
Subscript[l, 4] = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
nn2 = 2.0;(*为上子铺层层数*)
nn3 = h/t - nn2;(*下子铺层层数*)
. r9 G5 R4 i  K1 z. K8 gSubscript[h, 3] = t*nn2;(*梁的上部子厚度*)
Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
4 y& T2 T  P& |" ZSubscript[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};(*列表各铺层角度*)
2 _& O1 a6 Z/ g) \( E* R( a
% z7 B$ K! f8 z$ [. U2 u
ck = Cos[agk];(*列表求各铺层角度余弦*)
9 {* H: U' c# L5 H: O* qsk = Sin[agk];(*列表求各铺层角度正弦*)
5 Y1 r/ R7 f4 v! ss11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
, B9 P8 K& M% c% D* X             sk^4*Subscript[C, T] +
   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
A11 = \!\(
  V) I8 U* P- ~" n, s\*UnderoverscriptBox[\(\[Sum]\), \(k =
6 ^2 M$ _7 Z/ ?5 b* t' [    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
A31 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
6 K1 I* E1 T: W, m% o& C$ O     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
9 w) i7 `8 R4 H9 w4 C& QA41 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
3 O8 O! [; v/ U, F9 v# O     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)

$ R8 n6 j' e7 {, m3 K2 m  O
& P8 i) U( H# Z8 z# RPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
3 p1 {% t9 Z1 lD11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
D21 = D11
D31 = \!\(
4 [% P' ]( q3 L$ Q3 ^\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
" n7 z/ t# X. O+ L& ~) ]\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
2 v7 X; \2 j/ ]1 k\*UnderoverscriptBox[\(\[Sum]\), \(k =
8 b* j' r9 H& ~; P- I  n    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
k1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
  h- O! P: z% g; m3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
& m' _5 X) i3 m8 e# U: Q" C7 xa1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
+ H2 j; x: N! n: m  Y& u) o DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
4 }- J. w) A& t. n- D\[Beta]1 = h/Subscript[l, 1]; \[Beta]2 = h/Subscript[l, 2]; \[Beta]3 \
& E( ~" b/ m  h( _, L$ e* h3 c= h/Subscript[l, 3]; \[Beta]4 = h/Subscript[l, 4];
) Y; P# X% e/ _" l, ]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];
% n, u$ P- ?6 n' l/ m' RDD1 = 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[
: H9 y* u+ ~+ w$ k! W a2*\[Omega]^2, (4)^-1]; r13 = Power[
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},
, g& r0 X) Q) ~) F3 f  {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],
   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
+ j! z6 E1 |' e4 C' y   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},
0 @9 c+ h& [- l0 u% ?1 d  I  {-(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*
/ r* P. J3 z7 a; n    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},
, ^/ A0 K$ \/ @8 t& x! V0 _  {0, -(r12)^2*DD2*(\[Beta]2)^2, (r13)^2*DD3*(\[Beta]3)^2*
    Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
# A3 Q" m4 d7 S# w" J6 \0 Q% R6 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,
   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
8 w. t5 @: G, L    Cos[r14], 0, 0,
% k9 s1 c, ~+ b3 A: {1 Z   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
2 B6 V9 k4 h3 w$ |* e    Sin[r14]},
  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
    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],
8 C: u( x* }& d, F: J) a# C$ q+ a   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
/ W0 Z! W! Z! k, c4 X% B, C9 V" `   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,
' a2 C$ h' E) O2 G   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
. I8 p0 @: H: Y) m2 `   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,
3 z6 ^% M! q0 p1 D& |" j6 D   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
$ v/ x/ Z" y& Z' a' y   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
  \& M& K( O3 J6 l   0, 0, -r14*b2},
  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
    0, -Sin[r13]},
9 f* u  C, {- M  T) }7 m  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
6 l( Q1 @7 N6 z# Q. x6 D9 Q   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,
    0, -Sin[r14]},
) N; b6 s" x  x  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
- x- M! N( X, q, s! J" a   r14*b4*Sin[r14],
   0, r12, 0, -r14*b4*Cos[r14]}}
5 y' W' F5 N" u- D( naa1 = Det[BB]
Plot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

& ?9 a1 R% E7 f: W

百年孤独

这是什么语言编程？

游龙游龙

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

mathematical编的