MATLAB: Using Runge Kutta function to solve Rossler equations

chaotic dynamicschaotic systemsexplicit methodslorenz equationsrossler equationrunge kutta

I have created a runge kutta 4th order function by so

Best Answer

Write your f function as
x'(t)=-x(t)-y(t)
y'(t)=ax(t)+y(t)
z'(t)=b+z(t)(x(t)-c)
%where I fix parameters say c=5, b=3, a=0.2 Xˆ = (0, 0, 0.2)^T time interval [0,250] and step size h=0.01
f = @(t,y) [ -y(1)-y(2)
a*y(1)+y(2)
b+y(3)*(y(1)-c) ];