var x, y;
varexo e;
parameters alpha, beta;
alpha = 3;
beta = 3;

external_function(name = betainc, nargs = 3);
model;
y = exp(e)/20;
x = betainc(min(y, 1), alpha, beta);
end;

shocks;
var e; stderr 1.0;
end;

steady_state_model;
y = 1/20;
x = betainc(1/20, alpha, beta);
end;

steady;