loop in integral plot

function y=f(theta)

L=10;

n1=1-l^2/2;

nc=1;

Neff(theta)=L/tan(theta);

y=((Neff(theta)*n1^2*sin(theta)*cos(theta)/(1-n1^2*cos(theta)^2)^2))/ (n1^2*sin(theta)*cos(theta)/(1-n1^2*cos(theta)^2)^2);

I1=integral(f,asin(nc/n1),pi/2);

disp(I1)

end

plot between l vs I1

[shiv gaur]

pl help to plot by varying l from any value to I1