clc; clear all; close all;
figure('Position',[0 0 800 800]) hold on; x = -1:0.005:1; sz = size(x); for n = 1:15 for i = 1:sz(2) CK(i) = cos(n*x(i)*acos(x(i))); end plot(x,CK,'k'); axis off; end
Slightly modifing the Carotid-Kundalini function
x = -1:0.005:1; sz = size(x); for n = 1:25 for i = 1:sz(2) for j = 1:sz(2) CK(i,j) = cos(n*(x(i))*acos(x(i)))+... cos(n*(x(j))*acos(x(j))); end end end figure('Position',[0 0 800 800]) imagesc(CK); colormap bone; axis off;