numerical integration

I have code

   c=0.2; t =  0.00001:0.0001:20; trapz(t,normcdf(((log(q(1)/K(1))+t*c^2/2)/(c*sqrt(t))),0,1)*stblpdf(t,0.5,1,12,0,'quick'));

but then I want to use vectors for 'q' and 'K' so I coded

c=0.2; t =  0.00001:0.0001:20; z = trapz(t,normcdf(((log(q./K)+t.*c.^2/2)./(c.*sqrt(t))),0,1).*stblpdf(t,0.5,1,12,0,'quick'));

but get the error

??? Error using ==> plus
Matrix dimensions must agree.

how can this be done?

Answer 1

[Walter Roberson]

You would get that error if the length of q is not the same as the length of t. Are you attempting to calculate for each pair q(J), t(K) ? If you are then see bsxfun() or ndgrid()