a matrix with a dirac delta function

I am having trouble multiply a matrix with a dirac comb.

My matrix is the following

and I am to multiply? it by this Dirac function

to get this:

.

If I just multiply the two together, then I will just get 0's, but I shouldn't. This is my code currently.

for j = 1:n + 1
    u_approx_dirac = u_approx * dirac(x - c - j*delta);
end

x is a vector of 3001 elements,

t is a vector of 6001 elements.

Also, c = a for my code.

Please help me directly and not send me to other links. Thank you.

Answer 1

[Walter Roberson]

Guessing:

for j = 1:n + 1
    u_approx_dirac(:,j) = u_approx(:,j) .* dirac(x - c - j*delta);
end