我认为我已经找到答案了.在纸上做了一点数学后，我得到了:

I think that I have got the answer. After doing a bit of math on paper, I get this:

a=cumsum([a(1) ((b(1:end-1)+b(2:end)).*diff(c))./(d(1:end-1)+d(2:end))]);