LK

fun3 = @(t,x) [-5*x(1)+4*abs(x(2)); -3*x(1)*x(2)+2*sin(x(1))-x(2)+x(3); -x(1)*x(2) + 10*x(2)*cos(x(1))-5*x(3)];
options3 = odeset('Mass',[1 0 0 ; 0 1 0; 0 0 1]);
czas = [0 2];
x0 = [ 10 -10 5];
figure
[t,xa]=ode45(fun3,czas,x0,options3)
plot(t,xa)
hold on