function interp2z = Interp2InNonGrid(a, b)

%step 1
db_yoldmatrix = [

   -2.9295   -2.6113   -2.7456   -2.6449   -2.9099   -2.7162   -2.7606
   -2.1809   -0.9963   -1.6443   -0.7211   -2.6984   -0.3680   -1.0798
   -1.8052   -0.4590   -0.4026   -0.2334   -1.1700   -0.0101   -0.5283
   -1.2936   -0.1847    0.1789    0.4070   -0.5078    0.8609   -0.3604
   -0.2077    0.0931    0.4788    1.0851    0.6172    1.7332    1.3598
    0.9687    0.4968    0.8432    1.7000    1.6077    2.8251    1.4674
    2.8634    2.9300    2.9622    2.7653    2.8462    2.9405    2.3923 ] ;

db_zoldmatrix = [
  0.00010   0.00340   -0.02990	-0.24500	-0.11000  -0.00430  0.00000 ; ...
  0.00070   0.04680   -0.59210	-4.75960	-2.10240  -0.06160  0.00040 ; ...
 -0.00880  -0.13010    1.85590	-0.72390	-0.27290   0.49960  0.01300 ; ...
 -0.03650  -1.33270   -1.65230	 0.98100	 2.93690   1.41220  0.03310 ; ...
 -0.01370  -0.48080    0.22890	 3.68860	 2.43380   0.58050  0.01250 ; ...
  0.00080   0.07970    2.09670	 5.85910	 2.20990   0.13280  0.00130 ; ...
  0.00000   0.00530    0.10990	 0.29990	 0.11070   0.00570  0.00000 ] ;

%step 2
db_newcoly =  [-2.5, -2.0, -1.5, -1.0, -0.5, 0.0, 0.5, 1.0, 1.5, 2.0, 2.2 ] ;
db_oldvectorx= [ -2.19736, -1.75720, -1.24647, 0.77271, 1.09999, 2.03543, 2.60028 ] ;

%step 3
db_newmatrixz = zeros(11, 7) ;

for i = 1:7
db_newmatrixz(:,i) = interp1( db_yoldmatrix(:,i), db_zoldmatrix(:,i), db_newcoly ) ;
db_newmatrixy(:,i) = db_newcoly 
end

for i= 1:11 ;
db_newmatrixx(i ,:) = db_oldvectorx ;
end

%step 4
interp2z = interp2(db_newmatrixx, db_newmatrixy, db_newmatrixz, a, b, 'cubic') ;