MATLAB: How to plot a function obtained after integrating a two-parameter function with respect to one parameter

plot

y=0:0.01:1;
f=@(t,y) exp(-t.*y);
q=@(y) integral(@(t) f(t,y),0,2);
g=@(t,y) sin(t.*y);
p=@(y) integral(@(t) g(t,y),2,4);
k=q(y)+p(y);
M=cos(y);
L=M.*k
plot(y,L)
The actual functions are a bit different but they have the same parameters and I have followed the same procedure. I am neither getting values of L nor the plot. I also tried ezplot and fplot. ezplot shows error and fplot gives a blank plot.

Best Answer

Use symbolic toolbox
syms y t
f=exp(-t*y)
q=int(f,t,0,2)
p=int(g,t,2,4)
k=q+p
M=cos(y)
L=M*k
fplot(y,L,[0 1])