function [t,y]=rossler(y10,y20,y30,a,b,c);

%%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;

for i=1:3
    
[t,y]=ode45(func,[0.01,100],[y10*i,y20*i,y30*i]);

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,:));
grid on;
xlabel('x');ylabel('y');zlabel('z');
hold on;

subplot(2,2,2);
hold on;
k=plot(y(:,1),y(:,2),markers(i));
xlabel('x');ylabel('y');
%set(k,'color',cmap(color_index,:));

subplot(2,2,3);
hold on;
plot(y(:,2),y(:,3),markers(i));
xlabel('y');ylabel('z');

subplot(2,2,4);
hold on;
plot(y(:,1),y(:,3),markers(i));
xlabel('x');ylabel('z');


end