function [k2,k3,k4,k5,k6]=decaygen(n)
[x,y]=meshgrid(1:n,1:n);
k2=(x.^2+y.^2);
k3=(x.^2+y.^2).^1.5;
k4=(x.^2+y.^2).^2;
k5=(x.^2+y.^2).^2.5;
k6=(x.^2+y.^2).^3;