parabolic-arch code snippets in matlab

find tension, shear and moment on a pinned parabolic arch with vertical loading in matlab