xorig=load('check/validx.dat');
yorig=load('check/validy.dat');

x=load('data/validx.dat');
y=load('data/validy.dat');

nxdiff = abs((x - xorig) ./ xorig);
nydiff = abs((y - yorig) ./ yorig);

xdiff = abs((x - xorig));
ydiff = abs((y - yorig));

maximum = [max(max(xdiff(:)), max(ydiff(:))), max(max(nxdiff(:)), max(nydiff(:)))]
average = [max(norm(xdiff(:))/sqrt(length(xdiff(:))), norm(ydiff(:))/sqrt(length(ydiff(:)))), max(norm(nxdiff(:))/sqrt(length(nxdiff(:))), norm(nydiff(:))/sqrt(length(nydiff(:))))]

