You are an expert in writing code in many different languages. Your goal is to perform code completion for the following code, keeping in mind the rest of the code and the file meta data.
    Metadata details:
    filename: plot_increment_pdf.m
    path of file: ./repos/OPEN_FPE_IFT/Functions_v5-MATLAB
    
    the code of the file until where you have to start completion: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This function plots the probability density function (PDF) of the velocity increments 
% at the scale $r=L$, $r=\lambda$ and  $r=\eta$. The colored dashed line correspond to Castaing fits 
% (form factor $\Lambda_r$ \cite{Castaing_1990}) and grey dashed line to Gaussian fits.
%
% Arguments IN
% data = 1D array of data to which you would like to apply the low pass filter
% m_data = mean of the data
% nbins = The number of bins in which you would like to divide your data
% Fs = Acquisition/Sampling Frequency in Hz
% int_L = Integral length scale in meters
% taylor_L = Taylors length scale in meters
% diss_scale = Kolmogorv lenght scale in meters
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function plot_increment_pdf(data,nbins,save_path,save_name,int_L,taylor_L,m_data,Fs,norm_ur,varargin)
set(groot, 'defaultAxesTickLabelInterpreter','latex'); 
set(groot, 'defaulttextinterpreter','latex');   
set(groot, 'defaultLegendInterpreter','latex');
set(groot, 'defaultAxesT
    
function plots the probability density function (PDF) of the velocity increments 
% at the scale $r=L$, $r=\lambda$ and  $r=\eta$. The colored dashed line correspond to Castaing fits 
% (form factor $\Lambda_r$ \cite{Castaing_1990}) and grey dashed line to Gaussian fits.
%
% Arguments IN
% data = 1D array of data to which you would like to apply the low pass filter
% m_data = mean of the data
% nbins = The number of bins in which you would like to divide your data
% Fs = Acquisition/Sampling Frequency in Hz
% int_L = Integral length scale in meters
% taylor_L = Taylors length scale in meters
% diss_scale = Kolmogorv lenght scale in meters
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
function plot_increment_pdf(data,nbins,save_path,save_name,int_L,taylor_L,m_data,Fs,norm_ur,varargin)
set(groot, 'defaultAxesTickLabelInterpreter','latex'); 
set(groot, 'defaulttextinterpreter','latex');   
set(groot, 'defaultLegendInterpreter','latex');
set(groot, 'defaultAxesTitle','latex');

tau1            = round((int_L/m_data)*Fs);
tau2            = round((taylor_L/m_data)*Fs);

[incr1,incr2]   = Increment(tau1,tau2,data);

incr1           = incr1.';
incr2           = incr2.';

if nargin >= 9 
    tau3                = round((varargin{1}/m_data)*Fs);
    [incr1,incr2,incr3] = Increment_double_condition(tau1,tau2,tau3,data);
    incr1               = incr1.';
    incr2               = incr2.';
    incr3               = incr3.';
end