function [output] = IPmotionblurfilter(ui,vi)
a = 0.1;
b = 0.1;
T = 1;

%[u,v] = ndgrid([0:ui-1],[0:vi-1]);
[u,v] = ndgrid([-ui/2 + 1:ui/2],[-vi/2 + 1:vi/2]);
puavb = pi * (u * a + v * b);
%size(puavb,1)
%size(puavb,2)

% Moet gecentreerd naar het midden, net als gaussian
% Bovendien uitzondering voor puavb = 0, de uitkomst moet dan 1 zijn.
% want dit is een sync functie

output =  T * (puavb).^(-1) .* sin(puavb) .* exp(-j * puavb);
%T * (puavb).^(-1) .* sin (puavb) .* exp(-j * puavb);
%output = uint8(output);
%output = uint8(output);
%output = output ./ puavb