cancel
Showing results for 
Search instead for 
Did you mean: 
cancel
Showing results for 
Search instead for 
Did you mean: 

Community Tip - Have a PTC product question you need answered fast? Chances are someone has asked it before. Learn about the community search. X

Runge-Kutta program for system of ODE in Mathcad

vnovikov
1-Newbie

Runge-Kutta program for system of ODE in Mathcad

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?

1 ACCEPTED SOLUTION

Accepted Solutions

11 REPLIES 11
-MFra-
21-Topaz II
(To:vnovikov)

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

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-
21-Topaz II
(To:vnovikov)

Fourth-order Runge-Kutta method.jpg

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-
21-Topaz II
(To:vnovikov)

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.

See please

9-25.png

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

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

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

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

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

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

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

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

TOT

Vladimir Novikov wrote:

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

Персоналии

Top Tags