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

F='[-(y(2)+y(3));y(1)+a*y(2);b+y(3)*(y(1)-c)]';

a_eq=num2str(a);
b_eq=num2str(b);
c_eq=num2str(c);

F_a=regexprep(F,'a',a_eq);
F_b=regexprep(F_a,'b',b_eq);
F_c=regexprep(F_b,'c',c_eq);

markers=['--' '.' '*'];

func=inline(F_c,'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