function coefficients = load_coefficients(filename)
    delimiterIn = ' ';
    headerlinesIn = 3;
    data = importdata(filename, delimiterIn, headerlinesIn);
    coefficients = data.data;
end

function amplitude = get_amplitude(coefs, control_param1, control_param2)
    acc = 0;
    for i = 1:4
        for j = 1:4
            acc += a(i, j) * control_param1 ^ (i-1) * control_param2 ^ (j-1);

    amplitude = 10 ^ acc
end
