誘引なし - プログラム倉庫

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];
double[] th=new double[11];
int s;

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

pro(){

mlong=51;
nlong=91;

h=0.01;

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

t=0;
while(t<1000){

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/(th[s]-y[s]);
a[0][4*s]=1/(th[s]-y[s]);
}


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

for(s=1;s<11;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<52;s++){
a[s][40+s]=1;
}

for(s=1;s<11;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;


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]);
}


}
}