function g = likelihood(model,x,y,params)

switch model
    case 'stocVolatility'
        b = params(1);
        g = exp(-y.^2 .* exp(-x) ./ (2*b^2) - x/2) ./ sqrt(2*pi*b^2);
        
    otherwise
        error('Unknown model');
end