function pickover(a,b,c,d,n,mutation)
if nargin==5
	mutation='standard';
end

if strcmp(mutation,'standard')
	x=zeros(1,n);
	x(1)=0.1;
	y=x;
	for i=2:n
		xx=sin(y(i-1)*b)+c*sin(x(i-1)*b);
		yy=sin(x(i-1)*a)+d*sin(y(i-1)*a);
		x(i)=xx;
		y(i)=yy;
	end
	plot(x,y,'r.')
elseif strcmp(mutation,'alpha')
	x=zeros(1,n);
	x(1)=0.1;
	y=x;
	for i=2:n
		xx=sin(y(i-1)*b)+sin(x(i-1)*b).^2+sin(x(i-1)*b).^3;
		yy=sin(x(i-1)*a)+sin(y(i-1)*a).^2+sin(y(i-1)*c).^3;
		x(i)=xx;
		y(i)=yy;
	end
	plot(x,y,'.')
elseif strcmp(mutation,'beta')
	x=zeros(1,n);
	x(1)=0.1;
	y=x;
	for i=2:n
		xx=sin(y(i-1)*b)+sin(x(i-1)*b).^2;
		yy=sin(x(i-1)*a)+sin(y(i-1)*a).^2;
		x(i)=xx;
		y(i)=yy;
	end
	plot(x,y,'.')
elseif strcmp(mutation,'gamma')
	x=zeros(1,n);
	x(1)=0.1;
	y=x;
	for i=2:n
		xx=abs(sin(y(i-1)*b))+sin(x(i-1)*b).^2;
		yy=abs(sin(x(i-1)*a))+sin(y(i-1)*a).^2;
		x(i)=xx;
		y(i)=yy;
	end
	plot(x,y,'.')
end

