Решение системы ОДУ

[Гость]

Помогите пожалуйста!
Необходимо построить для данной системы ОДУ первое диффернциальное приближение, а затем с помощью этого приближения численно ее решить.
[tex]\varphi''' + 3\varphi \varphi'' - 2\varphi'^2 + \theta = 0, \\
\theta''+ 3 Pr \varphi \theta' = 0.[/tex]
Граничные условия:
[tex]\varphi(0) = \varphi'(0) = 0, \\
\theta(0) = 1.[/tex]

Ранее эту задачу я решил с помощью метода внешних и внутренних разложений используя численный метод Рунге-Кутты и получил графики во вложении.

[Prorab]

Первое дифференциальное приближение (конечно-разностная схема):