Re: [Maxima-discuss] dblquad fail

[Carlo de Falco]

Il 28 lug 2016 5:50 PM, "Ether Jones" <address@hidden> ha scritto:
>
>
> Thank you.
>
> This works:
> octave-3.6.4.exe:3>  I = dblquad(@(r,q)abs(
> r.*sqrt(r.^2+1-2*cos(2*pi*q))),0,1,0,1)
> I =  0.586626703381746
>
> ... but this doesn't:
> octave-3.6.4.exe:4>  I =
> dblquad(@(r,q)r.*sqrt(r.^2+1-2*cos(2*pi*q)),0,1,0,1)
> error: quadcc: integrand F must return a single, real-valued vector
> error: called from:
> error:   K:\OctaveVS\share\octave\3.6.4\m\general\dblquad.m at line 72,
> column 10
> error:    at line -1, column -1
> error: quadcc: integrand F must return a single, real-valued vector
> error:   K:\OctaveVS\share\octave\3.6.4\m\general\dblquad.m at line 65,
> column 5
>
> Why is the abs() necessary?

The error message you reported says:
error: quadcc: integrand F must return a single, real-valued vector

Is your integrand function real for any value of (r,q) in (0,1)?

c.