Alfred_F

Mr. Alfred Flaßhaar

590 Reputation

11 Badges

1 years, 331 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

@janhardo 

..."So I'm going to revise the calculation above." From your previous calculations, I could not determine the value of the fenced meadow. (Note: The value should be above 15200.) It should also be clear how the fence is closed at the eastern end along the x-axis – with or without a tangent. The proof that only symmetrical fence construction is possible is unnecessary here.

@dharr 

...the farmer's buying behavior seems odd – but he probably has his reasons.

When formulating the task, I had to compromise between precision and vivid clarity. That's probably why there was a misunderstanding, which @janhardo correctly cleared up - sorry.

@janhardo 

...but slightly more convex. Only on the oak tree is a corner allowed. The meadow should be as expansive as possible ;-).

@sand15 

...for Your detailed and helpful advice and the literature reference.

However, I'm not currently concerned with solving a specific practical problem. Rather, for other reasons, I want to construct a "pathological" example. It should demonstrate what theory has long known: Solving the variational problem as a minimum problem is equivalent to solving the resulting Euler differential equation as a necessary condition. My goal was to show, in an intuitive way, that the solution of Euler's differential equation corresponds, at least approximately, to the calculation of the functional minimum and to represent both in a single plot. I am particularly interested in calculating the minima of functionals, especially convex functions (see Rockafellar, Hiriart-Urruty/Lemarechal).

Since I've only recently begun to explore the world of Maple, there are several inaccuracies and errors in the command sequence of my file. Hence my plea for help.

Later I would like to try out so-called quadratic functionals (Lit.: Michlin "Variation Problems of Mathematical Physics") in Maple.

Addition:

Thank you again for your advice. The frequency analysis is particularly helpful. I still need to thoroughly understand the commands for it. Fifty years ago, this practical application (I know the theory) would have been very welcome.

@nm 

... is y(x)=const.=pi. Is this particular solution allowed?

@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

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