100たいぴ - プログラム倉庫

class pro{

int mlong,nlong,nx;
double[][] a=new double[502][902];
double[] b=new double[502];
double[] x=new double[502];
double h,c1,y1,c2,y2;
int t,mx;
double[] c=new double[101];
double[] y=new double[101];
double[] th=new double[101];
int s;

public static void main(String[] args) {
pro test=new pro();
}

pro(){

mlong=501;
nlong=901;

h=0.1;

for(s=1;s<101;s++){
th[s]=s;
y[s]=0.6*th[s];
c[s]=y[s];
}

t=0;
while(t<500){

for(mx=0;mx<mlong+1;mx++){
for(nx=1;nx<nlong+1;nx++){
a[mx][nx]=0;
}
b[mx]=0;
}

for(s=1;s<101;s++){
a[0][4*s-3]=1/c[s];
a[0][4*s-2]=-1/c[s];
a[0][4*s-1]=-1/(th[s]-y[s]);
a[0][4*s]=1/(th[s]-y[s]);
}


for(s=1;s<101;s++){
a[501][4*s-3]=1;
a[501][4*s-2]=-1;
a[501][4*s-1]=-1;
a[501][4*s]=1;
}

for(s=1;s<101;s++){
a[5*s-4][4*s-3]=1;
a[5*s-3][4*s-2]=1;
a[5*s-2][4*s-1]=1;
a[5*s-1][4*s]=1;
a[5*s][4*s]=1;
}

for(s=1;s<502;s++){
a[s][400+s]=1;
}

for(s=1;s<101;s++){
b[5*s-4]=h;
b[5*s-3]=h;
b[5*s-2]=h;
b[5*s-1]=h;
b[5*s]=y[s];
}

lp sub=new lp();
sub.mlong=mlong;
sub.nlong=nlong;
sub.a=a;
sub.b=b;
sub.makedata();
x=sub.x;

System.out.println(t);

for(s=1;s<101;s++){
c[s]=c[s]+x[4*s-3]-x[4*s-2];
y[s]=y[s]+x[4*s-1]-x[4*s];
}

t=t+1;
}

for(s=1;s<101;s++){
System.out.println(c[s]);
}


}
}