Runge-Kutta program for system of ODE in Mathcad

[vnovikov]

Hello!

So I have a program which finds solution for single ODE using Runga-Kutta method, but I need somehow to rewrite it to make it solve the system of two ODE. My question is, how I can make it done?

[-MFra-]

Tk(t) and Rk(t) are undefined. You can't use Runge-Kutta algorithm.

[vnovikov]

offcourse I can, there 2 equations and 2 unknowns, Tk(tau) and R(tau) and I'm trying to find them. Actualy I already did it, by using ODEsovle, but my goal is to do it with own program, I don't want to solve it entirely via runge kutta, I want to solve it with adams bashforth method, but it requires 4 starting points, when I have just two (initial conditions) so I need to start with runge kutta method to get additional points and then continue with adams-bashforth method.

[-MFra-]

[vnovikov]

my program might have some synthax issues, could you explain a little more what have you done? Tau is like time for this ode, I need to find how temperature and radius of water droplet changes with time.

[-MFra-]

If the arguments of D1 and D2 are three, the arguments of f and g must be three as well, right?

However, I do not know the physical problem that you are trying to solve.

It seems to me that you have tackled the problem in the wrong way.

[ValeryOchkov · Accepted solution]

See please

[vnovikov]

Спасибо, этот вариант выглядит правильным, я попробую когда буду дома, я вижу вы из МЭИ? Я тоже, я студент первого курса магистратуры ИТАЭ. Эта программа очень мне пригодится для научной работы, у меря еще вопрос, а нет ли у вас подобной программы для метода Адамса-Башфорта? Этот метод реализуется в функции Odesolve в Mathcad.

[ValeryOchkov]

Vladimir Novikov написал(а):

Спасибо, этот вариант выглядит правильным, я попробую когда буду дома, я вижу вы из МЭИ? Я тоже, я студент первого курса магистратуры ИТАЭ. Эта программа очень мне пригодится для научной работы, у меря еще вопрос, а нет ли у вас подобной программы для метода Адамса-Башфорта? Этот метод реализуется в функции Odesolve в Mathcad.

Я рад, что могу помочь вам. С какой вы кафедры?

У меня есть ещё программа этого метода, но с переменным шагом.

В библиотеке МЭИ можно взять книгу с этюдом

Этюд 9. Движения планет: расчет и визуализация в среде Mathcad или Часы Кеплера

[vnovikov]

Я с кафедры Низких Температур, а вы? Я что то не с мог найти метод Рунге-Кутта с переменным шагом в тех материалах, которые вы прислали, можете прислать отдельно? Вообще мне нужен метод Рунге-Кутта, для "разгона" метода Адамса-Башфорта, потому что методу Адамса-Башфорта необходимо четыре начальных точки, в то время как Рунге-Кутта всего две, именно поэтому я хочу использовать алгоритм Рунге-Кутта, но моя основная цель написать программу подобная которой реализуется в Given-Odesolve блоке, может вы знаете какую-нибудь литературу по поводу метода Адамса-Башфорта?

[ValeryOchkov]

TOT

[VladimirN]

Vladimir Novikov wrote:

Я с кафедры Низких Температур, а вы?...

Персоналии