ｆｘ - プログラム倉庫

class fou{
public static void main(String args[]){
int pulse[][]=new int[101][101];
double a[][]=new double[101][101];
double b[][]=new double[101][101];
double alpha[][]=new double[101][101];
double beta[][]=new double[101][101];
double gamma[][]=new double[101][101];
double theta[][]=new double[101][101];
int x,y,sx,sy;
double h;
double z;
double px,py;
int m,n;
for (x=0;x<101;x++){
for (y=0;y<101;y++){
pulse[x][y]=0;
if (Math.random()>0.5)pulse[x][y]=1;
}
}

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

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

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

for (m=1;m<101;m++){
for (n=0;n<101;n++){
z=0;
for (y=-100;y<101;y++){
py=y*h;
sy=y;
if (sy<0)sy=-sy;
z=z+a[sy][n]*Math.cos(n*py)*h;
}
alpha[m][n]=z/Math.PI;
}
}

for (n=0;n<101;n++){
z=0;
for (y=-100;y<101;y++){
py=y*h;
sy=y;
if (sy<0)sy=-sy;
z=z+a[sy][n]*h;
}
alpha[0][n]=z/Math.PI;
beta[0][n]=0;
}

for (m=1;m<101;m++){
for (n=0;n<101;n++){
z=0;
for (y=-100;y<101;y++){
py=y*h;
sy=y;
if (sy<0)sy=-sy;
z=z+a[sy][n]*Math.sin(n*py)*h;
}
beta[m][n]=z/Math.PI;
}
}

for (m=1;m<101;m++){
for (n=0;n<101;n++){
z=0;
for (y=-100;y<101;y++){
py=y*h;
sy=y;
if (sy<0)sy=-sy;
z=z+b[sy][n]*Math.cos(n*py)*h;
}
theta[m][n]=z/Math.PI;
}
}

for (n=0;n<101;n++){
z=0;
for (y=-100;y<101;y++){
py=y*h;
sy=y;
if (sy<0)sy=-sy;
z=z+b[sy][n]*h;
}
theta[0][n]=z/Math.PI;
gamma[0][n]=0;
}

for (m=1;m<101;m++){
for (n=0;n<101;n++){
z=0;
for (y=-100;y<101;y++){
py=y*h;
sy=y;
if (sy<0)sy=-sy;
z=z+b[sy][n]*Math.sin(n*py)*h;
}
gamma[m][n]=z/Math.PI;
}
}

}
}