seinzumtode https://blog.hatena.ne.jp/seinzumtode/ how to code something https://seinzumtode.hatenadiary.jp/ function naca(n1,n2,n3,n4) m = 0.01*n1 p = 0.1*n2 t = 0.1*n3 + 0.01*n4 n = 100; xs = linspace(0,1,n); xus = []; xls = []; yus = []; yls = []; for idx=1:n x = xs(idx); yt = 5*t*(0.2969*sqrt(x)-0.1260*x-0.3516*x^2+0.2843*x^3-0.1015*x^4); if x <= p yc = m/p^2*(2*p*x-x^2); dyc_dx = 2*m/p^2*(p-x); else y… 190 <iframe src="https://hatenablog-parts.com/embed?url=https%3A%2F%2Fseinzumtode.hatenadiary.jp%2Fentry%2F2025%2F05%2F12%2F201841" title="NACA４桁翼型のプロット - how to code something" class="embed-card embed-blogcard" scrolling="no" frameborder="0" style="display: block; width: 100%; height: 190px; max-width: 500px; margin: 10px 0px;"></iframe> https://cdn-ak.f.st-hatena.com/images/fotolife/s/seinzumtode/20250512/20250512201826.png Hatena Blog https://hatena.blog 2025-05-12 20:18:41 rich https://seinzumtode.hatenadiary.jp/entry/2025/05/12/201841 1.0 100%