Fitting a probit model to 200 observations of 3 standard normal covariates

Question

Background and Task: Consider a random sample of size n with a binary outcome Y_i. Assume Y_i ~ Bern(pi_i). Assume the probit link function pi_i=Phi(X_i^T beta)