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

游龙游龙

0 `2 }4 r' o& x* r, ]' {
/ w% z. ?: B7 g5 C' f麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
1 C! [/ ?7 S5 v. V/ @行列式求不出，画图找不到根
(*定义材料常数*)
Subscript[\[Nu], TL] = 0.25;
5 t/ \! ?% n, cSubscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
' m/ a$ ^* c3 L' _Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
. z6 f( r! N0 H/ @Subscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
Subscript[\[Nu], LT] =
+ R' t; K9 C( ?! D8 N% W& Y  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
0 `, [' {/ u( P9 r# |7 N" 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];
4 m0 P1 G) `+ |8 ~0 G9 n- {Subscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
' e0 D4 O1 A- x9 m(*脱层几何参数*)
\[Rho] = 1600;(*脱层密度*)
8 h' j, D* b/ C5 G; ?c = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
. W% ?( Z& {$ V3 _# N% c( S" fl = 1;(*为梁总的长度*)
! U! [8 G, B% H- P4 Z/ i1 cSubscript[l, 1] = 0.3;(*为梁1区的长度*)
Subscript[l, 2] = 0.3;(*为梁2区的长度*)
Subscript[l, 3] = 0.4;(*为梁3区的长度*)
4 H+ z8 A6 E& `: j( pSubscript[l, 4] = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
: ^0 [5 n2 b6 n5 tnn2 = 2.0;(*为上子铺层层数*)
, f& r; Q8 p, }; G( C+ X( |/ W0 Inn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
" c6 w5 ^/ l- {! o$ v9 Z' fSubscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
" Z" x& [* u) B) E2 j9 TSubscript[h, 2] = h;(*梁2区的厚度*)
(*求各单层刚度*)
# G+ p$ \5 H7 G% [6 g4 l( hagk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)
/ n4 k7 z, W2 C7 H. J( \
" K; T! y2 \- V: z1 |' c: J
ck = 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];(*求各区刚度*)
$ F9 u% l- {' CA11 = \!\(
6 Z( H& I4 _  {2 `\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
A21 = A11
- O% m; X% O. C3 Q# VA31 = (\!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
A41 = (\!\(
4 [+ N* Z# {5 |1 |$ F\*UnderoverscriptBox[\(\[Sum]\), \(k =
$ v& {& a7 J; L/ F     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
  K9 q9 \6 i: Z3 f3 U
* x" E- P4 z; B- g1 t4 N
  x4 ]0 h; y3 s2 g; T; cPrint["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
& v4 C: M$ c9 X. C$ Y1 o- E1 yD11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
0 M0 \& f) I  ^0 w0 V    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
7 c" I' ?$ J6 s7 g2 H. s' H$ l7 }\*FractionBox[\(t^3\), \(12\)])\)\)
( j+ Z2 a* Z1 ]9 {' @3 Y, ?+ vD21 = D11
D31 = \!\(
& O3 k  \( {( Z1 R* e+ E* A) l7 A0 p\*UnderoverscriptBox[\(\[Sum]\), \(k =
1 b' U5 R) i$ I& M0 G9 {' |: c    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
3 K$ M  T+ m# E- S) f3 i* j/ K\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
. j6 N  F# p- f3 |1 w" m4 Z' \1 z6 c+ XD41 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
2 M: }8 A/ w9 |8 H! I; F% Q% Z    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
- \+ N% P& H+ h3 o8 R5 \, H\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
k1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
) e/ G, Z& b' X+ b9 j. s0 g3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
* g2 l5 B' R& B- `! ca1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
, W& w. ?8 w: x. Q DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
1 a0 A- ~/ ^( ~. X9 L& z' E DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
0 ]5 F" I. `" }; [9 `: i\[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 = \
4 r# D/ A2 q; D' V/ q0 TSubscript[l, 2]/Subscript[l, 4];
/ Y. S; o# A- C" XDD1 = 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[
! E2 _8 u* V1 O; k9 l 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,
   Cos[r12], 0, 0},
- v) {- J# P0 V& v" ~! Q7 D  {-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
8 n- S& Y3 @: x& B/ n  e4 N   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
7 D2 x, C* n6 e4 ^+ \5 M   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
    DD1*(\[Beta]1)^2*Sin[r11], 0, 0, 0},
; _# t2 K/ j: c2 k9 k+ A8 W' U2 {  {-(r11)^3*DD1*(\[Beta]1)^3*Sinh[r11], 0, 0,
7 v$ G- V0 C" a9 N. ?; ^   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,
# j6 Z6 k0 z7 b0 f   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
2 ~. P5 h: w! r: S6 B8 N' g; N   0, -(r13)^3*DD3*(\[Beta]3)^3, -(r14)^3*DD4*(\[Beta]4)^3},
, |3 F% j# c% p5 N. v6 y  {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*
6 v0 m" \% s7 h    Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
8 J4 F- }3 x$ d% R2 Z- l   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
, V) v: a4 l" D    Cos[r14], 0, 0,
   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
& G0 n9 X- h4 n' f    Sin[r14]},
( D3 F+ ^8 L5 n5 n/ O% x" h9 H  {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,
& n% o, A0 _  @; e$ h# Q   (r13)^3*DD3*(\[Beta]3)^3*Cosh[r13], (r14)^3*DD4*(\[Beta]4)^3*
4 p/ x) F, H$ a3 L    Cosh[r14], 0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sin[r13],
9 V0 K& N5 P: U3 N   (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]},
  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
+ r1 \( y2 W  Q7 b   Sin[r11], 0, 0, 0},
; Q) j/ b5 x0 b. S  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
- e' r; i9 g( n! D   0, -r13*b1},
  {Cosh[r11], 0, 0, -1, Sinh[r11], 0, 0, 0, Cos[r11], 0, 0, -1,
. F1 ^: R! c& s- l) z   Sin[r11], 0, 0, 0},
% ^0 ?9 A2 n) K  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0,
% p7 E. S0 X5 _2 U) g; v" _   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
- i0 j8 p2 r9 x   0, 0, -r14*b2},
  {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 c; z% d) Q4 O1 |, x   0, 0, r12, -r13*b3*Cos[r13], 0},
+ T' R* @4 x9 ^# ^0 e+ N  {0, 1, 0, -Cosh[r14], 0, 0, 0, -Sinh[r14], 0, 1, 0, -Cos[r14], 0, 0,
    0, -Sin[r14]},
  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
! a0 c, h: Y9 j' a8 t   r14*b4*Sin[r14],
4 C& d4 ]; P3 Y' H- j$ n; T9 s# s   0, r12, 0, -r14*b4*Cos[r14]}}
aa1 = Det[BB]
9 T8 i( }3 J9 n! _6 V8 _, ?Plot[aa1, {\[Omega], -10, 10}]
ss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]

$ Y! o; X9 I2 Z1 j

百年孤独

这是什么语言编程？

游龙游龙

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

mathematical编的