How to obtain the numerical integration of a set of different values of prameter which appear in Xmin

(k,x,y) is a generated integrand as the following codes. v is a parameter in TuuV. I want to get a set of value of the integration "M=arrayfun(@(D) integral3(@(k,x,y) TuuV(k,x,y,D),kc.*sqrt(mu+delta-v),ku,0,pi/2,0,pi/4,'reltol',1e-4),v);" with different values of v. And thus I set v=0:0.01:0.5 in the codes as following. However the codes show the error "Xmin must be floating point scalar". (If I set v=0.5 in the codes, the problem disappears) How to solve this problem? Many thanks!

The codes are as following:

clc
clear
close all
mu=11;
delta=10;
vh=16;
HBAR=1.05457266e-34;
ME=9.1093897e-31;
ELEC=1.60217733e-19;
kc=sqrt(2.*ME.*ELEC./HBAR.^2);
ku=kc.*sqrt(mu+delta);
v=0:0.01:0.5;
M=arrayfun(@(D) integral3(@(k,x,y) TuuV(k,x,y,D),kc.*sqrt(mu+delta-v),ku,0,pi/2,0,pi/4,'reltol',1e-4),v);
fid=fopen('D:\ceshi.txt','wt');
fprintf(fid,'%g\n',M);
fclose(fid);

Answer 1

[Torsten]

M=arrayfun(@(D) integral3(@(k,x,y) TuuV(k,x,y,D),kc.*sqrt(mu+delta-D),ku,0,pi/2,0,pi/4,'reltol',1e-4),v);