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

游龙游龙

' c  X& C, o% [7 m, j& o7 n
- J) H' m) L1 ]- ?4 H4 W6 B麻烦大家帮忙看下这个程序，不知道哪里出错了，小弟是刚入门，请高手指教
. |2 D2 j9 d( _8 e" e8 a9 x9 U行列式求不出，画图找不到根
& [' Q1 B3 Y: x8 B' f9 c  y! r  t6 q(*定义材料常数*)
Subscript[\[Nu], TL] = 0.25;
Subscript[EE, L] = 53.8*10^9;(*L方向的弹性模量*)
Subscript[EE, T] = 17.93*10^9;(*T方向的弹性模量*)
0 J( h* C9 x8 ^& O6 a% YSubscript[G, LT] = 8.96*10^9;(*LT方向的剪切模量*)
Subscript[\[Nu], LT] =
( M: l# ?7 |/ L, b  Subscript[\[Nu], TL]*Subscript[EE, L]/Subscript[EE, T];
+ x/ ]! K  g* n8 m; p\[Mu] = 1 - Subscript[\[Nu], LT]*Subscript[\[Nu], TL];
- H/ S9 l" x" v6 R0 D$ r- E) C. USubscript[C, L] = Subscript[EE, L]/\[Mu];
Subscript[C, S] = Subscript[\[Nu], LT]*Subscript[EE, T]/\[Mu];
Subscript[C, T] = Subscript[EE, T]/\[Mu];
Subscript[C, LT] = Subscript[G, LT];
0 t' G; R7 l/ C! B4 t5 i# F* d(*脱层几何参数*)
\[Rho] = 1600;(*脱层密度*)
0 v6 D; |+ E& @  Y4 I0 M# u7 ]c = 0.06;(*阻尼*)
t = 0.01;(*t每层厚度*)
h = 0.05;(*h为梁总厚度*)
l = 1;(*为梁总的长度*)
, u' Z2 r8 C! K' V+ \3 ^1 MSubscript[l, 1] = 0.3;(*为梁1区的长度*)
Subscript[l, 2] = 0.3;(*为梁2区的长度*)
6 T* H& a7 I9 F; zSubscript[l, 3] = 0.4;(*为梁3区的长度*)
Subscript[l, 4] = 0.4;(*为梁4区的长度*)
nn1 = h/t;(*nn1为总铺层厚度*)
1 w$ T9 j( G7 ^; C& tnn2 = 2.0;(*为上子铺层层数*)
! w8 h8 M) F  enn3 = h/t - nn2;(*下子铺层层数*)
Subscript[h, 3] = t*nn2;(*梁的上部子厚度*)
" ]" `0 O6 p& ^2 rSubscript[h, 4] = h - Subscript[h, 3];(*梁的下部子厚度*)
Subscript[h, 1] = h;(*梁1区的厚度*)
Subscript[h, 2] = h;(*梁2区的厚度*)
9 H7 s! X; S: n* q(*求各单层刚度*)
agk = {0 Pi/180, 90 Pi/180, 0 Pi/180, 90 Pi/180, 0 Pi/180};(*列表各铺层角度*)
5 M; P: }$ a- Q, K, h0 @" m& i" e
' H& `5 Q- |' ~0 t; y) h
8 b' O3 {' P" y& g1 f. pck = Cos[agk];(*列表求各铺层角度余弦*)
& S! `% P& u( X" @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];(*求各区刚度*)
" [$ c% a3 h) `2 \" n0 n5 vA11 = \!\(
5 ~' o: G) x4 o3 G( Z; `# t\*UnderoverscriptBox[\(\[Sum]\), \(k =
$ [, E2 [+ r8 B& N    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*t)\)\)(*一区拉伸刚度*)
5 |- O' O4 n2 v9 `! yA21 = A11
; V- N8 K7 B& C/ d# IA31 = (\!\(
3 ^  z0 L! A, f1 @9 Y\*UnderoverscriptBox[\(\[Sum]\), \(k =
! Z( ?$ i* R( ?+ @6 T  `& `     1\), \(nn2\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*三区拉伸刚度*)
A41 = (\!\(
! ~* l4 o) J% i! J8 u\*UnderoverscriptBox[\(\[Sum]\), \(k =
7 ^5 t/ y. s+ P. P# m+ s) s$ \     nn2 + 1\), \(nn1\)]\(\((s11k[\([\)\(k\)\(]\)])\)*t\)\))(*四区拉伸刚度*)
" @6 o. V6 `/ k/ _& D; b- J
" G2 [1 Y. p7 U) z& E6 c( J! ^, x1 j
  R* y5 `: b2 g+ |2 [Print["A21=A11;D21=D11"];(*二区拉伸刚度、弯曲刚度*)
; {7 b- a5 L( r  B# U  r& |D11 = \!\(
\*UnderoverscriptBox[\(\[Sum]\), \(k =
& }& Q: x2 W+ f7 D0 {& M    1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)
7 h" h' \2 o- X6 ID21 = D11
D31 = \!\(
2 R$ U) ?+ V& T1 P/ L) \\*UnderoverscriptBox[\(\[Sum]\), \(k =
    1\), \(nn2\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*三区弯曲刚度*)
, `2 d& [; P6 l. t  i4 FD41 = \!\(
5 i2 l! P) c( f( M\*UnderoverscriptBox[\(\[Sum]\), \(k =
    nn2 + 1\), \(nn1\)]\((s11k[\([\)\(k\)\(]\)]*
\*FractionBox[\(t^3\), \(12\)])\)\)(*四区弯曲刚度*)
& Z7 `8 c# V, ~6 U% q4 qk1 = Subscript[l, 1]/l; k2 = Subscript[l, 2]/l; k3 = Subscript[l, \
3]/l; k4 = Subscript[l, 4]/l; \[Alpha] = h/l;
% y2 x, k/ J6 j) H& }a1 = (\[Alpha]*k1)/(DD1*(\[Beta]1)^3); a2 = (\[Alpha]*k2)/(
DD2*(\[Beta]2)^3); a3 = (\[Alpha]*k3)/(
DD3*(\[Beta]3)^3); a4 = (\[Alpha]*k4)/(DD4*(\[Beta]4)^3);
; M0 h8 y) T. C9 b\[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, \
) c- l6 r9 ^# ?$ w& h1 x/ ?1]/Subscript[l, 4]; b3 = Subscript[l, 2]/Subscript[l, 3]; b4 = \
Subscript[l, 2]/Subscript[l, 4];
DD1 = D11/(A11*h^2); DD2 = D21/(A11*h^2); DD3 = D31/(A11*h^2);
3 s3 Y* I8 T/ P7 xDD4 = D41/(A11*h^2);
r11 = Power[a1*\[Omega]^2, (4)^-1]; r12 = Power[
a2*\[Omega]^2, (4)^-1]; r13 = Power[
, f- ^. A+ ^* p% N) }! P 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 O) @4 ?7 G' j7 r# e" D" _  {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,
7 M) [, C/ y) s' M   Sin[r12], 0, 0},
  {0, Sinh[r12], 0, 0, 0, Cosh[r12], 0, 0, 0, -Sin[r12], 0, 0, 0,
   Cos[r12], 0, 0},
* R8 L& Y* ]" Z4 Z  {-(r11)^2*DD1*(\[Beta]1)^2*Cosh[r11],
0 ~, W* v6 b5 p) e; B   0, (r13)^2*DD3*(\[Beta]3)^2, (r14)^2*DD4*(\[Beta]4)^2,
   -(r11)^2*DD1*(\[Beta]1)^2*Sinh[r11], 0, 0,
   0, (r11)^2*DD1*(\[Beta]1)^2*Cos[r11],
! p) N" R. |& z   0, -(r13)^2*DD3*(\[Beta]3)^2, -(r14)^2*DD4*(\[Beta]4)^2, (r11)^2*
    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],
' e; `* `2 P: M# i* e2 t   0, (r13)^3*DD3*(\[Beta]3)^3, (r14)^3*DD4*(\[Beta]4)^3, -(r11)^3*
    DD1*(\[Beta]1)^3*Sin[r11], 0, 0, 0,
: ~- `$ Q) w' j$ ~- n) t; {$ i- U   (r11)^3*DD1*(\[Beta]1)^3*Cos[r11],
7 D3 A. e, _1 Q& A) R" s3 x5 @   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*
    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*
9 E; [1 X/ h# q) `' Y) ~    Sinh[r14], 0, (r12)^2*DD2*(\[Beta]2)^2,
   -(r13)^2*DD3*(\[Beta]3)^2*Cos[r13], -(r14)^2*DD4*(\[Beta]4)^2*
) d2 C+ s; }8 a! h" r    Cos[r14], 0, 0,
   -(r13)^2*DD3*(\[Beta]3)^2*Sin[r13], -(r14)^2*DD4*(\[Beta]4)^2*
    Sin[r14]},
! Z/ \" n# ?& M9 H- Q5 P8 u. _  {0, 0, (r13)^3*DD3*(\[Beta]3)^3*Sinh[r13], (r14)^3*DD4*(\[Beta]4)^3*
* y3 F7 h) t3 f    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],
0 w0 }' P% M3 F1 W& t' @: \& A9 u9 v8 E   (r14)^3*DD4*(\[Beta]4)^3*Sin[r14],
   0, (r12)^3*DD2*(\[Beta]2)^3, -(r13)^3*DD3*(\[Beta]3)^3*Cos[r13],
; {( ~& P5 {6 O, L   -(r14)^3*DD4*(\[Beta]4)^3*Cos[r14]},
  {Cosh[r11], 0, -1, 0, Sinh[r11], 0, 0, 0, Cos[r11], 0, -1, 0,
   Sin[r11], 0, 0, 0},
* n6 S2 f* {7 C" y% D  {r11*Sinh[r11], 0, 0, 0, r11*Cosh[r11], 0, -r13*b1,
   0, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
   0, -r13*b1},
. J4 [9 O8 U! {6 g  {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,
2 E0 b4 |' v5 K   0, -r14*b2, -r11*Sin[r11], 0, 0, 0, r11*Cos[r11],
   0, 0, -r14*b2},
  {0, 1, -Cosh[r13], 0, 0, 0, -Sinh[r13], 0, 0, 1, -Cos[r13], 0, 0, 0,
. W: i0 @/ m- q    0, -Sin[r13]},
  {0, 0, -r13*b3*Sinh[r13], 0, 0, r12, -r13*b3*Cosh[r13], 0, 0, 0,
   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]},
* c* H% A) o8 x' _" N6 ?  {0, 0, 0, -r14*b4*Sinh[r14], 0, r12, 0, -r14*b4*Cosh[r14], 0, 0, 0,
7 h+ }* B' r$ U+ ^, q& Z   r14*b4*Sin[r14],
   0, r12, 0, -r14*b4*Cos[r14]}}
aa1 = Det[BB]
Plot[aa1, {\[Omega], -10, 10}]
1 {3 q, [( X0 _3 T, F% wss = FindRoot[aa1 == 0, {\[Omega], 2}][[1]]
* u: d% A0 i  w2 [5 }$ j. s
! z1 Q7 }- d! Q

百年孤独

这是什么语言编程？

游龙游龙

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

mathematical编的