--- /dev/null
+function [ftot]=FrequencySpace1(mf, pixelsize)
+
+% /!\ pixel size in meter, to get rid of "le"
+fsamplex = 1/(pixelsize);
+
+%[ftot,gtot] = meshgrid([0:fsamplex/mf:fsamplex/2 -fsamplex/2+fsamplex/mf:fsamplex/mf:-fsamplex/mf],[0:fsampley/nf:fsampley/2 -fsampley/2+fsampley/nf:fsampley/nf:-fsampley/nf]); %originale version
+ftot = meshgrid([0:fsamplex/mf:fsamplex/2 -fsamplex/2+fsamplex/mf:fsamplex/mf:-fsamplex/mf], 1);
+end