%二維插值函數，輸入變量依次為：攻角數組，馬赫數數組，待插參數數組，插值點攻角值，插值點馬赫數值
%調用了一維插值函數oneinterp2

function F = twointerp(arg1,arg2,arg3,arg4,arg5)
[m n]=size(arg3);
%確定要用到哪二行馬赫數
if arg5<=arg2(1)%馬赫數小于或者等于最小給定馬赫數 
    i=1;
end    

if arg5>arg2(m)%馬赫數大于最大給定馬赫數
    i=m-1;
else
     i=1;
     while arg5>arg2(i+1)
          i=i+1;
     end    
end     

fy1=oneinterp2(arg1,arg3,arg4,i);
fy2=oneinterp2(arg1,arg3,arg4,i+1);
F=fy1+(fy1-fy2)*(arg5-arg2(i))/(arg2(i)-arg2(i+1));