プログラマ専用SNS ミクプラ

コード:

#include <stdio.h>
#include <math.h>

int main(void)
{
	int i,t,ni,nt;
	double dx,dt,a,pri,vis,time,pi();
	double u1[1000];
	double u2[1000];
	double D[100000][100];


	a=1.00;
	nt=100000;
	ni=100;
	dt=0.001;
	dx=0.0005;
	pri=10.0;
	vis=0.000001;

	for(i=1;i<=1000;i++){
		u1(i)=0;
	}
		
		for(t=1;t<=nt;t++){
			u1(1)=a*sin(((2*pi())/pri)*t*dt);
			u1(ni)=0.0;

			for(i=2;i<=ni-1;i++){
				u2(i)=vis*((u1(i+1)+u1(i-1)-2.0*u1(i))/(dx*dx))*dt+u1(i);
			}

			for(i=2;i<=ni-1;i++){
				u1(i)=u2(i);
			}
			for(i=1,i<=ni;i++){
				D(t,i)=u1(i);
			}
			if(t % 100 == 0){
				time=t*dt;
				printf("%f %f",time,u1(20));
			}
		}

		t=10000;
		printf("距離,(ωt=0,ωt=π/2,ωt=π,ωt=3π/2");

		for(i=1;i<=ni;i++){
			printf(i*dx,D(t,i),D(t+2500,i),D(t+5000,i),D(t+7500,i);
		}

		return 0;
}

コード:

#include <stdio.h>
#include <math.h>
#define PI 6*asin(0.5)

int main(void)
{
	int i,t,ni,nt;
	double dx,dt,a,pri,vis,time;
	double u1[1000];
	double u2[1000];
	double D[100000][100];


	a=1.00;
	nt=100000;
	ni=100;
	dt=0.001;
	dx=0.0005;
	pri=10.0;
	vis=0.000001;

	for(i=1;i<=1000;i++){
		u1[i]=0;
	}
		
		for(t=1;t<=nt;t++){
			u1[1]=a*sin(((2*PI)/pri)*t*dt);
			u1[ni]=0.0;

			for(i=2;i<=ni-1;i++){
				u2[i]=vis*((u1[i+1]+u1[i-1]-2.0*u1[i])/(dx*dx))*dt+u1[i];
			}

			for(i=2;i<=ni-1;i++){
				u1[i]=u2[i];
			}
			for(i=1,i<=ni;i++;){
				D[t][i]=u1[i];
			}
			if(t % 100 == 0){
				time=t*dt;
				printf("%f %f",time,u1[20]);
			}
		}

		t=10000;
		printf("距離,(ωt=0,ωt=π/2,ωt=π,ωt=3π/2");

		for(i=1;i<=ni;i++){
			printf("%f %f %f %f %f",i*dx,D[t,i],D[t+2500,i],D[t+5000,i],D[t+7500,i]);
		}

		return 0;
}

コード:

#include <stdio.h>
#include <math.h>
#define PI 6*asin(0.5)

int main(void)
{
	int i,t,ni,nt;
	double dx,dt,a,pri,vis,time;
	double u1[1000];
	double u2[1000];
	double D[100000][100];


	a=1.00;
	nt=100000;
	ni=100;
	dt=0.001;
	dx=0.0005;
	pri=10.0;
	vis=0.000001;

	/*一次元配列u1[i]の中身を初期化*/
	for(i=0;i<1000;i++){
		u1[i]=0;

	}
	    roop1:
		for(t=0;t<=nt;t++){         　　　　　　
			u1[0]=a*sin(((2*PI)/pri)*(t+1)*dt);      /*初期条件*/
			u1[ni]=0;                           　　　　　 /*初期条件*/
			u1[1]=u1[0]+0.0005;
			u1[2]=u1[0]+0.001;
			for(i=0;i<ni-1;i++){    
				u2[i+1]=vis*((u1[i+2]+u1[i]-2.0*u1[i+1])/(dx*dx))*dt+u1[i+1];  /*微分を行う*/
			}

			for(i=0;i<ni-1;i++){
				u1[i+1]=u2[i+1];       /*次の計算に備えて変数を更新*/
			}
			for(i=0,i<ni;i++;){    
				D[t][i]=u1[i];     
			}
			if(t % 100 == 0){      /*t/100の余りが0となったとき(tが100の倍数の時)以下を行う*/
				time=t*dt;         /*時間の更新*/
				printf("%f %f",time,u1[20]);
			}
		}
		
		t=70000;
		printf("距離,(ωt=0,ωt=π/2,ωt=π,ωt=3π/2");

		for(i=0;i<=ni;i++){
			printf("%f\n %f\n %f\n %f\n %f\n",i*dx,D[t][i],D[t+2500][i],D[t+5000][i],D[t+7500][i]);
		}

		return 0;
}

コード:

for(i=0,i<ni;i++;)

コード:

    for(i=0;i<ni;i++)