%
%This function is run by  dd
%
%  Last change:  10/2/2005 - Florian Knorn

function p = sf_reg_model(param,k)

% param(1) = scaling
% param(2) = kappa
% param(3) = gamma

p = param(1) * (k+param(3)).^( -param(2) );