rc = 4.0
Ec(1,2) = 2.45
re(1,2) = 2.75
alpha(1,2) = 4.10
delr = 0.1
rho0(2) = 3.0
Cmax(1,1,2) = 2.81
Cmax(2,2,1) = 2.20
Cmax(1,2,1) = 2.40
Cmax(1,2,2) = 2.40
Cmin(1,1,2) = 0.55
Cmin(2,2,1) = 0.35
Cmin(1,2,1) = 0.45
Cmin(1,2,2) = 0.45
Cmax(1,1,1) = 1.91
Cmin(1,1,1) = 0.31
Cmax(2,2,2) = 2.60
Cmin(2,2,2) = 0.75
nn2(1,1) = 1
nn2(2,2) = 1