## 20 Reputation

5 years, 206 days

## Thank you...

@Kitonum Thanks a lot for correcting the syntax.

## RE: integration over 3d vectors...

@Carl Love OK, I have updated the worksheet for a better example.

1-D and 3-D integrals of this type typically have different behaviours.

## RE: Integration involving dot products o...

@Arny, the following could be a start, but not much of an output from Maple.

 >
 >
 (1)

 >
 (2)

 >
 (3)

 >
 (4)

 >
 (5)

## @Carl Love Thanks, I corrected the ...

@Carl Love Thanks, I corrected the input.
Apparently, no closer to the result though.

## p, q, w being 3-d vectors...

Thank you tomleslie 3149 and Rouben Rostamian for taking interest.

The integration is is indeed over (dp,dq,dw) ∈ R3×R3×Ri.e. as indicated dp = d3p, where p is really \vec{p}, p ∈ R3 etc.

Is it possible to help Maple understand if dp were to be directly expressed in as a 3-vector differential?

It is possible to express the integrals in spherical co-ordinates, and without phi dependence, would be left with integration over six dimensions, with no vector variables left in the integral.

## Re: Integration (analytic or numeric) in...

@ecterrab
Thank you once again. The updated worksheet, with the actual integrand, variables of integration and limits is attached.

The integration is being attempted first analytically, to the extent possible. It is very likely that only numerical evaluations will be possible, at least for some of the integration variables. The ideal range of integration is -Infty to +Infty for each of the variables p1, p3 and p4, but it could be tested for convergence over limits -1 to +1, -10 to +10 etc. There is no integration over vector "r", which could be taken as 1 or its equivalent in Cartesian co-ordinates. In Cartesian co-ordinates, this is a nine dimensional integral. However, the dimensions might come down if polar co-ordinates are used, if they are assumed to be 3-vectors. This is not attempted here.

For comparison, a Mathematica output is also attached as an image. The Mathematica output seems to evaluate the scalar triple product directly, which Maple doesn't do right away (most likely because it's not built into the output of "CrossProduct[v1,v2] etc).

## Maple worksheet attached...

Attached is the Maple worksheet

VectorIntegral.mw

The integration (whether numeric, or analytic, hopefully both) needs to be done avoiding the singularities. There is no integration over r (i.e. \vec{r}).

Numerical integration, after any analytic integration, would need to be done for r = 1, 1.5, ....10 (so r={1,1,1} etc.) in cartesian coordinates over p1, p3, etc from -Infinity to + Infinity.
Later, there would be other integrations, over p1 etc. The full integral would have terms in the denominator like (p1 - p4)^2, (1+p4)^4, etc. and inetgrate over p4 etc. Due to reasons of angular dependencies, it's best if cartesian coordinates are used.

## extending the integration (using Maple 2...

Using the methodology you kindly provided, and using Maple 2017, I tried to extend the method to an integration that needs to be done (to the extent possible, in closed form, and then numerically) over the following expression over vectors p1, p3, p4, and p5. (My guess is this must be a combination of a rookie error, plus not knowing how exactly to hand the integral in Maple).

 >
 (1)
 >
 (2)
 >
 (3)
 >
 (4)
 >
 (5)
 >
 (6)
 >
 (7)
 >
 (8)
 >
 (9)
 >
 (10)
 >
 (11)
 >
 >

## Thanks - and a rookie question...

Thank you!

My Maple 17 worksheet is giving the following errors, is there something incorrect? (Tried correcting the p3 entry in the worksheet you had attached)

 >
 >
 (1)
 >
 (2)
 >
 (3)

This is the integratal to be evaluated

 > 'exp(i*(p__1_ . r_))/((p__3_ - p__1_)^2)';
 (4)

Note the ' ' delaying the evaluation, if you allow evaluation, the integrand is

 >
 (5)

Define iota

 >
 (6)

Because (4), has the value of i before being the imaginary unit, I need to input the integrand again

 > 'exp(i*(p__1_ . r_))/((p__3_ - p__1_)^2)';
 (7)

Now on the integration

 >
 (8)
 >
 (9)

In the above, there is the product of three Dirac delta functions, that can be represented as a single 3D Dirac delta

 >
 (10)
 >

 >
 >
 (1)
 >
 (2)
 >
 (3)

This is the integratal to be evaluated

 > 'exp(i*(p__1_ . r_))/((p__3_ - p__1_)^2)';
 (4)

Note the ' ' delaying the evaluation, if you allow evaluation, the integrand is

 >
 (5)

Define iota

 >
 (6)

Because (4), has the value of i before being the imaginary unit, I need to input the integrand again

 > 'exp(i*(p__1_ . r_))/((p__3_ - p__1_)^2)';
 (7)

Now on the integration

 >
 (8)
 >
 (9)

In the above, there is the product of three Dirac delta functions, that can be represented as a single 3D Dirac delta

 >
 (10)
 >

## Re: The Physics:-Vectors package...

with(Physics[Vectors]);
print(output redirected...); # input placeholder
[&x, +, ., ChangeBasis, Component, Curl, DirectionalDiff,

Divergence, Gradient, Identify, Laplacian, Nabla, Norm, Setup,

diff]

X := _i*x1+_j*x2+_k*x3;
print(output redirected...); # input placeholder
_i x1 + _j x2 + _k x3

P := _i*p1+_j*p2+_k*p3;
print(output redirected...); # input placeholder
_i p1 + _j p2 + _k p3

e^(X.P);
print(output redirected...); # input placeholder
(p1 x1 + p2 x2 + p3 x3)
e

int(e^(p1*x1+p2*x2+p3*x3), P);
%;
Error, (in int) integration range or variable must be specified in the second argument, got _i*p1+_j*p2+_k*p3


Not the solution expected?

 Page 1 of 1
﻿