how can I solve two equations with two unknown variable using matlab?

point_x = 2;
point_y = 4;
length_tang = 1
r = 5;

syms x y
E = [(-x)*(point_x - x) + (r - y)*(point_y - y) == 0, (point_x - x)^2 + (point_y - y)^2 == length_tang^2];
S = solve(E,x,y)
S = 
  struct with fields:

      x: [2×1 sym]
      y: [2×1 sym]

S.x
ans =
   2
 6/5

S.y
ans =
    5
 17/5

Too slow? It came back virtually instantaneously.

timeit(@() solve(E,x,y))
ans =
      0.19394

0.2 seconds is too slow?

Note that it is probably possible for SOME sets of inputs, there is no real valued solution at all.