quadl problem (Matrix dimensions must agree.).

hi all, that is my function that i want to integrate with quadl

function [f] = integranda3(p,px,t,x)

a=p(1)+(1i*p(5)*p(3)*px);                                       

b=sqrt(power(a,2)+((power(p(3),2))*((power(px,2))-(1i*px))));   

o11=(1i*px*x);

o12=(p(6)*((power(px,2)-(1i*px))/(a+(b.*(coth((b*t)/2))))));

o1=exp(o11-o12);

o2=exp((((-2*p(1)*p(2))/(power(p(3),2)))*log((cosh((b*t)/2))+((a/b)*sinh((b*t)/2))))+((p(1)*a*p(2)*t)/power(p(3),2)));

f=o1.*o2;

end

where p is a vector of 6 real numbers, px is de parameter of the integration (d/dpx), t is equal to 1 and x is a real number.

when i do

funzione=@(px)integranda3(p,px,t,x);

INT=quadl(funzione,-100,100);

i get an error:

Error using  - 
Matrix dimensions must agree.

Error in integranda3 (line 12)
o1=exp(o11-o12);

i tried with minus(o11,o12) or with other integration like quadcc but i cant manage to solve this errore.

while .* solve this problem for *, and power solve for ^, it seems that minus dont solve this problem for -

Thank you so much for any help.

V.M.

Answer 1

[bym]

You probably have to fully vectorize your calculations:

.*  % for multiply
./ % for divide
.^ % for power

[Valerio Matteucci]

it works thank you i forgot ./