Without solving the question completely, what I can suggest is that use the information about symmetry of a and b with respect to mean to reduce the problem to find just one unknown, say c, and then find the value of c theoretically as a function of a quantile of the standard normal distribution, and then use the qnorm function in R to calculate the value.