clear?all
close?all
clc

load?naca1.dat
X=naca1（:1）;
Y=naca1（:2）;
RE=6e6;
MACH=0;
alpha=3;
tic
[p]=xfoil（XYalphaREMACH）;
toc
CL=p.cl;
CD=p.cd;
CM=p.cm;
figure（1）
plot（p.xp.y）
axis?equal
figure（2）
if?（RE~=0?||?MACH~=0）
????plot（p.xp.y-2p.x-p.cpip.x-p.cpv）
????xlabel（‘x‘）
????ylabel（‘-C_p‘）
????legend（‘airfoil‘‘cpi‘‘cpv‘）
????xlim（[-0.1?2]）
????grid?on
else
????plot（p.xp.y-2p.x-p.cpi）
????xlabel（‘x‘）
????ylabel（‘-C_p‘）
????legend（‘airfoil‘‘cpi‘）
????xlim（[-0.1?2]）
????grid?on
end

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????500??2015-03-17?09:39??main.m
?????文件????????4320??2014-05-20?10:19??naca1.dat
?????文件????????4320??2014-06-15?15:49??NACA3320.dat
?????文件??????319168??2015-03-17?10:11??xfoil.m
?????文件????????1313??2015-03-17?09:11??license.txt