NACA0012翼型“C”型网格生成并解拉普拉斯方程

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