# Tag Info

Contrary to what @user21's answer, I don't think that you need to do anything special for point loads. Let's see why. A point load can be represented as a Dirac delta "function". So, in your case it would be something like $$R = \rho \delta(x - x_i)\, ,$$ where $\rho$ is the intensity of the source and $x_i$ is the position. If we use a weighted ...
Seen that you also have an account over at mathematica stackexchange I am going to show an implementation using Mathematica. This is certainly not the only way you can do it but hopefully gets you started. We start by creating a 1D mesh Needs["NDSolveFEM"] region = Line[{{0}, {1}}]; includePoints = {{1/3}, {2/3}}; mesh = ToElementMesh[region, &...