Alfred_F

Mr. Alfred Flaßhaar

540 Reputation

11 Badges

1 years, 258 days
Brandenburg, Germany
As a retired individual with degrees from German universities in mathematics/analysis and structural engineering, I spent my professional life in responsible positions in research, teaching, and practical application, working on the mathematical modeling of states and processes in real-world systems. Now I have the time to explore interesting mathematical problems using Maple. It is my professional curiosity that drives me.

MaplePrimes Activity


These are replies submitted by Alfred_F

@nm 

...So, this is about the famous pendulum equation. Then choosing the initial value y(∞)=π is physically meaningless, because this deflection angle describes an indifferent equilibrium state, which, according to the derivation of the equation, does not capture this state. In contrast, the pendulum equation for deflections smaller than π is derived from equilibrium conditions or conservation laws. Overall, the textbook problem therefore appears to me to be in a confusing form. It would be better to ask: what limiting value allows us to estimate the solution – and that is π.

edited:

From my point of view, only the reversal of the problem statement according to the textbook makes sense: With an initial displacement of "close" pi, the mathematical pendulum will swing infinitely long.

Apart from the fact that calculating the solution of a second-degree ODE also requires two initial values, solutions exist, for theoretical reasons, only for finite abscissas of the initial values. This becomes "experimentally" obvious when the second-degree ODE is treated as a system of two first-degree equations with the same initial value.

p. s.

test.mw

@acer 

I wasn't familiar with "convert". Last question: How can the lhs be converted to the rhs in (4) of the attached file test2, instead of using the command "is" in the argument of "is"?

restart

kernelopts(version)

`Maple 2024.2, X86 64 WINDOWS, Oct 29 2024, Build ID 1872373`

(1)

simplify(convert(exp(u)/(1+exp(u))^2, trigh))

(cosh(u)+sinh(u))/(1+cosh(u)+sinh(u))^2

(2)

is(convert(exp(u)/(1+exp(u))^2 = (cosh(u)+sinh(u))/(1+cosh(u)+sinh(u))^2, exp))

true

(3)

is((cosh(u)+sinh(u))/(1+cosh(u)+sinh(u))^2 = 1/(4*cosh((1/2)*u)^2))

true

(4)

plot([exp(u)/(1+exp(u))^2, 1/(4*cosh((1/2)*u)^2)], u)

 

NULL

Download test2.mw

@dharr 

The goal of my question was not to solve the specific problem, it's easy to solve and an exact solution function is known. The background of my question is: How can an integral equation with an unknown function at the integration limit be solved in Maple if the integrand doesn't have a conventional/closed-form antiderivative? How can the coefficients of a linear combination of appropriately chosen functions be determined to obtain an approximate function?

Using the example chosen here, whose solution is known, I tried to test this using the `minimize` function. The result not only showed the coefficients but also a special value of t corresponding to the minimum. What options need to be added to `minimize` to obtain coefficients and then an approximate function over a t-interval?

The problem I remember with the unknown function at the integration limit dates back to the 1990s. The integrand had a terrifying structure; it was nonlinear and had no antiderivative. Back then, my only option was to use the Ritz method in Mathcad. Since then, hardware and software options have improved. Therefore, I'm attempting to solve this old problem again using these newer methods.

@Rouben Rostamian

...

is interesting:

It is identical to every special solution to the initial value of the form y(u)=u, where the additional conditions for u and y(u) must be met.

@Carl Love 

I would never have thought of that. I need to take a closer look at the 1D and 2D input methods and their differences. Where can I find more information about that?

I have got it. :-)

@Carl Love 

Your method is, of course, much easier to apply than my experiments with parameterized real and imaginary parts. It's interesting to observe, using examples, how the zeros of polynomials "arise" and paths in the complex plane of polynomials are mapped into this plane.

@Christian Wolinski 

As a Maple beginner, I'm not yet familiar with all the commands You use and the logical structure of nested code. Therefore, I'd appreciate your advice regarding the error message in the attached file. It doesn't seem to be the parentheses...

test.mw

@acer @nm

Learned something new again. Now, playing with polynomials in the complex plane is enjoyable. I am currently interested in how polynomials map given paths in the complex plane to this plane.

@nm 

...how can the polynomial be factored using radicals instead of decimals?

restart

with(PolynomialTools)

term := factor(x^4-12*x-12, complex)

(x+.9360298293)*(x+.8098950470+2.088088515*I)*(x+.8098950470-2.088088515*I)*(x-2.555819923)

(1)

convert(term, radical)

(x+.9360298293)*(x+.8098950470+2.088088515*I)*(x+.8098950470-2.088088515*I)*(x-2.555819923)

(2)

NULL

Download test.mw

@Kitonum 

...are easily understood with pen and paper. It's enough to apply the rules for powers and logarithms. The original equation only looks dangerous, but after simplification it is not. I managed to do that without Maple before. Maple is, of course, much faster. A plot for variants of C also looks quite nice.

Currently, the "surd" command is very interesting to me, as I want to explore the complex plane with Maple, so thank You for that.

@dharr 

Your very first sentence, with the remark about "default behavior," has already solved the problem. I will begin with the "residue counter" (residue theorem) and refresh my memory of the good-old book by Hurwitz and Courant.

@dharr 

...gives in Maple2024.2: [o,undefined] also - strange. 

Assumption:

ic:=y(infinity)=5*pi/4 is during the insertion test interpreted in Maple as the limit x-->infinity. This limit is not uniquely determinable, analogous to lim(sin(x), x-->infinity). Therefore, y(infinity) can be anything according to pi/4+k*pi. In any case, ic can be varied mod pi without changing the solution.

BTW:

I forgot to mention that while a general solution exists in the domain for x>0 according to the Picard-Lindelöf theorem, the proof presupposes the finiteness of the abscissa for the initial value. And y(infinity), however, lies, in a sense, on the boundary, which is not covered by the existence and uniqueness theorems. It is therefore not surprising that Maple leads classical solution methods to strange results.

@dharr 

"For your solution, odetest(sol, [ode, ic], y(x)) also gives [0,undefined]. "

With Maple 2024.2 I have got the attached results (with which I had played around ;-) )test.mw
 

restart

with(RealDomain)

assume(x > 0)

#Die Beschränkung x<>0 wird nach Auflösung der DGL nach y´ offensichtlich. Wegen des hier gewählten Anfangswertes muß x>0 gelten. The restriction x<>0 becomes obvious after solving the differential equation for y'. Because of the initial value chosen here, x>0 must hold.

ode := x^2*(diff(y(x), x))+cos(2*y(x)) = 1NULL

x^2*(diff(y(x), x))+cos(2*y(x)) = 1

(1)

#Folgende Lösung "mysol" wurde extern berechnet zum Anfangswert y(infinity)=5*Pi/4. Dieser exotische Anfangswert ist als Grenzwert für x-->infinity zu verstehen. Da der Cotangens periodisch ist, führt dieser Anfangswert nicht zu einer eindeutigen Lösung. Jede Ordinate Pi/4+k*2*Pi erfüllt die Forderung zu "mysol". Daher wählt Maple z. B. für (4) die Ordinate Pi/4 (blaue Kurve).         The following solution, "mysol," was calculated externally for the initial value y(infinity) = 5π/4. This unusual initial value is to be understood as the limit as x approaches infinity. Since the cotangent is periodic, this initial value does not lead to a unique solution. Every ordinate, π/4 + k*2π, satisfies the requirement for "mysol". Therefore, Maple chooses, for example, the ordinate π/4 (blue curve) for (4).

mysol := cot(y(x)) = 2/x+1

cot(y(x)) = 2/x+1

(2)

solve(cot(y(x)) = 2/x+1, y(x))

arccot((x+2)/x)

(3)

limit(arccot((x+2)/x), x = infinity)

(1/4)*Pi

(4)

#Der folgende Befehl liefert die zum Anfangswert y(infinity)=5*Pi/4 die korrekte Lösung (rote Kurve).             The following command provides the correct solution (red curve) for the initial value      y(infinity)=5*Pi/4.

sol := dsolve({ode, y(infinity) = (5/4)*Pi}, y(x))

y(x) = (3/2)*Pi-arctan((x+2)/x)

(5)

NULL

``

odetest(sol, ode, y(x))

0

(6)

plot([3*Pi*(1/2)-arctan((x+2)/x), arccot((x+2)/x)], x = 0 .. 100, color = ["Red", "Blue"])

 

NULL


 

Download test.mw

.

@dharr 

...the initial value is the key. Every initial value ordinate pi/4+k*2*pi yields the same solution. Without a "special command," the initial value does not lead to a unique solution. It appears that "sol := dsolve({ode, y(infinity) = 5/4*Pi}, y(x))" is required.

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