lp10 - プログラム倉庫

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;

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

pro(){

h=0.01;
c1=0.3;
y1=0.4;
c2=0.4;
y2=0.3;

t=0;
while(t<100){

mlong=9;
nlong=17;

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

a[0][1]=1/c1;
a[0][2]=-1/c1;
a[0][3]=-1/(1-y1);
a[0][4]=1/(1-y1);
a[0][5]=1/c2;
a[0][6]=-1/c2;
a[0][7]=-1/(1-y2);
a[0][8]=1/(1-y2);

a[1][1]=1;  
a[2][2]=1;
a[3][3]=1;
a[4][4]=1;
a[5][5]=1;  
a[6][6]=1;
a[7][7]=1;
a[8][8]=1;

a[9][1]=1;  
a[9][2]=-1;
a[9][3]=-1;
a[9][4]=1;
a[9][5]=1;  
a[9][6]=-1;
a[9][7]=-1;
a[9][8]=1;


b[1]=h;
b[2]=h;
b[3]=h;
b[4]=h;
b[5]=h;
b[6]=h;
b[7]=h;
b[8]=h;
b[9]=0;

a[1][9]=1;
a[2][10]=1;
a[3][11]=1;
a[4][12]=1;
a[5][13]=1;
a[6][14]=1;
a[7][15]=1;
a[8][16]=1;
a[9][17]=1;




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

System.out.println(c1);
System.out.println(y1);

c1=c1+x[1]-x[2];
y1=y1+x[3]-x[4];
c2=c2+x[5]-x[6];
y2=y2+x[7]-x[8];

t=t+1;
}



}
}