intensity = function(data, ts, a, delta, beta) {
let index = 0;
const intensity = [];
for(let i = 0; i < ts.length; i++) {
const t = ts[i];
intensity[i] = {
t: t,
intensity: a + data.filter((x) => x.t < t).map((x) => x.g * Math.exp(-delta * (t - x.t))).reduce((a, b) => a + b, 0)
};
}
return intensity
}