clear all;
clc;
x=0:0.005:1;
t=0.12;
m=0.2900;
k1=6.643;
yt=t/0.2*(0.2969*sqrt(x)-0.126*x-0.3516*x.^2+0.2843*x.^3-0.1015*x.^4);
t1=x(find(x<0.2900));
yc1=1.0/6.0*k1*(t1.^3-3*m*t1.^2+m^2*(3-m)*t1);
t2=x(find(x>=0.2900));
yc2=1.0/6.0*k1*m^3*(1-t2);
yc=[yc1 yc2];
y1=yt+yc;
y2=yc-yt;
plot(x,yc,'r',x,y1,'g',x,y2,'b');
hold on
axis equal
title('NACA 24012');
xlabel('chord');