function [xcoord,ycoord] = make_naca(naca_code,num_coord,length)

xtemp = linspace(0,1,num_coord);
xu = linspace(0,length,num_coord); xl=xu;
max_camber = str2num(naca_code(1))/100;
camber_pos = str2num(naca_code(2))/10;
p = round(camber_pos*num_coord);
c_line(1:p) = max_camber/camber_pos^2*(2*camber_pos*xtemp(1:p)-xtemp(1:p).^2);
c_line(p+1:num_coord) = max_camber/(1-camber_pos)^2*(1+2*camber_pos*(xtemp(p+1:num_coord)-1)-xtemp(p+1:num_coord).^2);
c_line = c_line*length;
thick = str2num(naca_code(3:4))/100*length;
thick_dist = thick*5*(.2969*xtemp.^.5-.1260*xtemp-.3516*xtemp.^2+.2843*xtemp.^3-.1015*xtemp.^4);
thick_dist = thick_dist - xtemp.^5*thick_dist(num_coord);
yu = c_line+thick_dist;
yl = c_line-thick_dist;

xcoord=[fliplr(xl),xu(2:num_coord)];
ycoord=[fliplr(yl),yu(2:num_coord)];