There is problem in Index values in calculatepinleft

The pin left index values should be for ja = 2 and jb = 4, but in eq-1 (in attached screenshot) it is shown index 3 for ja and index 5 for jb . so please provide the clarity on this.

Thanks & Regards,
Ajay

Hey @jyskrwr, the indices are correct as shown.

Remember, these index values correspond to locations in the overall primary stiffness matrix. This overall matrix for the structure still allows for a moment at the pin, until we remove that DoF at the top of the next code block…

removedDoF = restrainedDoF + pinDoF

Seán