Dim Mat M
List Ans[1]->N
{N,2}->Dim Mat H
{N,1}->Dim Mat R
{N,6}->Dim Mat L
{N,N}->Dim Mat S
{N,N}->Dim Mat T
{N+1,N}->Dim Mat U
{N+1,N}->Dim Mat V
{N+1,N}->Dim Mat A
{N,1}->Dim Mat I
"height:"
For 1->I To N
?->Mat H[I,1]
Next
"PseudoAn:"
For 1->I To N
?->Mat H[I,2]
Next
"i factors:"
For 1->I To N
?->Mat I[I,1]
Next
"Lnh=Sum(fi.M.i)"
For 1->I To N
For 1->J To N
Mat F[J,I]->Mat R[J,1]
Next
Trn Mat R*Mat M
Mat Ans*Mat I
Mat Ans[1,1]->Mat L[I,1]
Mat L[I,1]_
Next
"Mn=Sum(fi.M.fi):"
For 1->I To N
For 1->J To N
Mat F[J,I]->Mat R[J,1]
Next
Trn Mat R*Mat M
Mat Ans*Mat R
Mat Ans[1,1]->Mat L[I,2]
Mat L[I,2]_
Next
"Gn=Lnh/Mn:"
For 1->I To N
Mat L[I,1]/Mat L[I,2]->Mat L[I,3]
Mat L[I,3]_
Next
"si=M.fn.Gn:"_
For 1->I To N
For 1->J To N
Mat F[J,I]->Mat R[J,1]
Next
(Mat M*Mat R)*Mat L[I,3]
For 1->J To N
Mat Ans[J,1]->Mat S[J,I]
Mat S[J,I]*Mat H[I,2]->Mat T[J,I]
Next
Next
Mat S_
ClrMat S
"Fi=si.An:"_
Mat T_
ClrMat T
"Mnstar=Gn*Lnh:"
For 1->I To N
Mat L[I,3]*Mat L[I,1]->Mat L[I,4]
Mat L[I,4]_
Next
"Lntheta=hi.M.fn:"
(Trn Mat H*Mat M)*Mat F
For 1->I To N
Mat Ans[1,I]->Mat L[I,5]_
Next
"hnstar=Lntheta/Lnh:"
For 1->I To N
Mat L[I,5]/Mat L[I,1]->Mat L[I,6]_
Next
"uj=Gn*fn*An/omega^2:"_
For 1->I To N
For 1->J To N
(Mat L[I,3]*Mat H[I,2]/Mat W[I,1]^2)*Mat F[J,I]->Mat U[I,J]
(Mat L[I,3]*Mat H[I,2]*Mat F[J,I]->Mat A[I,J]
Next
Next
For 1->I To N
0->S
0->T
For 1->J To N
Mat U[J,I]^2+S->S
Mat A[J,I]^2+T->T
Next
\sqrtS->Mat U[N+1,I]
\sqrtT->Mat A[N+1,I]
Next
Mat U_
"aj=Gn*fn*An:"_
Mat A_
ClrMat A
"Vb=Mnstar*An:"
0->S
For 1->I To N
Mat L[I,4]*Mat H[I,2]_
Ans^2+S->S
Next
\sqrtS_
"Mbn=hnstar*Vb:"
0->S
For 1->I To N
Mat L[I,6]*Mat L[I,4]*Mat H[I,2]_
Ans^2+S->S
Next
\sqrtS_
"Duj=uj-uj-1:"_
For 1->I To N
0->W
For 1->J To N
Mat U[J,I]->S
If I<>N:Then Mat U[J,I+1]->T:Else 0->T:IfEnd
S-T->Mat V[J,I]
Mat V[J,I]^2+W->W
Next
\sqrtW->Mat V[N+1,I]
Next
Mat V_
ClrMat H
ClrMat R
ClrMat L
ClrMat U
ClrMat A
"finish"