How to add limit condition to existing code

How to add lim h ->0 condition to the below code. How is tril, toeplitz in code related to the above equation? How to edit the sum condition? Is there any simple form to write the above equation?

    h = 0.01    
    n  = numel(y);
    m  = 0:(n-1);
    s  = (-1) .^ m;

    M  = tril(ones(n));
    R  = toeplitz(y(:)');
    T  = meshgrid((gamma(a+1)/(h^a)) * s ./ (factorial(m).*gamma(a-m+1)));
    Y  = reshape(sum(R .* M .* T, 2 ), size(y));   

Answer 1

[Torsten]

If you are not satisfied with the approximation for h=0.01 in the code, you could evaluate the function fgl_deriv for two or three different values for h and evaluate the line or the parabola through the results at h=0.

Best wishes

Torsten.

Answer 2

[KALYAN ACHARJYA]

How to add lim h ->0 condition to the below code-

limit(main_function code)

detail read here

[Elysi Cochin]

i too went through it. But as the functions, tril, toeplitz, meshgrid, reshape

i'm not able to relate it

It would be of great help if you could help me edit the code