
t=-10:0.001:25;
x1=10*(exp(-3*t)).*(heaviside(t-5)-heaviside(t-9));
x2=t.*(heaviside(t-6)-heaviside(t-14));
convolution=conv(x1,x2,'same');
 
figure;
subplot(3,1,1);
plot(t, x1);
ylabel('f(x1)');
subplot(3,1,2);
plot(t, x2);
ylabel('f(x2)');
subplot(3,1,3);
plot(t, convolution);
xlabel('Time');
ylabel('convolution');


