どんｂ - プログラム倉庫

class pro{

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

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

pro(){

mlong=41;
nlong=81;

h=0.01;

for(s=1;s<11;s++){
c[s]=0.4;
y[s]=0.4;
}

t=0;
while(t<100){


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<11;s++){
a[0][4*s-3]=1/c[s];
a[0][4*s-2]=-1/c[s];
a[0][4*s-1]=-1/(1-y[s]);
a[0][4*s]=1/(1-y[s]);
}


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

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

for(s=1;s<42;s++){
a[s][40+s]=1;
}

for(s=1;s<41;s++){
b[s]=h;
}

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


for(s=1;s<11;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<11;s++){
System.out.println(c[s]);
System.out.println(y[s]);
}


}
}