clear?all
clc
%物理域繪制
%?畫出左側和上側外界區域
x=-30:0.8:50;
for?i=1:length（x）
????if?i<=31
???????y（i）=（900-x（i）^2）^0.5;
????elseif?i<=length（x）
????????y（i）=30;
????end
end
x0=x;y0=y;

x22=0:0.0125:1;??%畫翼型（翼型單側一共80個點）
y22=0.1781*x22.^0.5-0.0756*x22-0.2122*x22.^2+0.1705*x22.^3-0.0609*x22.^4;

%翼型上的數據儲存（81個）
x_wall=[x22];
y_wall=[y22];
y_wall（:end）=0;
l_wall=length（y_wall）;
%畫翼面以后的網格
x_low=3.33:2.33:50;
y_low=zeros（size（x_low））;
l_low=length（y_low）;

%左側上側邊界數據儲存（暫時只存到上邊界結束）
x_left=[x0（1:l_wall）];
y_left=[y0（1:l_wall）];
l_left=length（y_left）;
%上側的后半部分
x_up=[x0（l_wall+1:l_wall+l_low-1）];
y_up=[y0（l_wall+1:l_wall+l_low-1）];
l_up=length（y_up）;

n=9;%插值層數

%儲存位置的大數組（暫時只存到上邊界結束）
X=zeros（n+1l_wall+l_low）;
Y=zeros（n+1l_wall+l_low）;

%給X和Y的第一層（也就是第一行）賦值（分成n層，實際有n+1