Informazioni sull'utilizzo di quad in MATLAB

Question

Si verifica uno strano problema con la funzione quad. Stavo usando quad per calcolare semplice integrale, ed ha funzionato per 10 a 20 volte, poi Matlab emette il seguente errore:

Error using quad (line 75) 
The integrand function must return an output vector of the same length as the input vector. 
yteor(k) = quad(@(q)(exp(-(q).^2).*q.^2/(k.^2+1)), 0, 1); 

Qui q e k sono scalari. Non riesco a capire cosa è sbagliato e perché ha funzionato diverse ore fa.

Modifica

Ecco il mio codice

for k=1:100, 
    xteor(k)=step*k; 
    yteor(k)=quad(@(q)(exp(-(q).^2).*q.^2/((step.*k+1).^2)),0,1); 
end plot(xteor,yteor,'r'); 

Answer 1

Il seguente frammento funziona per me su Octave (Matlab GNU clone)

step = 1; 
xteor = zeros(100,1); 
yteor = zeros(100,1); 

for k=1:100, 
    xteor(k)=step*k; 
    yteor(k)=quad(@(q)(exp(-(q).^2).*q.^2/((step.*k+1).^2)),0,1); 
end 
plot(xteor,yteor,'r'); 
pause 

La mia ipotesi è che il vostro errore è il conseguenza di qualcos'altro accadeva in precedenza nel tuo codice (forse relativo a step non essendo uno scalare?). Invece di concentrarsi su questa linea dove sorge l'errore. Prova a cercare cosa hai modificato appena prima che l'errore venga visualizzato.