    function r = log_resistance_model(X,t)
    k = 0.01;
    nu = 1;
    r = X + log((1-k*t.^nu));
    end