(This article was first published on The power of R, and kindly contributed to R-bloggers)
Hey all of you,I got such an interesting problem, a friend of mine is modelling something using bayesian statistics and she got an equation system to solve, but she's stucked and she showed to me such a long script to ty to solve it with numeric approximation, anyone has any idea about how to solve it?

As you can see we're talking about the normal distribution function, so.... let's see if any of you has any idea about how to solve this.
To leave a comment for the author, please follow the link and comment on his blog: The power of R.
R-bloggers.com offers daily e-mail updates about R news and tutorials on topics such as: visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...

Zero Inflated Models and Generalized Linear Mixed Models with R.
Zuur, Saveliev, Ieno (2012).