ふぉうおｐ - プログラム倉庫

class fou{
public static void main(String args[]){
double data[]=new double[101];
double delta;
double alpha[]=new double[21];
int m,n,s;

for (n=0;n<101;n++){
data[n]=50;
}

delta=Math.PI/100;


for (m=0;m<20;m++){
alpha[m]=seekalpha(m,data);
}


for(s=1;s<40;s++){
System.out.println(seekfou(alpha,s*delta));
}

}

public static double seekalpha(int m,double[] data){
int n,nx,n1;
double delta,x1,y1,z;
delta=Math.PI/2000;
z=0;
for (nx=-1999;nx<2000;nx++){
n=nx;
if (n<0)n=-n;
n1=(int)(n/20);
x1=nx*delta;  
y1=data[n1]*Math.cos(m*x1);  
z=z+delta*y1;
}
z=z/Math.PI;
return z;
}


public static double seekfou(double[] alpha,double x){
double z;
int m;
z=alpha[0]/2;
for (m=1;m<20;m++){
z=z+alpha[m]*Math.cos(m*x);
}
return z;
}


}