Axel Vogt

5153 Reputation

19 Badges

16 years, 307 days
Munich, Germany

MaplePrimes Activity

These are replies submitted by Axel Vogt

Please post code instead of images.


'eval(%, omega=-omega) + %';

simplify(%) assuming 0<= omega:
evalc(Re(%))  assuming 0<= omega:
#map(identify, %):
vRe:=unapply(%, omega, t);

Now one wants to integrate vRe from 0 to 40 (instead of infinity), giving a function it t.

vRe(omega, 1); plot(%, omega=14.5 .. 14.6); # (t=1 as an example)

That plotting shows something which probably needs a Cauchy Principal Value and it is no clear (for me) why summation is good enough (omega = 10 is ok, it cancels out)

I do not quite understand why one can cut off at 0 and 40 and why that summing is a reasonable or rough approximation.

The integrand has singularities - the denominator only depends on omega and has zeroes, omega = 10 gives a pole for example.

You should post it as code and care for correct brackets

I do not quite understand your sheet, but having Sum(f(n), n= a .. b) and setting g:= k -> f(5*k) then Sum(g(n), n= a .. b) might be what you are looking for. So write your summands as function f.

Likewise and not caring for the discontinuity:

G:=w -> Int(eval(A1+A2,omega=w),theta=0..2*Pi, method = _d01ajc, epsilon=1e-3);

plot(G, 0 .. 50);


My point is: Maple uses the conventions (or I get you wrong), I do not want to discuss conventions


Here it may be simple/elegant - but what about running up to Pi + Pi/4 or Pi + Pi/2 ?


Maple uses the conventions in Abramowitz & Stegun or DLMF,

I am not aware of a commitment though, say through Maple's documentation

@Axel Vogt 

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/ .


@Carl Love 

Maple Worksheet - Error

Failed to load the worksheet /maplenet/convert/ .


@maple123 you already have 2 examples on treating such a task, I will not do the 3rd.

@Carl Love I see what you mean: integral = antiderivative + constant(y),
though the tasks ignored any constants and therefore are not quite wellposed

But ignoring such presentations would be a restriction for parametric integrals (and I was viewing at it as such)

The second question was also asked later by maple123 and answered by a different user. That post vanished (, but here is a sketch:

Write it as Int(diff(u(x,y),x)^2,x) + Int(diff(u(x,y),x$2)*u(x,y),x) = A + B.

Then IntegrationTools:-Parts(A, (diff(u(x,y),x))^1) + B gives you the desired result

1 2 3 4 5 6 7 Last Page 1 of 197