a = 1;
s1 = .1;
s2 = .1;
s3= .03;

k = 0.001:0.01:20;
Y = [sin(a*k)./(k.^.5);cos(a*k)./(k.^.5);k];

load KFdata;

figure(1);clf;
rotate3d on;
plot3(X(1,:),X(2,:),X(3,:),'y.','linewidth',2);hold on;
plot3(Y(1,:),Y(2,:),Y(3,:),'m','linewidth',1);
drawnow;