ふぉうｙ - プログラム倉庫

class fou{
public static void main(String args[]){
int s,sx,n;
double x,y,z;
double h;
double a[]=new double[101];
double b[]=new double[101];
int pulse[]=new int[101];

for (s=0;s<101;s++){
pulse[s]=0;
if (Math.random()>0.5)pulse[s]=1;
}

h=Math.PI/100;
for (n=1;n<101;n++){
z=0;
for (s=-100;s<101;s++){
x=s*h;
sx=s;
if (sx<0)sx=-sx;
z=z+pulse[sx]*Math.cos(n*x)*h;
}
a[n]=z/Math.PI;
}

for (n=1;n<101;n++){
z=0;
for (s=-100;s<101;s++){
x=s*h;
sx=s;
if (sx<0)sx=-sx;
z=z+pulse[sx]*Math.sin(n*x)*h;
}
b[n]=z/Math.PI;
}

z=0;
for (s=-100;s<101;s++){
x=s*h;
sx=s;
if (sx<0)sx=-sx;
z=z+pulse[sx]*h;
}
a[0]=0.5*z/Math.PI;

for (s=-100;s<101;s++){
x=s*h;
sx=s;
if (sx<0)sx=-sx;
System.out.print(pulse[sx]);
System.out.print(",");
System.out.println(fou(a,b,x));
}

}

public static double fou(double[] a,double[] b,double x){
double f1;
int n;
f1=a[0];
for (n=1;n<101;n++){
f1=f1+a[n]*Math.cos(n*x)+b[n]*Math.sin(n*x);
}
return f1;
}

}