quad, Warning: Infinite or Not-a-Number function value encountered.

Good evening, Sorry, my previous message contains some french word

I would like to plot (t,f(t)), with

 f(t)=integral_0^ (2*pi*1.0e+12)*75.2933 y(t,w) dw

For that, i used 'quad' to compute my integral, but there is a problem when matlab run: Warning: Infinite or Not-a-Number function value encountered.

Thanks i advantage for you help.

function y=toltallA(w,T) It not necessary to see between begin parameter and end parameter.

%%Begin parameter L=2.9e-6; v= 21.6e+3; a=0.142e-9; N_0=6.022e+23; A= 3*sqrt(3)*a^2*N_0/4;% w_max=(2*pi*1.0e+12)*75.2933; delta_M=3; B_N=3.85e-25; %% B_U=7.7e-25;%% alpha=3; %% theta=1000; c_d=1.0e-6; h=6.626e-34; hbar= h/(2*pi); kB=1.38*1.0e-23; coef=hbar^2/(4*pi*kB*v^2); %end of my parameter %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

tau=1./(v/L+ (A*delta_M^2*2*pi*c_d)/(2*pi*v^2*w_max^2)*w.^4+(B_N+B _U*exp(-theta/(alpha*T)))*T^3*w.^2);%%

%% this my integrand %%%%%%%y(t,w)%%%%%%%%%%%%%

y= coef*v^2/T^2*tau.*w.^3.*exp(hbar*w./(kB*T))./(exp(hbar*w./(kB*T))-1).^2;

%%% MAin Program : I use quad to approximate my integral.

T=1:10:1000;

QlA=zeros(size(T));

for k=1:length(T)

QlA(k)= quad(@(Z)toltallA(Z,T(k)),1.e-20,(2*pi*1.0e+12)*75);

end

thanks

[Walter Roberson]

duplicates http://www.mathworks.co.uk/matlabcentral/answers/69727-infinite-or-not-a-number-function-value-encountered