%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%?Define?Parameters?（Node?amount）
n=100;
%%?Initial?Varibles
x=linspace（01n）;
ya=zeros（size（x））;
yb=zeros（size（x））;
z=ones（size（x））;
%%?main?loop?to?compute?Airfoil‘s?geometry?by?using???Mapping?????
????for?i=1:1:n
???????ya（i）=ya（i）+IterFunction（x（i））;
???????yb（i）=-ya（i）;
????end
????figure（‘Color‘[1?1?1]）
????plot（xya‘*-g‘‘linewidth‘2）;
????hold?on;
????plot（xyb‘+-r‘‘linewidth‘2）;
????hold?off;
????axis?equal;
????grid?on;
????title（‘NACA001