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

游龙游龙

麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
9 q1 V! `2 W  A# a, c0 ?行列式求不出，画图找不到根
(*定义材料常数*)
8 S, g0 y) I" \0 L" GSubscript[\[Nu], TL] = 0.25;
  L! \$ H9 y5 V) A2 @Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
8 _) ?" ?/ ^! g/ |/ w+ H- I1 \Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
* L2 w# {' [+ y/ ?% ^  a5 GSubscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
Subscript[\[Nu], LT] =
' b, n. j4 s. R1 x  g  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
$ o1 ], |% ~, ]; m4 {, fSubscript[C, L] = Subscript[EE, L]/\[Mu];
Subscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
3 ~/ W; V/ ?' iSubscript[C, T] = Subscript[EE, T]/\[Mu];
' v. ~* }/ x' f3 k# d8 ?, ^Subscript[C, LT] = Subscript[G, LT];
(*脱层几何参数*)
) _8 N: ]3 Y9 T! |\[Rho] = 1600;(*脱层密度*)
# R0 G, E9 Q) {  [- n5 Hc = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
2 G& ^& B& B1 c+ U" L9 bl = 1;(*为梁总的长度*)
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
, U: W' x0 x+ e7 W2 P  F9 x& v5 L2 FSubscript[l, 2] = 0.3;(*为梁2区的长度*)
Subscript[l, 3] = 0.4;(*为梁3区的长度*)
Subscript[l, 4] = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
: x) j0 W- h2 T4 c; |/ Qnn2 = 2.0;(*为上子铺层层数*)
nn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
Subscript[h, 2] = h;(*梁2区的厚度*)
(*求各单层刚度*)
  x& s7 v& a* U+ uagk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)
- ^+ U+ F" i" i, d0 S- I
7 Y6 a$ `) K8 q: @7 I8 J
0 e6 r% I" [  s1 T7 F0 ick = Cos[agk];(*列表求各铺层角度余弦*)
sk = Sin[agk];(*列表求各铺层角度正弦*)
s11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
             sk^4*Subscript[C, T] +
   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
  |: V/ A- B- eA11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
4 w( u2 ^6 {1 I+ ?0 ?: E9 Z    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
A31 = (\!\(
* D8 a5 N- i6 Y1 |7 I\*UnderoverscriptBox[\(\[Sum]\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
, S: W% q/ d7 M& n: _* YA41 = (\!\(
0 u- p6 M" e) s7 a5 `5 ?\*UnderoverscriptBox[\(\[Sum]\), \(k =
" [2 K3 v' L! ], T8 k$ ~     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
. _2 b0 p0 a; A& A- I

" P0 h' @* v6 a; k) l/ mPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
D11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
D21 = D11
/ T, C! O/ x4 i  e+ Q: {5 V% SD31 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
1 A* V- j5 ~; aD41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
6 B1 _5 Y9 C( d/ h8 D$ ?. G. X\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
k1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
; |1 m5 P7 A9 U' F$ P0 Y8 Ja1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
6 L. q0 W* t' Q% Z DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
; [+ r4 v- x) m. e- {  y; a$ U DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
' ?5 t% `9 Y" J\[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 = \
: g4 E- f# G) o; z8 v& D1 FSubscript[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);
r11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
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},
) ^2 b& {9 {, d! h  {0, 0, 0, 0, 1, 0, 0, 0, 0, 0, 0, 0, 1, 0, 0, 0},
5 n9 w* K, _) j/ D  {0, Cosh[r12], 0, 0, 0, Sinh[r12], 0, 0, 0, Cos[r12], 0, 0, 0,
- U) h& g% k" y: ?$ B  `   Sin[r12], 0, 0},
! w7 w% W' A/ ^6 E/ o4 d  {0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
   Cos[r12], 0, 0},
1 x. T5 H9 j4 ^  {-(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,
  {: v; f+ k3 O1 f; J% s   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},
* y4 h7 x% z# ^, z! X2 Q  {-(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,
2 O/ F0 U1 b) R+ F* r$ C3 I* l( }; `   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
5 C! R; l% g2 T9 Q  c; y1 i   0, -(r13)^3*DD3*(\[Beta]3)^3, -(r14)^3*DD4*(\[Beta]4)^3},
7 ^1 _2 s) I# ?2 H5 ^  {0, -(r12)^2*DD2*(\[Beta]2)^2, (r13)^2*DD3*(\[Beta]3)^2*
4 {6 q4 d) c7 `5 V! ^. @  f    Cosh[r13], (r14)^2*DD4*(\[Beta]4)^2*Cosh[r14],
, ]! G/ B! U# v: t3 k& J  I! [   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*
    Cos[r14], 0, 0,
   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
! Y' ]8 q- ]9 O9 F% G# n    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*
; k$ _& H' P  x; e$ q; M3 F' X    Cosh[r14], 0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sin[r13],
   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
9 ~! |! x0 k4 V" s% S   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
& e) @- @3 e& u1 O  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
2 ]8 B* ~+ C; x' i  ^, V% L   Sin[r11], 0, 0, 0},
+ t: ^# ~. o5 R% s/ C  l) x  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
: O7 K7 y+ g: ?# W* V  C  Z0 _) T) }/ g   0, -r13*b1},
. v, l) g( s2 J: m& n  {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,
* J/ I* Y/ x2 `# T, {   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
8 Z: s7 V% _; @! b+ r   0, 0, -r14*b2},
  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
% ~9 J, l# w' l: l    0, -Sin[r13]},
  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
   r13*b3*Sin[r13],
/ V4 u4 V* x+ V% o   0, 0, r12, -r13*b3*Cos[r13], 0},
. w) p8 Y- q2 c" Z' j$ q& h2 T2 B  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
    0, -Sin[r14]},
8 S8 i8 `, g1 k  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
+ x  J' D9 o' Y' @( ~5 W   r14*b4*Sin[r14],
   0, r12, 0, -r14*b4*Cos[r14]}}
aa1 = Det[BB]
Plot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

5 B! N' q( s8 M8 M* R; Z3 k

百年孤独

这是什么语言编程？

游龙游龙

百年孤独 发表于 2014-7-8 09:30
" k- _$ L: _0 w" E  g# d5 u这是什么语言编程？

3 ]- ^/ J7 t/ Tmathematical编的