function [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)-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,10],[y10,y20,y30,y40]);

%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,:));
gr