1% Copyright (C) 2009  VZLU Prague, a.s., Czech Republic
2%
3% Author: Jaroslav Hajek <highegg@gmail.com>
4%
5% This file is part of NLWing2.
6%
7% NLWing2 is free software; you can redistribute it and/or modify
8% it under the terms of the GNU General Public License as published by
9% the Free Software Foundation; either version 3 of the License, or
10% (at your option) any later version.
11%
12% This program is distributed in the hope that it will be useful,
13% but WITHOUT ANY WARRANTY; without even the implied warranty of
14% MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
15% GNU General Public License for more details.
16%
17% You should have received a copy of the GNU General Public License
18% along with this software; see the file COPYING.  If not, see
19% <http://www.gnu.org/licenses/>.
20%
21
22% -*- texinfo -*-
23% @deftypefn{Function File} plotplanform (ac, npan = 80)
24% @end deftypefn
25
26function plotplanform (wing)
27  zac = wing.zac;
28  xac = wing.xac;
29  cac = wing.cac;
30
31  x = [zac; flipud(zac)];
32  y = [xac - 1/4*cac; flipud(xac + 3/4*cac)];
33
34  if (wing.sym)
35    x = [x; -x];
36    y = [y; y];
37  else
38    x(end+1,:) = x(1,:);
39    y(end+1,:) = y(1,:);
40  endif
41
42  plot (x, y);
43  r = 0.05 * max(abs (x));
44  lims = [min(x) - r, max(x) + r];
45  set (gca, "xlim", lims);
46  set (gca, "ylim", lims);
47
48endfunction
49