Re: Integration with quad

[Przemek Klosowski]

On 12/09/2009 02:07 PM, Thomas Göbel wrote:
> function y = rect(t)
> if t<=0.5
>      rect=1
> else
>      rect=0
> end
> endfunction
>
> For integration i use
>
> foo = quad('rect',-1,2)
I think you want this:

function y=rect(t); if t<=0.5 y=1; else y=0; end; endfunction

Your function does not assign to the result variable 'y' so it fails
on what I believe is the first call to it from 'quad()'.