{"id":1618,"date":"2021-10-25T14:40:26","date_gmt":"2021-10-25T05:40:26","guid":{"rendered":"http:\/\/www.miyamoto-lab.com\/?page_id=1618"},"modified":"2024-11-13T10:15:24","modified_gmt":"2024-11-13T01:15:24","slug":"32iahsdfgja","status":"publish","type":"page","link":"https:\/\/www.miyamoto-lab.com\/ja\/32iahsdfgja\/","title":{"rendered":"11. 11\u7ae0\u7df4\u7fd2\u554f\u984c\u306e\u7b54\u3048\u30fb3"},"content":{"rendered":"\n<h4 class=\"wp-block-heading\">\u7df4\u7fd2\u554f\u984c\uff12<\/h4>\n\n\n\n<div class=\"hcb_wrap\"><pre class=\"prism line-numbers lang-c\" data-lang=\"C\"><code>#include&lt;stdio.h&gt;\n#include&lt;math.h&gt;\n#define T 0.5     \/\/ calculation duration\n#define max 10    \/\/ the maximum number of i\n#define dt 0.01   \/\/ delta t\n#define dx 0.1    \/\/ delta x\n#define omega 1.2 \/\/ omega\n#define a 1.0     \/\/ a\n#define eps 0.00001\n#define pi 3.141592\ndouble  u[max+1][2]; \/\/u_i^n is written as u[i][n], while n=0 or 1\n\nint main(void){\n  double uN,time=0;\n  double err,r=a*dt\/dx\/dx; int i,j;\n  for(i=0;i&lt;=max;i++)\n    u[i][0]=0.0;u[i][1]=0.0;\n  for(i=0;i&lt;=max;i++){     \/\/ define the initial condition\n    u[i][0]=sin(pi*(double)i*dx);\n  }\n  u[0][1]=0; u[max][1]=0;  \/\/ put the boundary condition\n  while (T&gt;time){\n   do {\n     err=0.0;\n     for(i=1;i&lt;max;i++){\n       uN=u[i][1];\n       u[i][1]=(1-(double)omega)*uN\n                + (double)omega*(1\/(1+2*r)*u[i][0]\n                + r\/(1+2*r)*(u[i+1][1]+u[i-1][1]));\n       if (err&lt;fabs(u[i][1]-uN))\n         err=fabs(u[i][1]-uN);\n     }\n   }while(err&gt;eps);\n   printf(&quot;%3.3f\\t&quot;,time);\n   time=time+dt;\n   for(i=0;i&lt;=max;i++)\n     printf(&quot;%3.3f\\t&quot;,u[i][0]);\n   printf(&quot;\\n&quot;);\n   for(i=1;i &lt; max;i++)\n     u[i][0]=u[i][1];\n  }\n  return (0);\n}<\/code><\/pre><\/div>\n","protected":false},"excerpt":{"rendered":"<p>\u7df4\u7fd2\u554f\u984c\uff12<\/p>\n","protected":false},"author":3,"featured_media":0,"parent":0,"menu_order":0,"comment_status":"closed","ping_status":"closed","template":"","meta":{"_locale":"ja","_original_post":"http:\/\/www.miyamoto-lab.com\/?page_id=1618","footnotes":""},"class_list":["post-1618","page","type-page","status-publish","hentry","ja"],"_links":{"self":[{"href":"https:\/\/www.miyamoto-lab.com\/wp-json\/wp\/v2\/pages\/1618","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.miyamoto-lab.com\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/www.miyamoto-lab.com\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/www.miyamoto-lab.com\/wp-json\/wp\/v2\/users\/3"}],"replies":[{"embeddable":true,"href":"https:\/\/www.miyamoto-lab.com\/wp-json\/wp\/v2\/comments?post=1618"}],"version-history":[{"count":4,"href":"https:\/\/www.miyamoto-lab.com\/wp-json\/wp\/v2\/pages\/1618\/revisions"}],"predecessor-version":[{"id":1876,"href":"https:\/\/www.miyamoto-lab.com\/wp-json\/wp\/v2\/pages\/1618\/revisions\/1876"}],"wp:attachment":[{"href":"https:\/\/www.miyamoto-lab.com\/wp-json\/wp\/v2\/media?parent=1618"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}