1 Pages (6 items)
[Custom Functions] Substitution in expressions with dummy variables (diff, sum etc.) - Messages
#1 Posted: 7/26/2023 4:03:12 PM
Consider the following example. Substitution using at() is applied to a symbolic expression. Once f(t) is not defined, the derivative can't be calculated. It is nice, that t is not just replaced by 0 but the result anyways is not correct. The return value should be the unmodified input expression if evaluation/simplification fails.
You see the problem if you evaluate input and return value for a particular f(x), the results aren't equal.

Blatt20.sm (5 KiB) downloaded 62 time(s).
You see the problem if you evaluate input and return value for a particular f(x), the results aren't equal.
Blatt20.sm (5 KiB) downloaded 62 time(s).
Martin Kraska
Pre-configured portable distribution of SMath Studio: https://en.smath.info/wiki/SMath%20with%20Plugins.ashx
#2 Posted: 7/26/2023 5:52:02 PM
I don't know what the intention of the author of the at() function is, but what I can observe in other CAS is that there are always something like two versions for the substitution, in which perhaps one maybe returns the expected results.


My opinion is that the only procedure about which it can be said that the substitution is correct or not is the limit operation

Maybe Maxima can redefine at() using limits.
Best regards.
Alvaro.
My opinion is that the only procedure about which it can be said that the substitution is correct or not is the limit operation
Maybe Maxima can redefine at() using limits.
Best regards.
Alvaro.
1 users liked this post
Martin Kraska 7/26/2023 6:13:00 PM
#3 Posted: 7/26/2023 6:10:09 PM
Whatever I do on Maxima side, if the result contains a symbolic non-simplificable substitution (can happen with laplace transforms) then I'd like to display that using the at() function with it's nice operator form to display the result. Yet, as of now, the function isn't robust enough for this.
In the example below you see that the return value from Maxima contains an expression with at(), which is removed on SMath postprocessing side even if optimization=none.
It is part of the test suite for Maxima in the interactive handbook or on github.
![2023-07-26 23_05_43-SMath Solver - [Plugin Maxima Test Suite.sm_].png](/en-US/file/AP5nCG/2023-07-26-23_05_43-SMath-Solver---_Plugin-Maxima-Test-Suite_sm___png)
The limit function, by the way is another construction site. Currently, it is not possible to use the single-sided limits. We have that on our todo list.
In the example below you see that the return value from Maxima contains an expression with at(), which is removed on SMath postprocessing side even if optimization=none.
It is part of the test suite for Maxima in the interactive handbook or on github.
The limit function, by the way is another construction site. Currently, it is not possible to use the single-sided limits. We have that on our todo list.
Martin Kraska
Pre-configured portable distribution of SMath Studio: https://en.smath.info/wiki/SMath%20with%20Plugins.ashx
1 users liked this post
Alvaro Diaz Falconi 7/26/2023 10:13:00 PM
#4 Posted: 7/27/2023 4:47:43 AM
Hi Martin,
I think I can handle it plugin-side, I'll look into it ASAP.
I think I can handle it plugin-side, I'll look into it ASAP.
If you like my plugins please consider to support the program buying a license; for personal contributions to me: paypal.me/dcprojects
1 users liked this post
sergio 7/27/2023 8:19:00 AM
#5 Posted: 7/27/2023 1:11:57 PM
WroteYou see the problem if you evaluate input and return value for a particular f(x), the results aren't equal.
Pretty obvious that f(x):= sin(0.5*x)
evaluation t=0 will not be same.
Vectorize over t[0,0.125..pi]
#6 Posted: 7/28/2023 12:48:23 PM
By same token of diff at(2)
1. 1rst derivative of erf(x) is the integrand.
that you derive D(n,x) at will, range/vectorize.
2. Derivative of Gamma(x) from infinitesimal analysis.
1. 1rst derivative of erf(x) is the integrand.
that you derive D(n,x) at will, range/vectorize.
2. Derivative of Gamma(x) from infinitesimal analysis.
1 Pages (6 items)
-
New Posts
-
No New Posts