Модель турбулентности Лоренц

Построить модель динамической системы турбулентности (Лоренц, 1963г.) по следующим детерминированным дифференциальным уравнениям.
(dx(t))/dt=s?(y(t)-x(t)), (1)
(dy(t))/dt=r?x(t)-y(t)-x(t)z(t), (2)
(dz(t))/dt=x(t)y(t)-b?z(t). (3)
Параметры модели:s=10;r=28; b=8/3.
Шаг dtи время моделирования Tвыбрать самостоятельно.
На графиках отобразить фазовые портреты по всем функциям (зависимости x(t) от y(t),x(t) от z(t) иy(t) от z(t)) и графики зависимостей всех функций от времени по отдельности.