The intersection function will rarely exist as a "function". At best, it is an implicit function, thus a problem with solution locus defined only as the set of solutions to an equation. Your implicit function is nonlinear, and sufficiently so that no simple algebraic form will exist to extract the solution locus for your problem.
To plot that implicit function, you use a contour plot. In fact, it appears you already know that fact, so there is little I need to do here.
You can extract the points along that contour using contourc, without creating the contour plot itself. However, contour also returns that same information as it does the contour plot.
M1 = contour(X1,Y1,Z1,v)
M1 =
Columns 1 through 9
0.9 1.7 1.675 1.65 1.625 1.6 1.575 1.55 1.525
17 0.17037 0.17213 0.17411 0.17635 0.17889 0.18175 0.18499 0.18867
Columns 10 through 18
1.5 1.475 1.4641 1.45 1.425 1.4 1.375 1.35 1.3381
0.19286 0.19765 0.2 0.20358 0.2109 0.2196 0.23005 0.24279 0.25
Here, M1(1,1) is the value of the level curve, thus 0.9.
M1(2,1) is the number of points it found along that curve. Note that sometimes a contour may be represented as disjoint, non-intersecting sets. so you may have multiple curves that represent a single contour. So these first two numbers are necessary to describe the contour.
