Suppose to have these data

The simulink model is calulating
you can calulate the integral in the discrete form with cumtrapz
% data
t = linspace(0,1,20).';
x = 2*t;
y = t.^2;
c = 3;
% function to integrate
dz = c*(x-y);
% function z
z = cumtrapz(t,dz);
% plot
figure,hold on
plot(t,z,'LineWidth',2)
% check
plot(t,3*t.^2-t.^3,'--','LineWidth',2)


