function [F,r_buz]=BuzoriusCorrection(Fm,u,tau_c)

% Buzorius model 
nm=0.085;
alpha=7/8;
% tau_c=0.4;
% u=8;
z=18;

% u=[0:20];
r_buz=1./(1+(2.*pi.*nm.*tau_c.*u./z).^alpha);

% figure;plot(u,r,'--');

% Fm=F.*r_buz
F=Fm./r_buz;

