Листинг 1

function L=lor(t,x);

global s r b;

L=[s*(x(2)-x(1));r*x(1)-x(2)-x(1).*x(3);-b*x(3)+x(1).*x(2)];

Листинг 2

global s r b;

s=10;

r=28;

b=8/3;

% задание начальных условий

x0=[0 0.2 0.6];

[T,X]=ode45('lor',[0 100],x0);

% lor – имя файла-функции, cодержащий определение

% функции; [0 100] – вектор, определяющий интервал

% интегрирования; x0 – вектор начальных условий

plot3(X(:,3),X(:,1),X(:,2))

Напишите программу вывода колебаний нелинейного математического маятника

dx/dt+w²sin(x)=0

Листинг 1

function F=mayat(t,y);

% определение функции, стоящей в правой части системы ОДУ (3)

global omega;

F=[y(2); - omega^2*sin(y(1))];