"STEPS:":?->S
Dim Mat M
List Ans[1]->N
{N,2}->Dim Mat W
{N,N}->Dim Mat F
{N,1}->Dim Mat R
"INVERSE K:"_
Mat K^-1_
Mat Ans*Mat M->Mat D
"D=K^-1*M:"_
For 1->I To N
Mat D_
Fill(1,Mat R)
For 1->J To S
Mat D*Mat R:Mat Ans[1,1]->L
Mat D*Mat R/L->Mat R
"Li:":L_
"Ri:"_
Mat R_
Next
\sqrt(1/L)->Mat W[I,1]
2*\pi/Mat W[I,1]->Mat W[I,2]
"Omega:":Mat W[I,1]_
For 1->K To N
Mat R[K,1]->Mat F[K,I]
Next
(Mat R*Trn Mat R)*Mat M->Mat X
(Trn Mat R)*Mat M*Mat R->Mat Y
(Identity N)-Mat X/Mat Y[1,1]->Mat S
"Ri.Tr(Ri).M:"_
Mat X_
"Tr(Ri).M.R:":Y_
"S=I-RRM/RMR:"_
Mat S_
Mat D*Mat S->Mat D
"D=D*S:"_
Mat D_
Next
ClrMat R
ClrMat X
ClrMat Y
ClrMat D
ClrMat S
"finish"