You are right.
On Wed, Oct 19, 2016 at 12:37 AM, Steven G. Johnson
wrote:
>
>
> On Tuesday, October 18, 2016 at 4:34:38 PM UTC-4, Michele Zaffalon wrote:
>>
>> quadgk(t -> cis(gamma(t)), 0, 1)
>>
>
> No, this is wrong because you forgot the Jacobian factor.
>
On Tuesday, October 18, 2016 at 7:05:17 PM UTC-4, Mosè Giordano wrote:
>
> In any case, I have to admit that quadgk is much more powerful than
> what I expected, at least because purely implemented in Julia, so can
> work with any Julia type.
>
quadgk also supports arbitrary-precision (BigFloat
On Tuesday, October 18, 2016 at 7:05:17 PM UTC-4, Mosè Giordano wrote:
>
> Hi Steven,
>
> 2016-10-19 0:36 GMT+02:00 Steven G. Johnson >:
> > For example:
> >
> > quadgk(z -> 1/z, 1, 1im, -1, -1im)
> >
> > integrates 1/z over a closed counter-clockwise diamond-shaped contour
> around
Hi Steven,
2016-10-19 0:36 GMT+02:00 Steven G. Johnson :
> For example:
>
> quadgk(z -> 1/z, 1, 1im, -1, -1im)
>
> integrates 1/z over a closed counter-clockwise diamond-shaped contour around
> the origin in the complex plane, returning 2πi by the residue theorem.
Did you mean
quadgk(z
On Tuesday, October 18, 2016 at 4:34:38 PM UTC-4, Michele Zaffalon wrote:
>
> quadgk(t -> cis(gamma(t)), 0, 1)
>
No, this is wrong because you forgot the Jacobian factor.
On Tuesday, October 18, 2016 at 4:27:22 PM UTC-4, digxx wrote:
>
> do u have an example for how to use a contour?
> quadgk(cis,0,1+1*im)=
>
probably integrates over the straight line so how can I integrate over the
> line gamma(t)=t+im*t^2
>
By contour, I just meant straight-line segments. Fo
quadgk(t -> cis(gamma(t)), 0, 1)
On Tue, Oct 18, 2016 at 10:27 PM, digxx wrote:
> do u have an example for how to use a contour?
> quadgk(cis,0,1+1*im)
> probably integrates over the straight line so how can I integrate over the
> line gamma(t)=t+im*t^2
>
do u have an example for how to use a contour?
quadgk(cis,0,1+1*im)
probably integrates over the straight line so how can I integrate over the
line gamma(t)=t+im*t^2
quadgk is fine for vector-valued functions (and in fact any integrand type
supporting +, -, *real, and norm, and can also integrate over infinite
intervals and contours in the complex plane.
But for multidimensional integrals you should use Cuba.jl or Cubature.jl or
GSL.jl or similar.
>
> Btw: Can quadgk also be used for complex functions? (the integration is
> still over a real range)
>
Just try it:
julia> quadgk(cis, 0, pi/2)
(1.0 + 0.im,2.482534153247273e-16)
However for more advanced integration (in particular for multidimensional
integration and vecto
Btw: Can quadgk also be used for complex functions? (the integration is
still over a real range)
thx
I think what you want is
g(s) = quadgk(t->f(s,t), 0, 1)
On Saturday, October 15, 2016 at 6:37:58 PM UTC-6, digxx wrote:
>
> having a function of the form f(s,t) defined is it possible to somehow
> tell quadgk to not evaluate until I supply a value for s while t should be
> the integration varia
Perhaps this is what you mean:
s = 1.0
quadgk(t -> f(s,t),0,1)
On Sunday, October 16, 2016 at 2:37:58 AM UTC+2, digxx wrote:
>
> having a function of the form f(s,t) defined is it possible to somehow
> tell quadgk to not evaluate until I supply a value for s while t should be
> the integratio
14 matches
Mail list logo