DATA_SECTION
init_int n
init_vector y(1,n)
PARAMETER_SECTION
init_bounded_number mu(-10,10)
init_bounded_number sigma(0,10)
PROCEDURE_SECTION
nll=-sum(dnorm(y,mu,sigma,1))
REPORT_SECTION
report << "mu: " << mu << endl;
report << "sigma: " << sigma << endl;