Cfunction [t,y]=rossler(y10,y20,y30,y40); %%This function will take in three initial conditions for the rossler %%attractor differential equation and then plot the solutions based on the %%initial conditions and the values a, b, and c %%Representing Cos(theta2-theta1) as c1 % Sin(theta2-theta1) as s1 y1='[y(2);'; y2='s1*(y(4)^2-y(2)^2*c1)+9.8*(sin(y(3))*c1-2*sin(y(1)))/(2-(c1)^2);'; y3='y(4);'; y4='s1*(c1*y(4)^2+2*y(2)^2)-2*9.8*(sin(y(1)*c1-sin(y(3)))/((c1)^2-2))]'; F=strcat(y1,y2,y3,y4) c1='cos(y(3)-y(1))'; s1='sin(y(3)-y(1))'; F_a=regexprep(F,'c1',c1); F_b=regexprep(F_a,'s1',s1) markers=['--' '.' '*']; func=inline(F_b,'t','y') clf reset; cmap=colormap; %hold on; [t,y]=ode45(func,[0.01,100],[y10,y20,y30,y40]); plot(y(:,1),y(:,3)) %hold on; %subplot(2,2,1); %h=plot3(y(:,1),y(:,2),y(:,3),markers(i)); %color_index=i*10; %set(h,'color',cmap(color_index,:));