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

游龙游龙

  [! x0 ]2 @# X1 d* v
麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
行列式求不出，画图找不到根
(*定义材料常数*)
( H& V$ I1 s, Z$ Z- y0 M0 {Subscript[\[Nu], TL] = 0.25;
0 }* {0 L7 f, }0 LSubscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
7 Z. Y+ L: l: o* uSubscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
' h% H4 n# e9 c5 XSubscript[\[Nu], LT] =
: b8 v- K2 U& r  P# p  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
\[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];
8 o2 c% N! B/ q. k' cSubscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
(*脱层几何参数*)
\[Rho] = 1600;(*脱层密度*)
: B& a: o: ~) W3 V) S3 Oc = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
7 ^1 _2 o* _2 O7 }. X8 b0 Wl = 1;(*为梁总的长度*)
Subscript[l, 1] = 0.3;(*为梁1区的长度*)
Subscript[l, 2] = 0.3;(*为梁2区的长度*)
Subscript[l, 3] = 0.4;(*为梁3区的长度*)
6 h- b7 g/ e, G5 kSubscript[l, 4] = 0.4;(*为梁4区的长度*)
( c4 P) k- I5 ?2 c8 p9 S- `% mnn1 = h/t;(*nn1为总铺层厚度*)
nn2 = 2.0;(*为上子铺层层数*)
nn3 = h/t - nn2;(*下子铺层层数*)
$ e* k. v# f, r: ^; T+ L( J) mSubscript[h, 3] = t*nn2;(*梁的上部子厚度*)
Subscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
4 B9 h1 u1 V$ L5 x, v) tSubscript[h, 2] = h;(*梁2区的厚度*)
(*求各单层刚度*)
6 S1 t" e& A) a- F, S: P4 \1 xagk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)


6 l: Y5 A/ ]; y( X2 Bck = Cos[agk];(*列表求各铺层角度余弦*)
* _' P( M+ w5 A: M; O% psk = Sin[agk];(*列表求各铺层角度正弦*)
0 d% Z8 |4 m8 |% b7 J  ?" As11k = ck^4*Subscript[C, L] + 2 ck^2*sk^2*Subscript[C, S] +
             sk^4*Subscript[C, T] +
, g1 W% A% I8 Q$ g8 I3 G# G5 e   4 ck^2*sk^2*Subscript[C, LT];(*求各区刚度*)
A11 = \!\(
3 n6 m  d0 P& m4 p4 P/ a; G\*UnderoverscriptBox[\(\[Sum]\), \(k =
* D/ d) V: ?/ Z' v    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
& S+ I) j4 C8 p# c4 e' uA31 = (\!\(
2 J) C; ?. W2 N- J' |% q\*UnderoverscriptBox[\(\[Sum]\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
. F4 q" [( W3 n1 V. `9 aA41 = (\!\(
- }; A  ]8 o9 x' J- T3 v\*UnderoverscriptBox[\(\[Sum]\), \(k =
: d, [5 ?. U/ d4 S4 F* @     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
8 `- Q2 D2 _2 G. P' N4 t

Print["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
+ E3 e" \% F' L3 S7 mD11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
3 p6 U& J+ O: f0 p2 j$ ED21 = D11
1 E2 t7 D. Z; \0 F8 m, nD31 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
& \) }0 h: @% P1 P    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
. V7 {5 k0 t8 `" m9 T: R\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
D41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
! l; D9 U4 W% Z7 p    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*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;
$ `# r* |% I0 u& _8 d; ~a1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
) \5 |% v3 _: m+ X1 E  P% W DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
. v4 M/ F  A& t) w; a& t, 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 = \
/ g' M5 E2 t5 b1 n* ySubscript[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[
" s& i. ~# _  {& T$ H a3*\[Omega]^2, (4)^-1]; r14 = Power[a4*\[Omega]^2, (4)^-1];
1 m. Z+ H: Y* g, ]  F: W) N+ HBB = {{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},
3 h6 R* _. P* o6 V  {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},
  {-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
, x7 Q$ z5 X: ]- j   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
4 Q; _6 e. R% D   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},
$ n% W$ P1 D) _; Q% D  {-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
+ [  h4 `' a; @   0, -(r11)^3*DD1*(\[Beta]1)^3*Cosh[r11],
+ ]# R- a: ?! \  c   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 @/ `+ K2 ~9 f) A: c   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
  Z# A, N7 P, {! ]9 j7 X# d3 K   0, -(r13)^3*DD3*(\[Beta]3)^3, -(r14)^3*DD4*(\[Beta]4)^3},
. k- l0 o" O! ]' ]  n  {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*
    Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
% S7 t. w0 V0 B( I$ S   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
9 v6 s# z1 ~  u0 [6 `' o- m    Cos[r14], 0, 0,
   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
    Sin[r14]},
' \: m0 R3 G0 U0 i- d& Y2 ^/ g" \  {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,
3 S% e# @8 Q, b0 N$ F2 i7 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],
) ~; z4 P% u; {( Y, S   (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]},
  L% K% A. s8 C& W. [& [  T  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
  I) C* e5 S5 O! i: F   Sin[r11], 0, 0, 0},
  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
! a; V5 y- A/ F( e( ~   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
! v- J0 k8 \- _$ |8 u8 |: g  s9 r+ Q   0, -r13*b1},
  {Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
6 G! H7 [% U. [$ ?4 R! z: x   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],
+ b% e( A& }" g9 w+ `5 X9 I   0, 0, -r14*b2},
% B# o8 a3 O. p$ V& ^5 q: k  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
    0, -Sin[r13]},
* V: J* ]1 W- d& {  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
' a' f0 y9 ^  ?# {   r13*b3*Sin[r13],
% L$ p% t* [# c: ~6 f! a" S   0, 0, r12, -r13*b3*Cos[r13], 0},
6 B4 r1 c: v- B4 F, e; h! G$ y  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
/ G7 F, u1 }! S, W    0, -Sin[r14]},
  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
% ]! l) m$ U; U% s   r14*b4*Sin[r14],
/ j( Z) _# ]: S0 U2 v# [! ]   0, r12, 0, -r14*b4*Cos[r14]}}
$ Z2 b( K, Y& ~0 j* @aa1 = Det[BB]
4 G1 ?0 x7 c, {/ RPlot[aa1, {\[Omega], -10, 10}]
1 D# r+ F2 Y: Q7 nss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]
9 A% ]% [; \+ O- ]% d0 Q( Z% l: n

百年孤独

这是什么语言编程？

游龙游龙

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

5 V+ O1 I7 m- C. Qmathematical编的