function [sample] = initializeSample(args)

arguments
    args.radius (1,1) double {mustBeNumeric, mustBePositive} = 0.1 % [m]
    args.height (1,1) double {mustBeNumeric, mustBePositive} = 0.3 % [m]
    args.mass   (1,1) double {mustBeNumeric, mustBePositive} = 50 % [kg]
    args.freq   (1,1) double {mustBeNumeric, mustBePositive} = 100 % [Hz]
    args.offset (1,1) double {mustBeNumeric} = 0 % [m]
    args.x0     (1,1) double {mustBeNumeric} = 0 % [m]
    args.y0     (1,1) double {mustBeNumeric} = 0 % [m]
    args.z0     (1,1) double {mustBeNumeric} = 0 % [m]
end

sample = struct();

sample.radius = args.radius; % [m]
sample.height = args.height; % [m]
sample.mass = args.mass; % [kg]
sample.offset = args.offset; % [m]

sample.k.x = sample.mass * (2*pi * args.freq)^2; % [N/m]
sample.k.y = sample.mass * (2*pi * args.freq)^2; % [N/m]
sample.k.z = sample.mass * (2*pi * args.freq)^2; % [N/m]

sample.c.x = 0.1*sqrt(sample.k.x*sample.mass); % [N/(m/s)]
sample.c.y = 0.1*sqrt(sample.k.y*sample.mass); % [N/(m/s)]
sample.c.z = 0.1*sqrt(sample.k.z*sample.mass); % [N/(m/s)]

sample.x0 = args.x0; % [m]
sample.y0 = args.y0; % [m]
sample.z0 = args.z0; % [m]

save('./mat/stages.mat', 'sample', '-append');