Long before Einstein, Reimann found that distance on any curved (or not) space could be expressed by a metric. Here are some examples:
ds^2=dx^2+dy^2
ds^2 = dr^2 +r^2 d\theta ^2
Hopefully, these two metrics look familiar to you Professor. They both represent flat space, one is in cartesian coordinates and one is in polar coordinates. It is noteworthy that although the two metrics look different, they both represent the same thing. Here is a metric that expresses distance on the surface of a sphere:
ds^2 = R^2d\phi ^2 + R^2 Sin^2(\phi) d\theta ^2
This surface now has a curvature (the first two metrics haves zero curvature they are flat). It is certainly not obvious by looking that this metric belongs to a curved surface where as the flat metric in polar coordinates does not. How can we tell the curvature of a space using the metric? Reimann and Gauss found that all you have to do is take a lot of derivatives!
Here is Reimann's general form of the metric:
ds^2 = \sum_{\mu,\nu} g_{\mu,\nu} dx^{\mu}dx^{\nu}
Hopefully you can see that this is:
ds^2 = g_{0,0} (dx^0)^2 + g_{1,0}dx^1 dx^0 +...
Notice that the metric coeffecients g_{i,j} (there are 16 of them in 4-d space time because we have four coordinates x^i where i = 0,1,2,3) contain all the information, for the reason we put them in 4x4 matrix called the metric tensor. Now like I said, once we have the metric tensor, we are only a few thousand derivatives away from finding the curvature.
From the metric tensor we define the affine connection:
\Gamma_{\sigma,\mu,\nu} = \frac{1}{2}\sum_{\mu,\nu} g^{\mu,\nu}(\partial _{\nu} g_{\mu,\sigma}-\partial_{\sigma}g_{\mu,\nu} + \partial_{\mu}g_{\sigma,\nu})
As you can see, any particular affine connection (Ex. \Gamma_{1,1,1}) contains 48 terms. So we already have 64 affine connections with 48 terms each.
From the affine connection we find the rank 4 reimann curvature tensor:
Reimann_{\lambda,\mu,\nu,\kappa} = \frac{1}{2}(\partial _{\kappa,\mu} g_{\lambda,\nu} +\partial_{\kappa,\lambda}g_{\mu,\nu} +\partial_{\nu,\lambda} g_{\kappa,\mu} + \sum _{\sigma,\epsilon} g_{\sigma,\epsilon}(\Gamma_{\epsilon, \nu, \kappa}\Gamma{\sigma,\mu,\kappa}-\Gamma_{\epsilon,\kappa,\lambda}\Gamma_{\sigma,\mu,\nu}
For General Relativity we contract this to the Ricci Tensor:
R_{\mu,\kappa} = \sum_{\nu,\lambda} g^{\lambda,\nu} Reimann_{\lambda, \nu,\mu,\kappa}
And the Ricci Scalar:
R = \sum_{\nu,\lambda} g^{\lambda,\nu} R_{\lambda, \nu}
In the most general case, computation of the Ricci Scalar involves over 15,000 terms. But at least now we have reached a point where we can express the EFE:
R_{\mu,\nu} -\frac{1}{2} g_{\mu,\nu}R = \frac{8 \pi G}{c^2} T_{\mu,\nu}
You should be able to interpret the left side of the equation as the curvature of space time. The right side is called the stress energy tensor, and it is basically energy density (and momentum densities in different directions).
The dependent variables in the system of PDEs are the metric coefficients, the independent variables are the four coordinates of spacetime; the components of the stress energy tensor are the inhomogeneous terms. A straightfoward method of approach would be to choose a stress energy tensor and solve the system of PDEs for the metric coefficients; but, in all but the simplest cases, the equations are so ultra-complex that nothing on Earth is capable of this method. A more much more tractable approach is to input a particular metric and stress energy tensor, and solve the EFE for the time evolution of that metric.