使用runge kutta 4交替集成rossler和lorenz吸引子的代码

Question

我需要一个matlab代码，通过引入耦合来使用RK4集成rossler和lorenz吸引子。

我写了罗斯勒吸引子和劳伦兹作为函数。然后，我在劳伦兹吸引子中引入耦合g，以便使罗斯勒成为推动者。这是我的代码化

a= 0.2; b= 0.2; c = 5.7;g =0;
% initial condition
xo  = [0.75;0.75;1];
yo  = [0.75;0.75;1];
zo  = [0.75;0.75;1];
% computing the trajectory
dt = 0.01;
tspan = (dt:dt:1000);
xinput = xo;
yinput = yo;
zinput = zo;
X = zeros(3,tspan(end)/dt);
Y = zeros(3,tspan(end)/dt);
Z = zeros(3,tspan(end)/dt);
 for i = 1: tspan(end)/dt
 time = i*dt;
 xoutput = rk4angelstepb(@rosslerb,dt,xinput,a,b,c);
 youtput = rk4angelstepb(@lorenzbb,dt,yinput,x,y,sigma,beta,rho,g);
 zoutput = rk4angelstepb(@lorenzbbb,dt,zinput,x,z,sigma,beta,rho,g);

 X(:,i) = xoutput;
 Y(:,i) = xoutput;
 Z(:,i) = xoutput;
 xinput = xoutput;
 yinput = youtput;
 zinput = zoutput;
 end
plot3(X(1,:),X(2,:),X(3,:),'b')

 for i = 1: tspan(end)/dt
 time = i*dt;
 xoutput = rk4angelstepb(@rosslerb,dt,xinput,a,b,c);
 youtput = rk4angelstepb(@lorenzbb,dt,yinput,x,y,sigma,beta,rho,g);
 zoutput = rk4angelstepb(@lorenzbbb,dt,zinput,x,z,sigma,beta,rho,g);

我希望吸引子聚在一起，但会收到此错误消息 未定义的函数或变量“ y”。

Error in sil (line 18)
 youtput = rk4angelstepb(@lorenzbb,dt,yinput,y,sigma,beta,rho,g);