| Author |
Comment/Response |
Forum Moderator
email me
 |
 09/28/98 3:00pm
>
> Hi,
> I have tried in vain to compile the following
> \!\(\[Nu][k_] := \
> \((v = Round[{9, 3, 1} . k\/K0] + 14;
> If[v \[GreaterEqual] 1\ && \ v \[LessEqual] 27, v, 14])\);
> S\_l = Table[k = \[Xi][i];
> If[EvenQ[Round[
> \(k\[LeftDoubleBracket]1\[RightDoubleBracket] +
> k\[LeftDoubleBracket]2\[RightDoubleBracket] +
> k\[LeftDoubleBracket]3\[RightDoubleBracket]\)\/K0]],
> \(-c\_l\)\ InvSq[k], 0. ], {i, 27}]; \ S\_H = \(- .9882\) S\_l\ ;
> \[Mu] = \ 44.7\ e;
\[Kappa][kn_] :=
> DiagonalMatrix[\(c\_2\) Table[ke = kn - \[Xi][i]; ke . ke, {i, 27}]];
> offD := Table[
> If[i \[Equal] j, 0. , v = \[Nu][\[Xi][i] - \[Xi][j]];
> S\_H\[LeftDoubleBracket]v\[RightDoubleBracket] +
> S\_l\[LeftDoubleBracket]v\[RightDoubleBracket]], {i, 27}, {j, 27}];
> NOffD = offD;
> G[kn_] :=
> Module[{Eval, Evec},
> \t\t{Eval, Evec} = Eigensystem[NOffD + \[Kappa][kn]];
\t\t
> fdDiag = DiagonalMatrix[
> Table[fd[Eval\[LeftDoubleBracket]i\[RightDoubleBracket], \[Mu]], {
> i, 27}]];
\t\
> 0.5 \((Transpose[Evec] . fdDiag . TakeColumns[Evec, {14, 14}])\)
\t];
>
calcSEn :=
> Module[{MK}, MK = Table[0. , {i, 27}];
\t\t
> For[h = 1, h \[LessEqual] qN, \(h++\),
> For[i = 1, i \[LessEqual] qN, \(i++\),
> For[j = 1, j \[LessEqual] qN, \(j++\), \
\t\t\t\t\t
> MK += \(xw\[LeftDoubleBracket]h
> \[RightDoubleBracket]\)\[LeftDoubleBracket]2
> \[RightDoubleBracket]
> \(xw\[LeftDoubleBracket]i
> \[RightDoubleBracket]\)\[LeftDoubleBracket]2
> \[RightDoubleBracket]
> \(xw\[LeftDoubleBracket]j
> \[RightDoubleBracket]\)\[LeftDoubleBracket]2
> \[RightDoubleBracket]
> G[{\(xw\[LeftDoubleBracket]h
> \[RightDoubleBracket]\)\[LeftDoubleBracket]1
> \[RightDoubleBracket],
> \(xw\[LeftDoubleBracket]i
> \[RightDoubleBracket]\)\[LeftDoubleBracket]1
> \[RightDoubleBracket],
> \(xw\[LeftDoubleBracket]j
> \[RightDoubleBracket]\)\[LeftDoubleBracket]1
> \[RightDoubleBracket]}]]]];
\t\t\t\t\
> S\_H = Table[
> If[i == 14, c\_n\ MK\[LeftDoubleBracket]i\[RightDoubleBracket], \
> MK\[LeftDoubleBracket]i\[RightDoubleBracket] InvSq[\[Xi][i]]\ 2.
> \ c\_0], {i, 27}];
\t\t\ \ \ S\_H = Flatten[S\_H];
> S\_H = 1. \ Chop[S\_H, 10\^\(-30\)];
\t\ \ \ \ NOffD = offD;
\t]; \)
>
> Also for the last function, I tried using Map to do the integration of
> a vector function and had problems.
> Any suggestions?
> regards
> Mukund
>
=====
I have been unable to parse your code into sensible Mathematica input. Generally, inline notation such as the code you pasted above is difficult to
communicate with. Please repost your code in InputForm or as a notebook. If it does not fit in the message box, you can email it to student-support@wolfram.com
Forum Moderator
URL: , |
|
