A small suggestion

[BasemRajjoub]

That's a nice project, but I noticed that you are using the "three moments method" and I thought that you can benefit a lot from using "singularity functions" and it's already been implemented in Sympy.
Please take a look here:

https://docs.sympy.org/latest/modules/physics/continuum_mechanics/beam_problems.html

Greetings

[buddyd16]

Thanks for the info, right now I am working thru various approaches to the analysis so when I eventually get to a stiffness analysis I have a few ways to validate results. Singularity functions would be another good method to review.