ClrText
"Newmark method"
"Make sure P=fn1(T)"
"K stiffness="?->K
"C damping="?->C
"M mass="?->M
"gamma="?->\r
"beta="?->\th
"Dt="?->D
"NumOfSteps="?->N
"V0="?->V
"U0="?->U
K+\r/(\th*D)*C+1/(\th*D^2)*M->K
1/(\th*D)*M+\r/\th*C->A
1/(2*\th)*M+D*(\r/(2*\th)-1)*C->B
ClrText
"K^=":K_
"A=":A_
"B=":B_
{N+1,10}->Dim Mat A
For 1->I To N
(I-1)*D->T
T->Mat A[I,1]
\fn1->P
P->Mat A[I,2]
Next
U->Mat A[1,10]
V->Mat A[1,9]
(Mat A[1,2]-C*Mat A[1,9]-K*Mat A[1,10])/M->Mat A[1,3]
'main loop
"1.Ti 2.Pi 3.U2i 4.DPi"
"5.DPi^ 6.DUi 7.DU1i"
"8.DU2i 9.U1i 10.Ui"_
For 1->I To N
(I-1)*D->T
Mat A[I+1,2]-Mat A[I,2]->Mat A[I,4]
Mat A[I,4]+A*Mat A[I,9]+B*Mat A[I,3]->Mat A[I,5]
Mat A[I,5]/K->Mat A[I,6]
\r/(\th*D)*Mat A[I,6]-\r/\th*Mat A[I,9]+D*(1-\r/(2*\th))*Mat A[I,3]->Mat A[I,7]
1/(\th*D^2)*Mat A[I,6]-1/(\th*D)*Mat A[I,9]-1/(2*\th)*Mat A[I,3]->Mat A[I,8]
Mat A[I,10]+Mat A[I,6]->Mat A[I+1,10]
Mat A[I,9]+Mat A[I,7]->Mat A[I+1,9]
Mat A[I,3]+Mat A[I,8]->Mat A[I+1,3]
Next
"finish"