{"html":"<iframe src=\"https://hatenablog-parts.com/embed?url=https%3A%2F%2Fseinzumtode.hatenadiary.jp%2Fentry%2F2025%2F05%2F12%2F201841\" title=\"NACA\uff14\u6841\u7ffc\u578b\u306e\u30d7\u30ed\u30c3\u30c8 - 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>","url":"https://seinzumtode.hatenadiary.jp/entry/2025/05/12/201841","type":"rich","version":"1.0","published":"2025-05-12 20:18:41","blog_url":"https://seinzumtode.hatenadiary.jp/","categories":[],"image_url":"https://cdn-ak.f.st-hatena.com/images/fotolife/s/seinzumtode/20250512/20250512201826.png","height":"190","width":"100%","provider_url":"https://hatena.blog","title":"NACA\uff14\u6841\u7ffc\u578b\u306e\u30d7\u30ed\u30c3\u30c8","description":"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\u2026","blog_title":"how to code something","author_name":"seinzumtode","provider_name":"Hatena Blog","author_url":"https://blog.hatena.ne.jp/seinzumtode/"}