System of 2 stiff, nonlinear, second order differential equations. (A catalysis problem)

So I am trying to model a system of 2 equations.

I AM NOT ASKING YOU TO DO MY HOMEWORK FOR ME!!! haha Just push me in the right direction.

The equations are as follows:

1. d^2(y)/dx^2 + 2/x*dy/dx = O^2*y*exp(-2*g*(1/t-1))
2. d^2(y)/dx^2 + 2/x*dt/dx = -O^2*B*y*exp(-2*g*(1/t-1))

O,g, and B are constants with y,t, and x being the variables.

They model the profile of the dimensionless concentration y and the dimensionless temperature t along the dimensionless radius from the catalyst center x.

I know that I can break the second derivatives down to the derivative of an auxiliary variable that represents the first derivative of the root variable. My issue is that I have never really modeled something like this in MatLab and would like some pointers.

The ranges for the variables are:

0 < y < 1

? < t < 1

0 < x < 1

If I break apart the second order equations, I get:

1. Y = dy/dx
2. dY/dx = -2*Y/x + O^2*y*exp(-2*g*(1/t-1))
3. T = dt/dx
4. dT/dx = -(2*T/x + O^2*B*y*exp(-2*g*(1/t-1))

But this is where I get stuck. What do I do with these things?

I really appreciate any kind of help with this. This is not an urgent matter, but I would certainly like to learn how to do this stuff.

Thanks guys!