Rkadapt Bug

Rkadapt Bug - Rkadapt Bug - Сообщения

#1 Опубликовано: 03.10.2016 02:42:28
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Page9 Scrap RK.sm (41 КиБ) скачан 61 раз(а).
1 пользователям понравился этот пост
Radovan Omorjan 03.10.2016 04:28:00
#2 Опубликовано: 03.10.2016 04:32:44
Radovan Omorjan

Radovan Omorjan

325 сообщений из 2052 понравились пользователям.

Группа: Moderator

Hmm...Interesting example, but it does not work for me in Mcad 15 either?

Regards,
Radovan
oderkadapt.png[attach]
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
#3 Опубликовано: 03.10.2016 04:57:56
Вячеслав Мезенцев

Вячеслав Мезенцев

1402 сообщений из 1708 понравились пользователям.

Группа: Moderator

I don't know how to convert implicit form into explicit. If you show me the rules for the numeric case I'll add such feature.
Russia ☭ forever, Viacheslav N. Mezentsev
#4 Опубликовано: 03.10.2016 14:31:25
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Viacheslav, Radovan: as simple as that
add a dummy 0 to the vector of RK solution
to meet the length of the IC !!!

Forum RK Bug Solved.gif
#5 Опубликовано: 03.10.2016 16:00:44
Radovan Omorjan

Radovan Omorjan

325 сообщений из 2052 понравились пользователям.

Группа: Moderator

Wrote

Viacheslav, Radovan: as simple as that
add a dummy 0 to the vector of RK solution
to meet the length of the IC !!!


Ah yes, but this is the standard "trick" in Mathcad when you want to give some parameter to the ODE. The derivative of that parameter is zero - like your dummy second value.
Here, if I understood this well, you might use SMath in a better way - giving parameters as third, fourth,...argument (not possible in Mathcad).
Your example might be this way as well

rkadapt.PNG
When Sisyphus climbed to the top of a hill, they said: "Wrong boulder!"
1 пользователям понравился этот пост
#6 Опубликовано: 04.10.2016 05:45:23
Вячеслав Мезенцев

Вячеслав Мезенцев

1402 сообщений из 1708 понравились пользователям.

Группа: Moderator

Functions from ODE Solvers plugin use the rules from Mathcad help. If you have two initial conditions this means that you have the second order ODE and f(t,y) must be a vector with two rows.

Цитата

Arguments:

init is either a vector of n real initial values, where n is the number of unknowns, or a single scalar initial value, in the case of a single ODE.


Rewriting a Higher-Order ODE as a System of First-Order ODEs.pdf (142 КиБ) скачан 64 раз(а).
Russia ☭ forever, Viacheslav N. Mezentsev
#7 Опубликовано: 04.10.2016 11:13:35
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

Thanks Viacheslav,

Very useful for those not familiar [never had Mathcad].
I have added the automatic RK vector extractor D(x,Y).

Jean

ODE Euler_rkf4 EDU.sm (31 КиБ) скачан 51 раз(а).
#8 Опубликовано: 04.10.2016 22:21:15
Jean Giraud

Jean Giraud

983 сообщений из 6866 понравились пользователям.

Группа: User

...just change variables for new names.
or copy/paste the section in blank sheet.
Smath is very reactive wrt variables.

RK Vector.gif
  • Новые сообщения Новые сообщения
  • Нет новых сообщений Нет новых сообщений