I suggest the best way to verify this is to write the equation that you expect and check that ksdensity does what you want.
According to my interpretation of what you want, ksdensity does what you ask. But your formula has the bandwidth b in it. You would have to provide that to ksdensity if you want to control it.
Also, your paper seems to talk about truncating the density to some interval (say [0,1]) and then dividing by the integral over that interval. I would expect that to yield a biased estimate. Points in the interior of the interval get their full kernel included. Points near the edge get only a part of their kernel included.
The ksdensity function uses a different method to account for finite support of the distribution. You may want to look more at that. I know there are different ways to approach it, but the function implements one of those approaches.
