Error

 7 Pages«<567
 Previous Topic Next Topic
 Razonar #121 Posted : 10 December 2021 18:35:23(UTC) Rank: Advanced MemberGroups: Registered Joined: 28/08/2014(UTC)Posts: 998Was thanked: 548 time(s) in 357 post(s) Hi Viacheslav. Thanks again. I guess that with this new feature SMath must to be very close of version 1.0. Just a question about derivatives. Yours ode solvers returns n-1 columns for the n-degree differential equation, so: which is the better way to implement y (x) : diff(y(x),x,n) ? Because SMath diff doesn't recognize lspline. I think in 3 or 4 ways: return n columns for the differential equation of degree n, return functions y'(x), y''(x), ... and y(x) as you do for y(x), a function like lspline(3) but for the derivative (something like dlspline) since cubic splines have information about the derivative or a diff function like matlab's, which simply calculates the difference in the values of a vector, for each column of a matrix.But the truth is, I don't know which method could be the best, nor do I have any idea what is the best way to obtain the function derived from the result of any of your many solvers at this time.And one last question: could it be possible to implement the graph of parametric functions in the XY Plot plugin? Something like ParamtricPlot (x (t), y (t), a, b ) and that takes the number of points from the plugin configuration, or ParamtricPlot (x (t), y (t), a, b, N ) for N subdivisions (N + 1 points) of the interval a < t < b.Best regards.Alvaro.
 uni #122 Posted : 10 December 2021 20:42:11(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Originally Posted by: Razonar And one last question: could it be possible to implement the graph of parametric functions in the XY Plot plugin? Something like ParamtricPlot (x (t), y (t), a, b ) and that takes the number of points from the plugin configuration, or ParamtricPlot (x (t), y (t), a, b, N ) for N subdivisions (N 1 points) of the interval a I think it is possible to mark a "vector" variable with a tilde "~": ~t. It might look like this. I need to think about whether this will work. Russia ☭ forever Viacheslav N. Mezentsev 1 user thanked uni for this useful post. on 11/12/2021(UTC)
 uni #124 Posted : 13 December 2021 13:52:46(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Plugin updated.Changes:- added support for ODE systems with many variables;- Adams() ode solver added (explicit Adams-Bashforth 5);- CWD and TOL variables added;- project restructured.I don't know what AB method is used in Mathcad, but it clearly uses an adaptive step, since the usual explicit AB method is unstable (the small stiffness of the task makes it so). Edited by user 13 December 2021 18:22:33(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev 2 users thanked uni for this useful post. on 13/12/2021(UTC),  on 14/12/2021(UTC)
 uni #125 Posted : 13 December 2021 21:48:49(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Plugin updated.Changes:- Adams variable added (allows you to choose a method for the Adams() solver): "AB1".."AB8", "AB5" - default;- CWD variable value changed.To calculate the first steps of the method, Runge-Kutta of the 4th order is used, so while the order is not higher than the 5th it makes sense.Edited by user 15 December 2021 15:50:36(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev 1 user thanked uni for this useful post. on 14/12/2021(UTC)
 uni #126 Posted : 15 December 2021 15:49:35(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Plugin updated.Changes:- rkfixed variable added (allows you to choose a method for the rkfixed() solver: "RK1".."RK6", "RK4" - default);- to calculate the first steps for the Adams() solver now Runge-Kutta of the 6th order is used.Edited by user 16 December 2021 11:16:47(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev 1 user thanked uni for this useful post. on 15/12/2021(UTC)
 AleksMoisei #127 Posted : 17 December 2021 18:39:58(UTC) Rank: NewbieGroups: Registered Joined: 15/05/2020(UTC)Posts: 7Location: Voronezh Скажите пожалуйста, можно ли с помощью данного расширения производить символьные расчеты как в Маткаде или только это касается графиков. Пробую посчитать символьно передаточную функцию, задействовал Maple Tools, simplify выдал что-то не то что должен был по моим разумениям. Raschjot17.sm (25kb) downloaded 10 time(s).Edited by user 17 December 2021 18:41:52(UTC)  | Reason: Not specified
 uni #128 Posted : 17 December 2021 19:30:28(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Originally Posted by: AleksMoisei Скажите пожалуйста, можно ли с помощью данного расширения производить символьные расчеты как в Маткаде или только это касается графиков. Пробую посчитать символьно передаточную функцию, задействовал Maple Tools, simplify выдал что-то не то что должен был по моим разумениям.]Это расширение большей частью реализует численные аналоги некоторых функций из Mathcad. Символьные же вычисления можно попробовать выполнять при помощи MapleTools и функции maple(). Если она работает как-то не так, то можно посмотреть о "переговорах" между движками при помощи файла лога MapleTools.log, который находится в папке с плагином, если в файле конфигурации MapleTools.config включена опция отладки (True).Путь такой: %appdata%\SMath\extensions\plugins\32dfd679-8cfd-483a-b79a-19d5ea838750\Его нужно скопировать и вставить в адресную строку Проводника. Дальше выбрать подпапку с последний версией плагина, отредактировать файл MapleTools.config при закрытой программе и запустить её снова. После этого при старте будет создаваться файл лога MapleTools.log, в котором можно увидеть что именно передаёт функция maple() в символьный движок дополнения и как тот отвечает. Ну, а дальше уже искать ошибки.Следует также отметить, что синтаксис команды скорее всего будет отличаться, т.к. Mathcad имеет специальную обёртку, "скрывающую" некоторые дополнительные действия, чтобы перевести его запрос на понятный символьному движку "диалект". Поэтому нельзя в общем случае писать также как в Mathcad. Нужно понимать как работает Maple и быть знакомым с его командами, их форматами вызова.Edited by user 17 December 2021 19:39:14(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev
 Razonar #129 Posted : 17 December 2021 22:45:31(UTC) Rank: Advanced MemberGroups: Registered Joined: 28/08/2014(UTC)Posts: 998Was thanked: 548 time(s) in 357 post(s) Originally Posted by: AleksMoisei ... Maple Tools, simplify выдал что-то не то что должен был по моим разумениям. Raschjot17.sm (25kb) downloaded 10 time(s).Hi. Notice that SMath does not operate automatically expands matrix operations. Also, the maple's linalg package have some issues too, and sometimes needs to call evalm function, which means "eval as matrix". If D is a scalar I suppose that the correct result is (2), but if it is a 2x2 matrix I think it is (3).In both cases you cannot calculate A ^ (- 1) because the second column of A is zero.Could you comment what is the expression you expect for A? Also perhaps something from this post may be useful to you: https://en.smath.com/for...d-laplace-transform.aspxBest regards.Alvaro.
 uni #130 Posted : 21 December 2021 12:26:33(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Updated test for the ODE solvers. Section Math ODE precision (uni).sm (35kb) downloaded 13 time(s).Original document: Section Math ODE precision.sm (24kb) downloaded 49 time(s). Edited by user 22 December 2021 22:38:38(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev 2 users thanked uni for this useful post. on 21/12/2021(UTC),  on 22/12/2021(UTC)
 uni #131 Posted : 27 December 2021 23:04:34(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) I finally figured out how it works Mathcad Block. Spring pendulum.sm (49kb) downloaded 40 time(s). Edited by user 28 December 2021 01:19:04(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev 2 users thanked uni for this useful post. on 27/12/2021(UTC),  on 28/12/2021(UTC)
 uni #132 Posted : 24 April 2022 20:42:52(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Plugin updated.Changes:- MathcadBlock improved (doesn't perform substitutions);- File Access functions refactored;- project restructured.Here the variable x.min is defined above.Edited by user 24 April 2022 20:46:13(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev 1 user thanked uni for this useful post. on 25/04/2022(UTC)
 uni #133 Posted : 21 June 2022 18:44:53(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Plugin updated.Changes:- CWD and image processing functions fixed.P.S. Found a critical issue while trying to resolve this kind of system (x'(t) and x''(t) in one equation):The conversion of a complex system to a matrix form is not correct. A more complex algorithm is required. Only simple types of ODEs will be computed correctly.Edited by user 21 June 2022 20:48:22(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev 2 users thanked uni for this useful post. on 21/06/2022(UTC),  on 21/06/2022(UTC)
 Razonar #134 Posted : 22 June 2022 05:50:18(UTC) Rank: Advanced MemberGroups: Registered Joined: 28/08/2014(UTC)Posts: 998Was thanked: 548 time(s) in 357 post(s) Originally Posted by: uni ...P.S. Found a critical issue while trying to resolve this kind of system (x'(t) and x''(t) in one equation):...Hi Viacheslav. This is a way to handle some of coupled systems. Or at least that's what I think, I hope I didn't solve it wrong. odeCoupledExample.pdf (117kb) downloaded 4 time(s). odeCoupledExample.sm (88kb) downloaded 4 time(s).Best regards.Alvaro. 1 user thanked Razonar for this useful post. on 22/06/2022(UTC)
 uni #135 Posted : 22 June 2022 06:54:35(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Originally Posted by: Razonar I hope I didn't solve it wrong.These ODEs is called the Arenstorf Orbit. lsoda.Arenstorf Orbit.sm (16kb) downloaded 46 time(s). lsoda.Arenstorf Orbit.pdf (94kb) downloaded 48 time(s).The problem is that I use an universal algorithm that independently searches for derivatives in equations and arranges them in a matrix. I started with low order derivatives, so the equations were solved with respect to them first. It's wrong, but it worked in simple cases. When I redid the enumeration of derivatives of higher orders, I found that I did not take into account the case of different naming of variables. In this regard, an incorrect matrix is ​​obtained and the numerical algorithm cannot find a solution.I need to organize a double loop that first takes into account different variables, and then their derivatives.It is strange that I didn't do this earlier, because the rest of the code is designed for this.I had a very simple set of examples that I experimented with. As soon as I tried something more complicated - it immediately revealed the problem.Edited by user 22 June 2022 08:51:43(UTC)  | Reason: Not specified Russia ☭ forever Viacheslav N. Mezentsev 1 user thanked uni for this useful post. on 22/06/2022(UTC)
 Razonar #136 Posted : 22 June 2022 09:36:15(UTC) Rank: Advanced MemberGroups: Registered Joined: 28/08/2014(UTC)Posts: 998Was thanked: 548 time(s) in 357 post(s) Hi Viacheslav. I guess that there are not a general form for get the ode system from a coupled system with more than one variable derivative. For instance, this is the maple function for get the system if there are only one variable: https://www.maplesoft.co...ath=DEtools%2FconvertsysBut this other is for odes which "contains more than one function that is differentiated": https://www.maplesoft.co...aspx?path=convert%2FODEsThere it seems that the Abel equations have some very interesting properties. The relevant point here for me is that I know about nothing about that.But this method of isolating the second derivative for decoupling some useful equations is what I see that is useful for some usual system, for example, from lagrangian dynamics. This have some few other examples with coupled and not coupled systems with more than one dependent variable. Examples Lagrangian.sm (1,007kb) downloaded 8 time(s). Examples Lagrangian.pdf (2,187kb) downloaded 7 time(s).Best regards.Alvaro.Edited by user 22 June 2022 09:38:34(UTC)  | Reason: Not specified 1 user thanked Razonar for this useful post. on 22/06/2022(UTC)
 uni #137 Posted : 22 June 2022 11:09:06(UTC) Rank: Advanced MemberGroups: Registered, Advanced MemberJoined: 10/11/2010(UTC)Posts: 1,324Was thanked: 1157 time(s) in 668 post(s) Originally Posted by: Razonar I guess that there are not a general form for get the ode system from a coupled system with more than one variable derivative.Yes I agree, Alvaro. I meant that my algorithm is trying to convert the system of equations into a form that is understandable by numerical solvers. In our case, this is the "general form". I didn't know Maple had special features for this. Russia ☭ forever Viacheslav N. Mezentsev 1 user thanked uni for this useful post. on 22/06/2022(UTC)
 Razonar #138 Posted : 22 July 2022 08:25:20(UTC) Rank: Advanced MemberGroups: Registered Joined: 28/08/2014(UTC)Posts: 998Was thanked: 548 time(s) in 357 post(s) Originally Posted by: uni ... my algorithm is trying to convert the system of equations into a form that is understandable by numerical solvers ... Hi. This is a way for use rk-solvers for very hard coupled systems, calling al_nleqsol with a guess for the derivatives of high order at the origin. RKA.sm (51kb) downloaded 15 time(s). RKA.pdf (310kb) downloaded 12 time(s).Best regards. 2 users thanked Razonar for this useful post. on 22/07/2022(UTC),  on 22/07/2022(UTC)
 Users browsing this topic
 7 Pages«<567
Forum Jump
You cannot post new topics in this forum.
You cannot reply to topics in this forum.
You cannot delete your posts in this forum.
You cannot edit your posts in this forum.
You cannot create polls in this forum.
You cannot vote in polls in this forum.