Solving 3x1 matrix, two variables

I am trying to solve a 3x1 matrix which looks like follows:

C = [(2*cos(theta2)*sin(theta3))/25; (2*cos(theta2)*cos(theta3))/25 ; (2*sin(theta2))/25 + 81/1000];

The goal is get a answer for theta2 and theta3 for which C is equal to C = [0.0862;0.2370;0.1136];

I tried to it like this:

syms = theta2 theta3

AC = 0.08

OA = 0.081

C = [AC*cos(theta2)*sin(theta3); AC*cos(theta2)*cos(theta3); OA + AC*sin(theta2)];

solutions1 = solve(C == [0.0862;0.2370;0.1136])

however, doing this it get the following:

solutions1 =

struct with fields:

theta2: [0×1 sym], in which val =

Empty sym: 0-by-1

theta3: [0×1 sym], in which val =

Empty sym: 0-by-1

Can you please help me to solve this problem.

Kind regards.