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

游龙游龙

- F" S2 ~/ K/ {6 a; G
2 F4 Q7 x, Q8 x2 g/ \  J麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
" X6 c! q& ]7 U- z9 L行列式求不出，画图找不到根
9 ?1 w5 L- N/ _3 g3 d# h2 H4 t(*定义材料常数*)
5 ]; U& F) }) n# T7 aSubscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
( Q9 Y# T% h$ l; l; `+ R. \Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
7 I( i3 O% K2 y% C/ TSubscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
9 O0 C2 c! q) V+ I  B4 tSubscript[\[Nu], LT] =
  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
. F$ b/ |" N2 }7 N7 e. u# n; A/ @\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
Subscript[C, L] = Subscript[EE, L]/\[Mu];
2 A  o5 C( g( i  ^" tSubscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
) I2 J9 m) Z) X9 B) qSubscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
: @1 t& C" \4 M% j- \' q( y(*脱层几何参数*)
\[Rho] = 1600;(*脱层密度*)
c = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
l = 1;(*为梁总的长度*)
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
Subscript[l, 2] = 0.3;(*为梁2区的长度*)
/ W) {1 b$ P. k2 S$ z/ dSubscript[l, 3] = 0.4;(*为梁3区的长度*)
Subscript[l, 4] = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
" w/ H* ]1 D  i, U: R- a, mnn2 = 2.0;(*为上子铺层层数*)
nn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
: s+ S6 {& d' o9 jSubscript[h, 1] = h;(*梁1区的厚度*)
Subscript[h, 2] = h;(*梁2区的厚度*)
' i! q5 j; l3 d$ v' X3 w(*求各单层刚度*)
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)

  O! g7 @" n( N( ?" h
ck = Cos[agk];(*列表求各铺层角度余弦*)
sk = Sin[agk];(*列表求各铺层角度正弦*)
s11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
1 }/ {- K) ^, T  P+ u/ r9 p* y             sk^4*Subscript[C, T] +
. ]4 g7 p7 M: |% [6 ]* i. F1 x   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
/ U& E8 X5 W. y0 tA11 = \!\(
2 ?7 V/ \9 G6 j* y3 U( _3 e, J\*UnderoverscriptBox[\(\[Sum]\), \(k =
9 U& c* A! F. d- o* {3 o# H    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
A31 = (\!\(
0 z* a& C5 }4 u" M4 J' F\*UnderoverscriptBox[\(\[Sum]\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
" h* V1 Y& U8 @5 Y. t" J# OA41 = (\!\(
: S  |6 ?* Z. m/ e& v\*UnderoverscriptBox[\(\[Sum]\), \(k =
. j# e) P) R, G; D. q3 p# r     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)

& I. p, w* M# j! M9 L9 p2 b! i* y
) \" \& H2 r. p. v4 E5 qPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
D11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
) L6 V  t, ^; ?) Q1 y' l+ ^    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
: N1 t9 O6 L9 |. B6 MD21 = D11
D31 = \!\(
, I3 N* Y0 J/ F& ^* F: E\*UnderoverscriptBox[\(\[Sum]\), \(k =
7 {; T8 W2 e1 H: L: D' b( d, ^    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
4 d5 d( P! u. lD41 = \!\(
! Y) ^- r8 J) S# M: v\*UnderoverscriptBox[\(\[Sum]\), \(k =
9 H( E( a& j# U1 z- ^    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
& V) b4 W0 |9 h& l0 hk1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
0 j( b( G3 K$ M2 d: `6 @3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
. J" {7 i8 n6 x; H8 \; n# ca1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
* Y1 b# ]9 ?( H$ Q DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
% b" c$ ]' Y& t; }& g. w. a, x DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
" M0 q, w4 \, z/ j2 r2 N\[Beta]1 = h/Subscript[l, 1]; \[Beta]2 = h/Subscript[l, 2]; \[Beta]3 \
- O) m5 y% N+ k1 F= h/Subscript[l, 3]; \[Beta]4 = h/Subscript[l, 4];
* E7 i! Q1 ]8 \b1 = Subscript[l, 1]/Subscript[l, 3]; b2 = Subscript[l, \
1]/Subscript[l, 4]; b3 = Subscript[l, 2]/Subscript[l, 3]; b4 = \
: l" d' \0 X% }& [* t% B5 Y1 T9 jSubscript[l, 2]/Subscript[l, 4];
DD1 = D11/(A11*h^2); DD2 = D21/(A11*h^2); DD3 = D31/(A11*h^2);
" e" H: ^% L. S$ d8 o! b' n1 @DD4 = D41/(A11*h^2);
# {2 F, Q( a; k7 I5 f$ W% G9 wr11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
a2*\[Omega]^2, (4)^-1]; r13 = Power[
. l) B# C. L) a+ o8 Y 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},
  {7 Y5 j- b3 }! @. y* \  {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,
6 V% K$ S; A, |   Sin[r12], 0, 0},
" ^: L% r- M9 X1 s$ t) d& m  {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],
/ M4 q3 l& ?4 d& r  T   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
. e* n; S3 ^' O: }! V* T   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
& @% s2 L# W8 u. Z/ y7 T4 B5 V    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*
    DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
- B3 p0 B: u  a  w0 z% [   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*
% O& Z3 s/ z( w& |; d1 S    Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
, j2 A: Q" m, R7 O3 l   0, 0, (r13)^2*DD3*(\[Beta]3)^2*Sinh[r13], (r14)^2*DD4*(\[Beta]4)^2*
- S+ `+ r" {' Z9 g' N    Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
& t7 B" j; S" d2 g: p5 y( W7 w! z' h   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
    Cos[r14], 0, 0,
8 P1 N7 @# X" V; y8 t   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
    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,
/ ^$ o6 j6 j0 z! Q1 j   (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],
/ _7 G& ]( S* L: _   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
; M! K* ]2 M5 h6 B; M+ d   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
$ g5 E. y6 X4 N6 U  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
2 u1 Y5 P2 x: B- ^/ y   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
  i2 R& N/ o( P; i6 w0 B   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
" v# @5 `! v0 C1 M+ i9 A; f   0, -r13*b1},
  {Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
$ `* g5 r4 ~/ E) ?+ e, b. e4 D   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],
   0, 0, -r14*b2},
  e, O0 c' u' M9 X( f  {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],
1 c9 O4 }+ @6 j! a6 ^   0, 0, r12, -r13*b3*Cos[r13], 0},
7 v: G7 H( F: t0 q  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
$ a1 R! F9 W( x7 L    0, -Sin[r14]},
  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
7 _, c" |1 e7 _. ?, {' u& [   r14*b4*Sin[r14],
   0, r12, 0, -r14*b4*Cos[r14]}}
- Q- Z3 T, i8 Q& Baa1 = Det[BB]
Plot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]
" t+ X9 L" L, `- J, B0 ~: q. u

游龙游龙

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

- S4 ]+ S1 u5 R; m7 P( c9 ^mathematical编的

百年孤独

这是什么语言编程？