conv2(Dx,h,<span class="string">'same'</span>);
Dy = conv2(Dy,h,<span class="string">'same'</span>);

Dx = alpha*(Dx./norm(Dx));
Dy = alpha*(Dy./norm(Dy));
<span class="keyword">end</span>

<span class="keyword">function</span>[X2,y2] = create_distorted_dataset(X,y,sig,alpha,n,m)
N = (n+1)*size(X,1);
X2 = zeros(N,size(X,2));
y2 = zeros(N,1);
index = 1;
<span class="keyword">for</span> i=1:size(X,1)
    fprintf(<span class="string">'distorting data point: %d\r'</span>,i);
    <span class="keyword">for</span>