clear all %clear memory
global alpha
k=1;
for alpha=1:0.1:4
s(k)=quadl('myfunc',0,1);
a(k)=alpha;
k=k+1;
end
plot(a,s)
