Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- emptwob = {};
- For[j = 5, j < 7, j++,
- For[hone = 1, hone < j, hone++,
- For[htwo = 1, htwo < j, htwo++,
- If[hone == (j - 1)/2, Continue,
- If[htwo == (j - 1)/2, Continue,
- If[hone != htwo, R = 6.0000000000000005*(10)^(15);
- Q = 1.01*(10)^(32); SuM = Q + R^2; P = (10)^(50);
- B = LatticeReduce[{ {1, 0, Round[P*Log[j]]}, {0, 1,
- Round[P*Log[2*hone/(j - 1)]]}, {0, 0,
- Round[P*Log[2*htwo/(j - 1)]]} }]; D1 = B[[1]];
- D2 = B[[2]]; D3 = B[[3]]; F2 = D1;
- F3 = D2 - ((D2.F2)/(F2.F2))*F2;
- F4 = D3 - ((D3.F2)/(F2.F2))*F2 - ((D3.F3)/(F3.F3))*F3;
- epsloc = {Norm[D1]/Norm[F2], Norm[D1]/Norm[F3],
- Norm[D1]/Norm[F4] }; ceoone = Max[epsloc];
- teeta = Norm[D1]/ceoone;
- If[teeta^2 > SuM, LB = (Sqrt[teeta^2 - Q] - R)/P,
- Print[{j, hone, htwo}]]; AppendTo[emptwob, LB],
- Continue] ] ] ]]]; Print[Min[emptwob]]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement