ecterrab

14540 Reputation

24 Badges

20 years, 22 days

MaplePrimes Activity


These are replies submitted by ecterrab

@NOh 

This seems to be a different question, although I do not see the images, probably a problem in Mapleprimes website. Could you please upload a worksheet instead so that I give this a look? Thanks.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Age 

I'm out of town in this moment, and unable to check things further, but I suggest you to update your Physics library from the R&D Maplesoft Physics page and try your examples again. I will return to this topic in one week.

Best

Edgardo

@w-Teilchen 

Unfortunately, I didn't have time to look at your post yet, I've been coding the tetrad formalism of General Relativity, it is mostly finished. Your post, together with 2 other posts made here in Mapleprimes that, like yours, I couldn't address faster, are first in the list of things I plan to give a look starting Oct/15.

Edgardo S. Cheb-Terrab 
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Alejandro Jakubi , @rlopez

The problem is in the use of alias, as expressed by nm, or otherwise in the use of 2D input mode versus 1D input mode, as indirectly expressed by you two.

If you use alias, Robert, Typesetting:-Supress is of no use and things happen as mentioned by nm in this post. If on the other hand you do not use alias, AND input in 2D mode, things happen as you mentioned, Robert, and diff(g, t) returns not 0; otherwise if you use 1D input mode as I understand you use, Alejandro, then diff(g, t) returns 0 unless you use alias, in which case we are back to square one: things happen as said by nm. Alternatively, PDEtools:-declare is free of all these nuances, working OK with or without alias, and regardless of using 1D or 2D input - I mentioned it for that reason; but then we do not have the dot to indicate derivatives.

This requires a fix in Typesetting, as suggested by nm.

Edgardo S. Cheb-Terrab 
Physics, Differential Equations and Mathematical Functions, Maplesoft

@nm 

Note that, unfortunately, the steps you would like to see automated are however not possible to be automated. When dsolve returns ODESolStruc it is because the computer algebra system as a whole does not know how to solve the reduced ODE. Therefore, one cannot automate 'returning a final solution' as you suggest in your third and fourth paragraphs.

Likewise, odetest can only test the solution you give to it. Any automatic replacement of your solution or of your ODE by something else would be inappropriate. And the option useInt is an option of dsolve, not odetest, so this too cannot be automated. Think more about: you pass a solution and instead of testing it, odetest, automatically, would first compute another solution using dsolve's option useInt, test this other solution, then return 0. Would you accept that result? I wouldn't. It didn't test the solution I passed to it.

All this points to a relevant thing: computer algebra is of enormous value but you should not think of it as a "complete solution". There is no complete solution. There will always be cases that require your human-skilled approach. Whether you have the skillls to do better than the computer is another story, with which I agree with you, few people have those additional skills, but that doesn't change the facts, there will always be situations where you will need to guide the computer.

Yes, more examples in the help pages are always useful. I will see if I can add some of these of this post to the help page of odetest.

Edgardo S. Cheb-Terrab 
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Markiyan Hirnyk 

Hi Markiyan, y=0 is not a general solution, because it has not the necessary number of arbitrary constants. Therefore is not the solution expected from dsolve according to its help page. The fact that dsolve, in addition, returns essential singular solutions doesn't change that fact. By the way y = 0 is also not an essential singular solution of the ode posted, nor a singular solution of any kind, but just a particular solution that can be obtained taking one of the integration constants equal to 0 in the general solution.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Mac Dude 
I am sorry but the update for Maple 18 does not work in previous releases. 

Edgardo S. Cheb-Terrab 
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Markiyan Hirnyk 

Look closer into ?dsolve,details. In the Section entitled "Arbitrary Constants, General, Particular and Singular ODE solutions", you read: "By default, dsolve computes the general and essentially singular solutions of a nonlinear ODE; the latter are the singular solutions that cannot be obtained from the general solution by specializing the integration constants in any way."

In the same section, you read: "A general solution to an ODE depends on as many arbitrary constants as the differential order, say N. Conversely, any solution depending on N arbitrary constants is a general solution. Depending on N constants here means there is no possible redefinition of the _Cn in terms of other N independent constants that results in a form of the solution with fewer than N constants."

I suppose you missed this section when reading the page, but there is no possible misunderstanding: a solution of the form y(x) = 0, without arbitrary constants for a 2nd order ODE does not match these statements quoted from the help page ?dsolve,details.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Markiyan Hirnyk 

You say "The only closed form solution of the nonlinear ODE under consideration it can find is y(x)=0", but if you read the help page (?dsolve or ?dsolve,details) you see that dsolve is expected to always return a general solution , or otherwise return NULL. So this solution y(x) = 0 is not according to what is written in the help page, and so it is wrong even if mathematically a correct particular solution.

The actual problem was in one of the internal routines for computing solutions to nonlinear 2nd order ODEs in terms of elliptic JacobiSN functions, a condition being tested required an additional simplification, it is fixed now, the fix is available as usual in the Maplesoft R&D webpage for Differential Equations and Mathematical Functions.

Besides that, you could see that dsolve can compute a general solution for this equation if you try the Lie symmetries option, as in

 

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

Hi

I'm rather busy in this moment for several weeks ahead;let's see if I can give this one a look carefully and return useful feedback sooner than that.

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Michael_Watson 

Maple_Question_7.8.14_(reviewed).mw

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Alejandro Jakubi 

Just to clarify: the R&D udpates for Physics, Differential equations and Mathematical functions, are official Maplesoft updates. These  R&D updates pass through all the related Maplesoft internal test suites before being posted in the R&D pages for download. These R&D updates are however different from the more general dot updates in that R&D ones have not passed through beta testing as the dot updates do, and for that reason the R&D updates mechanism is so much more agile. In Physics, for instance, we frequently have one or more updates per day.

Edgardo S. Cheb-Terrab 
Physics, Differential Equations and Mathematical Functions, Maplesoft

@oldstudent 

I am not sure about your question: if you refer to the "DEs and Mathematical Functions Updates" mentioned in the title of this post, this update is already available for download for everybody in the Maplesoft R&D Differential Equations and Mathematical Functions webpage. This update, together with the update of the Maple Physics package are currently distributed only as a download from these two webpages, so you do not get them via tools -> check for updates.

Edgardo S. Cheb-Terrab 
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Alejandro Jakubi 

The status is: work in progress. I preferred to start where the gaps were bigger, clearly the 42 elliptic functions. Not what you asked but anyway some details: these "elliptic functions" include the 13 JacobiPQ, 13 InverseJacobiPQ, JacobiZeta, the 4 JacobiTheta, the 4 Weierstrass, and then the more familiar F, E, E', K, K', Pi and Pi' - 7 functions for which, in addition, we have an issue with their definition implemented in Maple in the 90's: it is non-standard, making things more difficult, because there is basically no literature to consult that uses these definitions.

There is always the issue of 'backwards compatibility" to consider but, generally speaking, the natural thing would be to go ahead redefining these F, E, E', K, K', Pi and Pi' elliptic functions according to the literature, resolve any differences by always following the modern and thorough NIST Digital Library of Mathematical Functions project, and finish with this historical issue in one go. The NIST project is the XXI century continuation of the work by Abramowitz & Stegun, Gradshteyn & Ryzhik, Bateman, Byrd & Friedman, P.B.M, etc..

Edgardo S. Cheb-Terrab
Physics, Differential Equations and Mathematical Functions, Maplesoft

@Alejandro Jakubi 

Thanks; it's fixed now. Regards.

Edgardo

First 47 48 49 50 51 52 53 Last Page 49 of 64