The method of Macaulay makes use of singularity functions to describe the integrations in the Euler-Bernoulli beam equation in a single equation, instead of splitting the structure into parts with equal load situations. In this report, an expansion for this method is derived, which makes it possible to solve structures with multiple stiffnesses and varying stiffness slopes using this method as well.