Cody

# Problem 1462. Dots in a Sphere

Solution 1427325

Submitted on 29 Jan 2018 by cokakola
This solution is locked. To view this solution, you need to provide a solution of the same size or smaller.

### Test Suite

Test Status Code Input and Output
1   Pass
user_solution = fileread('dots_in_sphere.m'); assert(isempty(strfind(user_solution,'regexp'))); assert(isempty(strfind(user_solution,'2str'))); assert(isempty(strfind(user_solution,'str2'))); assert(isempty(strfind(user_solution,'interp'))); assert(isempty(strfind(user_solution,'printf'))); assert(isempty(strfind(user_solution,'assert')));

2   Pass
r = 0; n_correct = 1; assert(isequal(dots_in_sphere(r),n_correct));

3   Pass
r = 0.5; n_correct = 1; assert(isequal(dots_in_sphere(r),n_correct));

4   Pass
r = 1; n_correct = 7; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = 3 2 3 2 1 2 3 2 3 RR(:,:,2) = 2 1 2 1 0 1 2 1 2 RR(:,:,3) = 3 2 3 2 1 2 3 2 3

5   Pass
r = 1.5; n_correct = 19; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = 3 2 3 2 1 2 3 2 3 RR(:,:,2) = 2 1 2 1 0 1 2 1 2 RR(:,:,3) = 3 2 3 2 1 2 3 2 3

6   Pass
r = 2; n_correct = 33; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = 12 9 8 9 12 9 6 5 6 9 8 5 4 5 8 9 6 5 6 9 12 9 8 9 12 RR(:,:,2) = 9 6 5 6 9 6 3 2 3 6 5 2 1 2 5 6 3 2 3 6 9 6 5 6 9 RR(:,:,3) = 8 5 4 5 8 5 2 1 2 5 4 1 0 1 4 5 2 1 2 5 8 5 4 5 8 RR(:,:,4) = 9 6 5 6 9 6 3 2 3 6 5 2 1 2 5 6 3 2 3 6 9 6 5 6 9 RR(:,:,5) = 12 9 8 9 12 9 6 5 6 9 8 5 4 5 8 9 6 5 6 9 12 9 8 9 12

7   Pass
r = 2.5; n_correct = 81; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = 12 9 8 9 12 9 6 5 6 9 8 5 4 5 8 9 6 5 6 9 12 9 8 9 12 RR(:,:,2) = 9 6 5 6 9 6 3 2 3 6 5 2 1 2 5 6 3 2 3 6 9 6 5 6 9 RR(:,:,3) = 8 5 4 5 8 5 2 1 2 5 4 1 0 1 4 5 2 1 2 5 8 5 4 5 8 RR(:,:,4) = 9 6 5 6 9 6 3 2 3 6 5 2 1 2 5 6 3 2 3 6 9 6 5 6 9 RR(:,:,5) = 12 9 8 9 12 9 6 5 6 9 8 5 4 5 8 9 6 5 6 9 12 9 8 9 12

8   Pass
r = 3; n_correct = 123; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = 27 22 19 18 19 22 27 22 17 14 13 14 17 22 19 14 11 10 11 14 19 18 13 10 9 10 13 18 19 14 11 10 11 14 19 22 17 14 13 14 17 22 27 22 19 18 19 22 27 RR(:,:,2) = 22 17 14 13 14 17 22 17 12 9 8 9 12 17 14 9 6 5 6 9 14 13 8 5 4 5 8 13 14 9 6 5 6 9 14 17 12 9 8 9 12 17 22 17 14 13 14 17 22 RR(:,:,3) = 19 14 11 10 11 14 19 14 9 6 5 6 9 14 11 6 3 2 3 6 11 10 5 2 1 2 5 10 11 6 3 2 3 6 11 14 9 6 5 6 9 14 19 14 11 10 11 14 19 RR(:,:,4) = 18 13 10 9 10 13 18 13 8 5 4 5 8 13 10 5 2 1 2 5 10 9 4 1 0 1 4 9 10 5 2 1 2 5 10 13 8 5 4 5 8 13 18 13 10 9 10 13 18 RR(:,:,5) = 19 14 11 10 11 14 19 14 9 6 5 6 9 14 11 6 3 2 3 6 11 10 5 2 1 2 5 10 11 6 3 2 3 6 11 14 9 6 5 6 9 14 19 14 11 10 11 14 19 RR(:,:,6) = 22 17 14 13 14 17 22 17 12 9 8 9 12 17 14 9 6 5 6 9 14 13 8 5 4 5 8 13 14 9 6 5 6 9 14 17 12 9 8 9 12 17 22 17 14 13 14 17 22 RR(:,:,7) = 27 22 19 18 19 22 27 22 17 14 13 14 17 22 19 14 11 10 11 14 19 18 13 10 9 10 13 18 19 14 11 10 11 14 19 22 17 14 13 14 17 22 27 22 19 18 19 22 27

9   Pass
r = 5; n_correct = 515; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = 75 66 59 54 51 50 51 54 59 66 75 66 57 50 45 42 41 42 45 50 57 66 59 50 43 38 35 34 35 38 43 50 59 54 45 38 33 30 29 30 33 38 45 54 51 42 35 30 27 26 27 30 35 42 51 50 41 34 29 26 25 26 29 34 41 50 51 42 35 30 27 26 27 30 35 42 51 54 45 38 33 30 29 30 33 38 45 54 59 50 43 38 35 34 35 38 43 50 59 66 57 50 45 42 41 42 45 50 57 66 75 66 59 54 51 50 51 54 59 66 75 RR(:,:,2) = 66 57 50 45 42 41 42 45 50 57 66 57 48 41 36 33 32 33 36 41 48 57 50 41 34 29 26 25 26 29 34 41 50 45 36 29 24 21 20 21 24 29 36 45 42 33 26 21 18 17 18 21 26 33 42 41 32 25 20 17 16 17 20 25 32 41 42 33 26 21 18 17 18 21 26 33 42 45 36 29 24 21 20 21 24 29 36 45 50 41 34 29 26 25 26 29 34 41 50 57 48 41 36 33 32 33 36 41 48 57 66 57 50 45 42 41 42 45 50 57 66 RR(:,:,3) = 59 50 43 38 35 34 35 38 43 50 59 50 41 34 29 26 25 26 29 34 41 50 43 34 27 22 19 18 19 22 27 34 43 38 29 22 17 14 13 14 17 22 29 38 35 26 19 14 11 10 11 14 19 26 35 34 25 18 13 10 9 10 13 18 25 34 35 26 19 14 11 10 11 14 19 26 35 38 29 22 17 14 13 14 17 22 29 38 43 34 27 22 19 18 19 22 27 34 43 50 41 34 29 26 25 26 29 34 41 50 59 50 43 38 35 34 35 38 43 50 59 RR(:,:,4) = 54 45 38 33 30 29 30 33 38 45 54 45 36 29 24 21 20 21 24 29 36 45 38 29 22 17 14 13 14 17 22 29 38 33 24 17 12 9 8 9 12 17 24 33 30 21 14 9 6 5 6 9 14 21 30 29 20 13 8 5 4 5 8 13 20 29 30 21 14 9 6 5 6 9 14 21 30 33 24 17 12 9 8 9 12 17 24 33 38 29 22 17 14 13 14 17 22 29 38 45 36 29 24 21 20 21 24 29 36 45 54 45 38 33 30 29 30 33 38 45 54 RR(:,:,5) = 51 42 35 30 27 26 27 30 35 42 51 42 33 26 21 18 17 18 21 26 33 42 35 26 19 14 11 10 11 14 19 26 35 30 21 14 9 6 5 6 9 14 21 30 27 18 11 6 3 2 3 6 11 18 27 26 17 10 5 2 1 2 5 10 17 26 27 18 11 6 3 2 3 6 11 18 27 30 21 14 9 6 5 6 9 14 21 30 35 26 19 14 11 10 11 14 19 26 35 42 33 26 21 18 17 18 21 26 33 42 51 42 35 30 27 26 27 30 35 42 51 RR(:,:,6) = 50 41 34 29 26 25 26 29 34 41 50 41 32 25 20 17 16 17 20 25 32 41 34 25 18 13 10 9 10 13 18 25 34 29 20 13 8 5 4 5 8 13 20 29 26 17 10 5 2 1 2 5 10 17 26 25 16 9 4 1 0 1 4 9 16 25 26 17 10 5 2 1 2 5 10 17 26 29 20 13 8 5 4 5 8 13 20 29 34 25 18 13 10 9 10 13 18 25 34 41 32 25 20 17 16 17 20 25 32 41 50 41 34 29 26 25 26 29 34 41 50 RR(:,:,7) = 51 42 35 30 27 26 27 30 35 42 51 42 33 26 21 18 17 18 21 26 33 42 35 26 19 14 11 10 11 14 19 26 35 30 21 14 9 6 5 6 9 14 21 30 27 18 11 6 3 2 3 6 11 18 27 26 17 10 5 2 1 2 5 10 17 26 27 18 11 6 3 2 3 6 11 18 27 30 21 14 9 6 5 6 9 14 21 30 35 26 19 14 11 10 11 14 19 26 35 42 33 26 21 18 17 18 21 26 33 42 51 42 35 30 27 26 27 30 35 42 51 RR(:,:,8) = 54 45 38 33 30 29 30 33 38 45 54 45 36 29 24 21 20 21 24 29 36 45 38 29 22 17 14 13 14 17 22 29 38 33 24 17 12 9 8 9 12 17 24 33 30 21 14 9 6 5 6 9 14 21 30 29 20 13 8 5 4 5 8 13 20 29 30 21 14 9 6 5 6 9 14 21 30 33 24 17 12 9 8 9 12 17 24 33 38 29 22 17 14 13 14 17 22 29 38 45 36 29 24 21 20 21 24 29 36 45 54 45 38 33 30 29 30 33 38 45 54 RR(:,:,9) = 59 50 43 38 35 34 35 38 43 50 59 50 41 34 29 26 25 26 29 34 41 50 43 34 27 22 19 18 19 22 27 34 43 38 29 22 17 14 13 14 17 22 29 38 35 26 19 14 11 10 11 14 19 26 35 34 25 18 13 10 9 10 13 18 25 34 35 26 19 14 11 10 11 14 19 26 35 38 29 22 17 14 13 14 17 22 29 38 43 34 27 22 19 18 19 22 27 34 43 50 41 34 29 26 25 26 29 34 41 50 59 50 43 38 35 34 35 38 43 50 59 RR(:,:,10) = 66 57 50 45 42 41 42 45 50 57 66 57 48 41 36 33 32 33 36 41 48 57 50 41 34 29 26 25 26 29 34 41 50 45 36 29 24 21 20 21 24 29 36 45 42 33 26 21 18 17 18 21 26 33 42 41 32 25 20 17 16 17 20 25 32 41 42 33 26 21 18 17 18 21 26 33 42 45 36 29 24 21 20 21 24 29 36 45 50 41 34 29 26 25 26 29 34 41 50 57 48 41 36 33 32 33 36 41 48 57 66 57 50 45 42 41 42 45 50 57 66 RR(:,:,11) = 75 66 59 54 51 50 51 54 59 66 75 66 57 50 45 42 41 42 45 50 57 66 59 50 43 38 35 34 35 38 43 50 59 54 45 38 33 30 29 30 33 38 45 54 51 42 35 30 27 26 27 30 35 42 51 50 41 34 29 26 25 26 29 34 41 50 51 42 35 30 27 26 27 30 35 42 51 54 45 38 33 30 29 30 33 38 45 54 59 50 43 38 35 34 35 38 43 50 59 66 57 50 45 42 41 42 45 50 57 66 75 66 59 54 51 50 51 54 59 66 75

10   Pass
r = 7.5; n_correct = 1791; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = 147 134 123 114 107 102 99 98 99 102 107 114 123 134 147 134 121 110 101 94 89 86 85 86 89 94 101 110 121 134 123 110 99 90 83 78 75 74 75 78 83 90 99 110 123 114 101 90 81 74 69 66 65 66 69 74 81 90 101 114 107 94 83 74 67 62 59 58 59 62 67 74 83 94 107 102 89 78 69 62 57 54 53 54 57 62 69 78 89 102 99 86 75 66 59 54 51 50 51 54 59 66 75 86 99 98 85 74 65 58 53 50 49 50 53 58 65 74 85 98 99 86 75 66 59 54 51 50 51 54 59 66 75 86 99 102 89 78 69 62 57 54 53 54 57 62 69 78 89 102 107 94 83 74 67 62 59 58 59 62 67 74 83 94 107 114 101 90 81 74 69 66 65 66 69 74 81 90 101 114 123 110 99 90 83 78 75 74 75 78 83 90 99 110 123 134 121 110 101 94 89 86 85 86 89 94 101 110 121 134 147 134 123 114 107 102 99 98 99 102 107 114 123 134 147 RR(:,:,2) = 134 121 110 101 94 89 86 85 86 89 94 101 110 121 134 121 108 97 88 81 76 73 72 73 76 81 88 97 108 121 110 97 86 77 70 65 62 61 62 65 70 77 86 97 110 101 88 77 68 61 56 53 52 53 56 61 68 77 88 101 94 81 70 61 54 49 46 45 46 49 54 61 70 81 94 89 76 65 56 49 44 41 40 41 44 49 56 65 76 89 86 73 62 53 46 41 38 37 38 41 46 53 62 73 86 85 72 61 52 45 40 37 36 37 40 45 52 61 72 85 86 73 62 53 46 41 38 37 38 41 46 53 62 73 86 89 76 65 56 49 44 41 40 41 44 49 56 65 76 89 94 81 70 61 54 49 46 45 46 49 54 61 70 81 94 101 88 77 68 61 56 53 52 53 56 61 68 77 88 101 110 97 86 77 70 65 62 61 62 65 70 77 86 97 110 121 108 97 88 81 76 73 72 73 76 81 88 97 108 121 134 121 110 101 94 89 86 85 86 89 94 101 110 121 134 RR(:,:,3) = 123 110 99 90 83 78 75 74 75 78 83 90 99 110 123 110 97 86 77 70 65 62 61 62 65 70 77 86 97 110 99 86 75 66 59 54 51 50 51 54 59 66 75 86 99 90 77 66 57 50 45 42 41 42 45 50 57 66 77 90 83 70 59 50 43 38 35 34 35 38 43 50 59 70 83 78 65 54 45 38 33 30 29 30 33 38 45 54 65 78 75 62 51 42 35 30 27 26 27 30 35 42 51 62 75 74 61 50 41 34 29 26 25 26 29 34 41 50 61 74 75 62 51 42 35 30 27 26 27 30 35 42 51 62 75 78 65 54 45 38 33 30 29 30 33 38 45 54 65 78 83 70 59 50 43 38 35 34 35 38 43 50 59 70 83 90 77 66 57 50 45 42 41 42 45 50 57 66 77 90 99 86 75 66 59 54 51 50 51 54 59 66 75 86 99 110 97 86 77 70 65 62 61 62 65 70 77 86 97 110 123 110 99 90 83 78 75 74 75 78 83 90 99 110 123 RR(:,:,4) = 114 101 90 81 74 69 66 65 66 69 74 81 90 101 114 101 88 77 68 61 56 53 52 53 56 61 68 77 88 101 90 77 66 57 50 45 42 41 42 45 50 57 66 77 90 81 68 57 48 41 36 33 32 33 36 41 48 57 68 81 74 61 50 41 34 29 26 25 26 29 34 41 50 61 74 69 56 45 36 29 24 21 20 21 24 29 36 45 56 69 66 53 42 33 26 21 18 17 18 21 26 33 42 53 66 65 52 41 32 25 20 17 16 17 20 25 32 41 52 65 66 53 42 33 26 21 18 17 18 21 26 33 42 53 66 69 56 45 36 29 24 21 20 21 24 29 36 45 56 69 74 61 50 41 34 29 26 25 26 29 34 41 50 61 74 81 68 57 48 41 36 33 32 33 36 41 48 57 68 81 90 77 66 57 50 45 42 41 42 45 50 57 66 77 90 101 88 77 68 61 56 53 52 53 56 61 68 77 88 101 114 101 90 81 74 69 66 65 66 69 74 81 90 101 114 RR(:,:,5) = 107 94 83 74 67 62 59 58 59 62 67 74 83 94 107 94 81 70 61 54 49 46 45 46 49 54 61 70 81 94 83 70 59 50 43 38 35 34 35 38 43 50 59 70 83 74 61 50 41 34 29 26 25 26 29 34 41 50 61 74 67 54 43 34 27 22 19 18 19 22 27 34 43 54 67 62 49 38 29 22 17 14 13 14 17 22 29 38 49 62 59 46 35 26 19 14 11 10 11 14 19 26 35 46 59 58 45 34 25 18 13 10 9 10 13 18 25 34 45 58 59 46 35 26 19 14 11 10 11 14 19 26 35 46 59 62 49 38 29 22 17 14 13 14 17 22 29 38 49 62 67 54 43 34 27 22 19 18 19 22 27 34 43 54 67 74 61 50 41 34 29 26 25 26 29 34 41 50 61 74 83 70 59 50 43 38 35 34 35 38 43 50 59 70 83 94 81 70 61 54 49 46 45 46 49 54 61 70 81 94 107 94 83 74 67 62 59 58 59 62 67 74 83 94 107 RR(:,:,6) = 102 89 78 69 62 57 54 53 54 57 62 69 78 89 102 89 76 65 56 49 44 41 40 41 44 49 56 65 76 89 78 65 54 45 38 33 30 29 30 33 38 45 54 65 78 69 56 45 36 29 24 21 20 21 24 29 36 45 56 69 62 49 38 29 22 17 14 13 14 17 22 29 38 49 62 57 44 33 24 17 12 9 8 9 12 17 24 33 44 57 54 41 30 21 14 9 6 5 6 9 14 21 30 41 54 53 40 29 20 13 8 5 4 5 8 13 20 29 40 53 54 41 30 21 14 9 6 5 6 9 14 21 30 41 54 57 44 33 24 17 12 9 8 9 12 17 24 33 44 57 62 49 38 29 22 17 14 13 14 17 22 29 38 49 62 69 56 45 36 29 24 21 20 21 24 29 36 45 56 69 78 65 54 45 38 33 30 29 30 33 38 45 54 65 78 89 76 65 56 49 44 41 40 41 44 49 56 65 76 89 102 89 78 69 62 57 54 53 54 57 62 69 78 89 102 RR(:,:,7) = 99 86 75 66 59 54 51 50 51 54 59 66 75 86 99 86 73 62 53 46 41 38 37 38 41 46 53 62 73 86 75 62 51 42 35 30 27 26 27 30 35 42 51 62 75 66 53 42 33 26 21 18 17 18 21 26 33 42 53 66 59 46 35 26 19 14 11 10 11 14 19 26 35 46 59 54 41 30 21 14 9 6 5 6 9 14 21 30 41 54 51 38 27 18 11 6 3 2 3 6 11 18 27 38 51 50 37 26 17 10 5 2 1 2 5 10 17 26 37 50 51 38 27 18 11 6 3 2 3 6 11 18 27 38 51 54 41 30 21 14 9 6 5 6 9 14 21 30 41 54 59 46 35 26 19 14 11 10 11 14 19 26 35 46 59 66 53 42 33 26 21 18 17 18 21 26 33 42 53 66 75 62 51 42 35 30 27 26 27 30 35 42 51 62 75 86 73 62 53 46 41 38 37 38 41 46 53 62 73 86 99 86 75 66 59 54 51 50 51 54 59 66 75 86 99 RR(:,:,8) = 98 85 74 65 58 53 50 49 50 53 58 65 74 85 98 85 72 61 52 45 40 37 36 37 40 45 52 61 72 85 74 61 50 41 34 29 26 25 26 29 34 41 50 61 74 65 52 41 32 25 20 17 16 17...

11   Pass
r = 10; n_correct = 4169; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = 300 281 264 249 236 225 216 209 204 201 200 201 204 209 216 225 236 249 264 281 300 281 262 245 230 217 206 197 190 185 182 181 182 185 190 197 206 217 230 245 262 281 264 245 228 213 200 189 180 173 168 165 164 165 168 173 180 189 200 213 228 245 264 249 230 213 198 185 174 165 158 153 150 149 150 153 158 165 174 185 198 213 230 249 236 217 200 185 172 161 152 145 140 137 136 137 140 145 152 161 172 185 200 217 236 225 206 189 174 161 150 141 134 129 126 125 126 129 134 141 150 161 174 189 206 225 216 197 180 165 152 141 132 125 120 117 116 117 120 125 132 141 152 165 180 197 216 209 190 173 158 145 134 125 118 113 110 109 110 113 118 125 134 145 158 173 190 209 204 185 168 153 140 129 120 113 108 105 104 105 108 113 120 129 140 153 168 185 204 201 182 165 150 137 126 117 110 105 102 101 102 105 110 117 126 137 150 165 182 201 200 181 164 149 136 125 116 109 104 101 100 101 104 109 116 125 136 149 164 181 200 201 182 165 150 137 126 117 110 105 102 101 102 105 110 117 126 137 150 165 182 201 204 185 168 153 140 129 120 113 108 105 104 105 108 113 120 129 140 153 168 185 204 209 190 173 158 145 134 125 118 113 110 109 110 113 118 125 134 145 158 173 190 209 216 197 180 165 152 141 132 125 120 117 116 117 120 125 132 141 152 165 180 197 216 225 206 189 174 161 150 141 134 129 126 125 126 129 134 141 150 161 174 189 206 225 236 217 200 185 172 161 152 145 140 137 136 137 140 145 152 161 172 185 200 217 236 249 230 213 198 185 174 165 158 153 150 149 150 153 158 165 174 185 198 213 230 249 264 245 228 213 200 189 180 173 168 165 164 165 168 173 180 189 200 213 228 245 264 281 262 245 230 217 206 197 190 185 182 181 182 185 190 197 206 217 230 245 262 281 300 281 264 249 236 225 216 209 204 201 200 201 204 209 216 225 236 249 264 281 300 RR(:,:,2) = 281 262 245 230 217 206 197 190 185 182 181 182 185 190 197 206 217 230 245 262 281 262 243 226 211 198 187 178 171 166 163 162 163 166 171 178 187 198 211 226 243 262 245 226 209 194 181 170 161 154 149 146 145 146 149 154 161 170 181 194 209 226 245 230 211 194 179 166 155 146 139 134 131 130 131 134 139 146 155 166 179 194 211 230 217 198 181 166 153 142 133 126 121 118 117 118 121 126 133 142 153 166 181 198 217 206 187 170 155 142 131 122 115 110 107 106 107 110 115 122 131 142 155 170 187 206 197 178 161 146 133 122 113 106 101 98 97 98 101 106 113 122 133 146 161 178 197 190 171 154 139 126 115 106 99 94 91 90 91 94 99 106 115 126 139 154 171 190 185 166 149 134 121 110 101 94 89 86 85 86 89 94 101 110 121 134 149 166 185 182 163 146 131 118 107 98 91 86 83 82 83 86 91 98 107 118 131 146 163 182 181 162 145 130 117 106 97 90 85 82 81 82 85 90 97 106 117 130 145 162 181 182 163 146 131 118 107 98 91 86 83 82 83 86 91 98 107 118 131 146 163 182 185 166 149 134 121 110 101 94 89 86 85 86 89 94 101 110 121 134 149 166 185 190 171 154 139 126 115 106 99 94 91 90 91 94 99 106 115 126 139 154 171 190 197 178 161 146 133 122 113 106 101 98 97 98 101 106 113 122 133 146 161 178 197 206 187 170 155 142 131 122 115 110 107 106 107 110 115 122 131 142 155 170 187 206 217 198 181 166 153 142 133 126 121 118 117 118 121 126 133 142 153 166 181 198 217 230 211 194 179 166 155 146 139 134 131 130 131 134 139 146 155 166 179 194 211 230 245 226 209 194 181 170 161 154 149 146 145 146 149 154 161 170 181 194 209 226 245 262 243 226 211 198 187 178 171 166 163 162 163 166 171 178 187 198 211 226 243 262 281 262 245 230 217 206 197 190 185 182 181 182 185 190 197 206 217 230 245 262 281 RR(:,:,3) = 264 245 228 213 200 189 180 173 168 165 164 165 168 173 180 189 200 213 228 245 264 245 226 209 194 181 170 161 154 149 146 145 146 149 154 161 170 181 194 209 226 245 228 209 192 177 164 153 144 137 132 129 128 129 132 137 144 153 164 177 192 209 228 213 194 177 162 149 138 129 122 117 114 113 114 117 122 129 138 149 162 177 194 213 200 181 164 149 136 125 116 109 104 101 100 101 104 109 116 125 136 149 164 181 200 189 170 153 138 125 114 105 98 93 90 89 90 93 98 105 114 125 138 153 170 189 180 161 144 129 116 105 96 89 84 81 80 81 84 89 96 105 116 129 144 161 180 173 154 137 122 109 98 89 82 77 74 73 74 77 82 89 98 109 122 137 154 173 168 149 132 117 104 93 84 77 72 69 68 69 72 77 84 93 104 117 132 149 168 165 146 129 114 101 90 81 74 69 66 65 66 69 74 81 90 101 114 129 146 165 164 145 128 113 100 89 80 73 68 65 64 65 68 73 80 89 100 113 128 145 164 165 146 129 114 101 90 81 74 69 66 65 66 69 74 81 90 101 114 129 146 165 168 149 132 117 104 93 84 77 72 69 68 69 72 77 84 93 104 117 132 149 168 173 154 137 122 109 98 89 82 77 74 73 74 77 82 89 98 109 122 137 154 173 180 161 144 129 116 105 96 89 84 81 80 81 84 89 96 105 116 129 144 161 180 189 170 153 138 125 114 105 98 93 90 89 90 93 98 105 114 125 138 153 170 189 200 181 164 149 136 125 116 109 104 101 100 101 104 109 116 125 136 149 164 181 200 213 194 177 162 149 138 129 122 117 114 113 114 117 122 129 138 149 162 177 194 213 228 209 192 177 164 153 144 137 132 129 128 129 132 137 144 153 164 177 192 209 228 245 226 209 194 181 170 161 154 149 146 145 146 149 154 161 170 181 194 209 226 245 264 245 228 213 200 189 180 173 168 165 164 165 168 173 180 189 200 213 228 245 264 RR(:,:,4) = 249 230 213 198 185 174 165 158 153 150 149 150 153 158 165 174 185 198 213 230 249 230 211 194 179 166 155 146 139 134 131 130 131 134 139 146 155 166 179 194 211 230 213 194 177 162 149 138 129 122 117 114 113 114 117 122 129 138 149 162 177 194 213 198 179 162 147 134 123 114 107 102 99 98 99 102 107 114 123 134 147 162 179 198 185 166 149 134 121 110 101 94 89 86 85 86 89 94 101 110 121 134 149 166 185 174 155 138 123 110 99 90 83 78 75 74 75 78 83 90 99 110 123 138 155 174 165 146 129 114 101 90 81 74 69 66 65 66 69 74 81 90 101 114 129 146 165 158 139 122 107 94 83 74 67 62 59 58 59 62 67 74 83 94 107 122 139 158 153 134 117 102 89 78 69 62 57 54 53 54 57 62 69 78 89 102 117 134 153 150 131 114 99 86 75 66 59 54 51 50 51 54 59 66 75 86 99 114 131 150 149 130 113 98 85 74 65 58 53 50 49 50 53 58 65 74 85 98 113 130 149 150 131 114 99 86 75 66 59 54 51 50 51 54 59 66 75 86 99 114 131 150 153 134 117 102 89 78 69 62 57 54 53 54 57 62 69 78 89 102 117 134 153 158 139 122 107 94 83 74 67 62 59 58 59 62 67 74 83 94 107 122 139 158 165 146 129 114 101 90 81 74 69 66 65 66 69 74 81 90 101 114 129 146 165 174 155 138 123 110 99...

12   Pass
r = 15; n_correct = 14147; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = Columns 1 through 30 675 646 619 594 571 550 531 514 499 486 475 466 459 454 451 450 451 454 459 466 475 486 499 514 531 550 571 594 619 646 646 617 590 565 542 521 502 485 470 457 446 437 430 425 422 421 422 425 430 437 446 457 470 485 502 521 542 565 590 617 619 590 563 538 515 494 475 458 443 430 419 410 403 398 395 394 395 398 403 410 419 430 443 458 475 494 515 538 563 590 594 565 538 513 490 469 450 433 418 405 394 385 378 373 370 369 370 373 378 385 394 405 418 433 450 469 490 513 538 565 571 542 515 490 467 446 427 410 395 382 371 362 355 350 347 346 347 350 355 362 371 382 395 410 427 446 467 490 515 542 550 521 494 469 446 425 406 389 374 361 350 341 334 329 326 325 326 329 334 341 350 361 374 389 406 425 446 469 494 521 531 502 475 450 427 406 387 370 355 342 331 322 315 310 307 306 307 310 315 322 331 342 355 370 387 406 427 450 475 502 514 485 458 433 410 389 370 353 338 325 314 305 298 293 290 289 290 293 298 305 314 325 338 353 370 389 410 433 458 485 499 470 443 418 395 374 355 338 323 310 299 290 283 278 275 274 275 278 283 290 299 310 323 338 355 374 395 418 443 470 486 457 430 405 382 361 342 325 310 297 286 277 270 265 262 261 262 265 270 277 286 297 310 325 342 361 382 405 430 457 475 446 419 394 371 350 331 314 299 286 275 266 259 254 251 250 251 254 259 266 275 286 299 314 331 350 371 394 419 446 466 437 410 385 362 341 322 305 290 277 266 257 250 245 242 241 242 245 250 257 266 277 290 305 322 341 362 385 410 437 459 430 403 378 355 334 315 298 283 270 259 250 243 238 235 234 235 238 243 250 259 270 283 298 315 334 355 378 403 430 454 425 398 373 350 329 310 293 278 265 254 245 238 233 230 229 230 233 238 245 254 265 278 293 310 329 350 373 398 425 451 422 395 370 347 326 307 290 275 262 251 242 235 230 227 226 227 230 235 242 251 262 275 290 307 326 347 370 395 422 450 421 394 369 346 325 306 289 274 261 250 241 234 229 226 225 226 229 234 241 250 261 274 289 306 325 346 369 394 421 451 422 395 370 347 326 307 290 275 262 251 242 235 230 227 226 227 230 235 242 251 262 275 290 307 326 347 370 395 422 454 425 398 373 350 329 310 293 278 265 254 245 238 233 230 229 230 233 238 245 254 265 278 293 310 329 350 373 398 425 459 430 403 378 355 334 315 298 283 270 259 250 243 238 235 234 235 238 243 250 259 270 283 298 315 334 355 378 403 430 466 437 410 385 362 341 322 305 290 277 266 257 250 245 242 241 242 245 250 257 266 277 290 305 322 341 362 385 410 437 475 446 419 394 371 350 331 314 299 286 275 266 259 254 251 250 251 254 259 266 275 286 299 314 331 350 371 394 419 446 486 457 430 405 382 361 342 325 310 297 286 277 270 265 262 261 262 265 270 277 286 297 310 325 342 361 382 405 430 457 499 470 443 418 395 374 355 338 323 310 299 290 283 278 275 274 275 278 283 290 299 310 323 338 355 374 395 418 443 470 514 485 458 433 410 389 370 353 338 325 314 305 298 293 290 289 290 293 298 305 314 325 338 353 370 389 410 433 458 485 531 502 475 450 427 406 387 370 355 342 331 322 315 310 307 306 307 310 315 322 331 342 355 370 387 406 427 450 475 502 550 521 494 469 446 425 406 389 374 361 350 341 334 329 326 325 326 329 334 341 350 361 374 389 406 425 446 469 494 521 571 542 515 490 467 446 427 410 395 382 371 362 355 350 347 346 347 350 355 362 371 382 395 410 427 446 467 490 515 542 594 565 538 513 490 469 450 433 418 405 394 385 378 373 370 369 370 373 378 385 394 405 418 433 450 469 490 513 538 565 619 590 563 538 515 494 475 458 443 430 419 410 403 398 395 394 395 398 403 410 419 430 443 458 475 494 515 538 563 590 646 617 590 565 542 521 502 485 470 457 446 437 430 425 422 421 422 425 430 437 446 457 470 485 502 521 542 565 590 617 675 646 619 594 571 550 531 514 499 486 475 466 459 454 451 450 451 454 459 466 475 486 499 514 531 550 571 594 619 646 Column 31 675 646 619 594 571 550 531 514 499 486 475 466 459 454 451 450 451 454 459 466 475 486 499 514 531 550 571 594 619 646 675 RR(:,:,2) = Columns 1 through 30 646 617 590 565 542 521 502 485 470 457 446 437 430 425 422 421 422 425 430 437 446 457 470 485 502 521 542 565 590 617 617 588 561 536 513 492 473 456 441 428 417 408 401 396 393 392 393 396 401 408 417 428 441 456 473 492 513 536 561 588 590 561 534 509 486 465 446 429 414 401 390 381 374 369 366 365 366 369 374 381 390 401 414 429 446 465 486 509 534 561 565 536 509 484 461 440 421 404 389 376 365 356 349 344 341 340 341 344 349 356 365 376 389 404 421 440 461 484 509 536 542 513 486 461 438 417 398 381 366 353 342 333 326 321 318 317 318 321 326 333 342 353 366 381 398 417 438 461 486 513 521 492 465 440 417 396 377 360 345 332 321 312 305 300 297 296 297 300 305 312 321 332 345 360 377 396 417 440 465 492 502 473 446 421 398 377 358 341 326 313 302 293 286 281 278 277 278 281 286 293 302 313 326 341 358 377 398 421 446 473 485 456 429 404 381 360 341 324 309 296 285 276 269 264 261 260 261 264 269 276 285 296 309 324 341 360 381 404 429 456 470 441 414 389 366 345 326 309 294 281 270 261 254 249 246 245 246 249 254 261 270 281 294 309 326 345 366 389 414 441 457 428 401 376 353 332 313 296 281 268 257 248 241 236 233 232 233 236 241 248 257 268 281 296 313 332 353 376 401 428 446 417 390 365 342 321 302 285 270 257 246 237 230 225 222 221 222 225 230 237 246 257 270 285 302 321 342 365 390 417 437 408 381 356 333 312 293 276 261 248 237 228 221 216 213 212 213 216 221 228 237 248 261 276 293 312 333 356 381 408 430 401 374 349 326 305 286 269 254 241 230 221 214 209 206 205 206 209 214 221 230 241 254 269 286 305 326 349 374 401 425 396 369 344 321 300 281 264 249 236 225 216 209 204 201 200 201 204 209 216 225 236 249 264 281 300 321 344 369 396 422 393 366 341 318 297 278 261 246 233 222 213 206 201 198 197 198 201 206 213 222 233 246 261 278 297 318 341 366 393 421 392 365 340 317 296 277 260 245 232 221 212 205 200 197 196 197 200 205 212 221 232 245 260 277 296 317 340 365 392 422 393 366 341 318 297 278 261 246 233 222 213 206 201 198 197 198 201 206 213 222 233 246 261 278 297 318 341 366 393 425 396 369 344 321 300 281 264 249 236 225 216 209 204 201 200 201 204 209 216 225 236 249 264 281 300 321 344 369 396 430 401 374 349 326 305 286 269 254 241 230 221 214 209 206 205 206 209 214 221 230 241 254 269 286 305 326 349 374 401 437 408 381 356 333 312 293 276 261 248 237 228 221 216 213 212 213 216 221 228 237 248 261 276 293 312 333 356 381 408 446 417 390 365 342 321 302 285 270 257 246 237 230 225 222 221 222 225 230 237 246 257 270 285 302 321 342 365 390 417 457 428 401 376 353 332 313 296 281 268 257 248 241 236 233 232 233 236 241 248 257 268 281 296 313 332 353 376 401 428 470 441 414 389 366 345 326 309 294 281 270 261 254 249 246 245 2...

13   Pass
r = 20; n_correct = 33401; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = Columns 1 through 15 1200 1161 1124 1089 1056 1025 996 969 944 921 900 881 864 849 836 1161 1122 1085 1050 1017 986 957 930 905 882 861 842 825 810 797 1124 1085 1048 1013 980 949 920 893 868 845 824 805 788 773 760 1089 1050 1013 978 945 914 885 858 833 810 789 770 753 738 725 1056 1017 980 945 912 881 852 825 800 777 756 737 720 705 692 1025 986 949 914 881 850 821 794 769 746 725 706 689 674 661 996 957 920 885 852 821 792 765 740 717 696 677 660 645 632 969 930 893 858 825 794 765 738 713 690 669 650 633 618 605 944 905 868 833 800 769 740 713 688 665 644 625 608 593 580 921 882 845 810 777 746 717 690 665 642 621 602 585 570 557 900 861 824 789 756 725 696 669 644 621 600 581 564 549 536 881 842 805 770 737 706 677 650 625 602 581 562 545 530 517 864 825 788 753 720 689 660 633 608 585 564 545 528 513 500 849 810 773 738 705 674 645 618 593 570 549 530 513 498 485 836 797 760 725 692 661 632 605 580 557 536 517 500 485 472 825 786 749 714 681 650 621 594 569 546 525 506 489 474 461 816 777 740 705 672 641 612 585 560 537 516 497 480 465 452 809 770 733 698 665 634 605 578 553 530 509 490 473 458 445 804 765 728 693 660 629 600 573 548 525 504 485 468 453 440 801 762 725 690 657 626 597 570 545 522 501 482 465 450 437 800 761 724 689 656 625 596 569 544 521 500 481 464 449 436 801 762 725 690 657 626 597 570 545 522 501 482 465 450 437 804 765 728 693 660 629 600 573 548 525 504 485 468 453 440 809 770 733 698 665 634 605 578 553 530 509 490 473 458 445 816 777 740 705 672 641 612 585 560 537 516 497 480 465 452 825 786 749 714 681 650 621 594 569 546 525 506 489 474 461 836 797 760 725 692 661 632 605 580 557 536 517 500 485 472 849 810 773 738 705 674 645 618 593 570 549 530 513 498 485 864 825 788 753 720 689 660 633 608 585 564 545 528 513 500 881 842 805 770 737 706 677 650 625 602 581 562 545 530 517 900 861 824 789 756 725 696 669 644 621 600 581 564 549 536 921 882 845 810 777 746 717 690 665 642 621 602 585 570 557 944 905 868 833 800 769 740 713 688 665 644 625 608 593 580 969 930 893 858 825 794 765 738 713 690 669 650 633 618 605 996 957 920 885 852 821 792 765 740 717 696 677 660 645 632 1025 986 949 914 881 850 821 794 769 746 725 706 689 674 661 1056 1017 980 945 912 881 852 825 800 777 756 737 720 705 692 1089 1050 1013 978 945 914 885 858 833 810 789 770 753 738 725 1124 1085 1048 1013 980 949 920 893 868 845 824 805 788 773 760 1161 1122 1085 1050 1017 986 957 930 905 882 861 842 825 810 797 1200 1161 1124 1089 1056 1025 996 969 944 921 900 881 864 849 836 Columns 16 through 30 825 816 809 804 801 800 801 804 809 816 825 836 849 864 881 786 777 770 765 762 761 762 765 770 777 786 797 810 825 842 749 740 733 728 725 724 725 728 733 740 749 760 773 788 805 714 705 698 693 690 689 690 693 698 705 714 725 738 753 770 681 672 665 660 657 656 657 660 665 672 681 692 705 720 737 650 641 634 629 626 625 626 629 634 641 650 661 674 689 706 621 612 605 600 597 596 597 600 605 612 621 632 645 660 677 594 585 578 573 570 569 570 573 578 585 594 605 618 633 650 569 560 553 548 545 544 545 548 553 560 569 580 593 608 625 546 537 530 525 522 521 522 525 530 537 546 557 570 585 602 525 516 509 504 501 500 501 504 509 516 525 536 549 564 581 506 497 490 485 482 481 482 485 490 497 506 517 530 545 562 489 480 473 468 465 464 465 468 473 480 489 500 513 528 545 474 465 458 453 450 449 450 453 458 465 474 485 498 ...

14   Pass
r = 25; n_correct = 65267; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = Columns 1 through 15 1875 1826 1779 1734 1691 1650 1611 1574 1539 1506 1475 1446 1419 1394 1371 1826 1777 1730 1685 1642 1601 1562 1525 1490 1457 1426 1397 1370 1345 1322 1779 1730 1683 1638 1595 1554 1515 1478 1443 1410 1379 1350 1323 1298 1275 1734 1685 1638 1593 1550 1509 1470 1433 1398 1365 1334 1305 1278 1253 1230 1691 1642 1595 1550 1507 1466 1427 1390 1355 1322 1291 1262 1235 1210 1187 1650 1601 1554 1509 1466 1425 1386 1349 1314 1281 1250 1221 1194 1169 1146 1611 1562 1515 1470 1427 1386 1347 1310 1275 1242 1211 1182 1155 1130 1107 1574 1525 1478 1433 1390 1349 1310 1273 1238 1205 1174 1145 1118 1093 1070 1539 1490 1443 1398 1355 1314 1275 1238 1203 1170 1139 1110 1083 1058 1035 1506 1457 1410 1365 1322 1281 1242 1205 1170 1137 1106 1077 1050 1025 1002 1475 1426 1379 1334 1291 1250 1211 1174 1139 1106 1075 1046 1019 994 971 1446 1397 1350 1305 1262 1221 1182 1145 1110 1077 1046 1017 990 965 942 1419 1370 1323 1278 1235 1194 1155 1118 1083 1050 1019 990 963 938 915 1394 1345 1298 1253 1210 1169 1130 1093 1058 1025 994 965 938 913 890 1371 1322 1275 1230 1187 1146 1107 1070 1035 1002 971 942 915 890 867 1350 1301 1254 1209 1166 1125 1086 1049 1014 981 950 921 894 869 846 1331 1282 1235 1190 1147 1106 1067 1030 995 962 931 902 875 850 827 1314 1265 1218 1173 1130 1089 1050 1013 978 945 914 885 858 833 810 1299 1250 1203 1158 1115 1074 1035 998 963 930 899 870 843 818 795 1286 1237 1190 1145 1102 1061 1022 985 950 917 886 857 830 805 782 1275 1226 1179 1134 1091 1050 1011 974 939 906 875 846 819 794 771 1266 1217 1170 1125 1082 1041 1002 965 930 897 866 837 810 785 762 1259 1210 1163 1118 1075 1034 995 958 923 890 859 830 803 778 755 1254 1205 1158 1113 1070 1029 990 953 918 885 854 825 798 773 750 1251 1202 1155 1110 1067 1026 987 950 915 882 851 822 795 770 747 1250 1201 1154 1109 1066 1025 986 949 914 881 850 821 794 769 746 1251 1202 1155 1110 1067 1026 987 950 915 882 851 822 795 770 747 1254 1205 1158 1113 1070 1029 990 953 918 885 854 825 798 773 750 1259 1210 1163 1118 1075 1034 995 958 923 890 859 830 803 778 755 1266 1217 1170 1125 1082 1041 1002 965 930 897 866 837 810 785 762 1275 1226 1179 1134 1091 1050 1011 974 939 906 875 846 819 794 771 1286 1237 1190 1145 1102 1061 1022 985 950 917 886 857 830 805 782 1299 1250 1203 1158 1115 1074 1035 998 963 930 899 870 843 818 795 1314 1265 1218 1173 1130 1089 1050 1013 978 945 914 885 858 833 810 1331 1282 1235 1190 1147 1106 1067 1030 995 962 931 902 875 850 827 1350 1301 1254 1209 1166 1125 1086 1049 1014 981 950 921 894 869 846 1371 1322 1275 1230 1187 1146 1107 1070 1035 1002 971 942 915 890 867 1394 1345 1298 1253 1210 1169 1130 1093 1058 1025 994 965 938 913 890 1419 1370 1323 1278 1235 1194 1155 1118 1083 1050 1019 990 963 938 915 1446 1397 1350 1305 1262 1221 1182 1145 1110 1077 1046 1017 990 965 942 1475 1426 1379 1334 1291 1250 1211 1174 1139 1106 1075 1046 1019 994 971 1506 1457 1410 1365 1322 1281 1242 1205 1170 1137 1106 1077 1050 1025 1002 1539 1490 1443 1398 1355 1314 1275 1238 1203 1170 1139 1110 1083 1058 1035 1574 1525 1478 1433 1390 1349 1310 1273 1238 1205 1174 1145 1118 1093 1070 1611 1562 1515 1470 1427 1386 1347 1310 1275 1242 1211 1182 1155 1130 1107 1650 1601 1554 1509 1466 1425 1386 1349 1314 1281 1250 1221 1194 1169 1146 1691 1642 1595 1550 1507 1466 1427 1390 1355 1322 1291 1262 1235 1210 1187 1734 1685 1638 1593 1550 1509 1470 1433 1398 1365 1334 1305 1278 1253 1230 1779 1730 1683 1638 1595 1554 1515 1478 1443 1410 1379 1350 1323 1298 1275 1826 1777 1730 1685 1642 1601 1562 1525 1490 1457 1426 1397 1370 1345 1322 1875 1826 1779 1734 1691 1650 1611 1574 1539 1506 1475 1446 1419 1394 1371 Columns 16 through 30 1350 1331 1314 1299 1286 1275 1266 1259 1254 1251 1250 1251 1254 1259 1266 1301 1282 1265 1250 1237 1226 1217 1210 1205 1202 1201 1202 1205 1210 1217 1254 1235 1218 1203 1190 1179 1170 1163 1158 1155 1154 1155 1158 1163 1170 1209 1190 1173 1158 1145 1134 1125 1118 1113 1110 1109 1110 1113 ...

15   Pass
r = 50; n_correct = 523305; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = Columns 1 through 15 7500 7401 7304 7209 7116 7025 6936 6849 6764 6681 6600 6521 6444 6369 6296 7401 7302 7205 7110 7017 6926 6837 6750 6665 6582 6501 6422 6345 6270 6197 7304 7205 7108 7013 6920 6829 6740 6653 6568 6485 6404 6325 6248 6173 6100 7209 7110 7013 6918 6825 6734 6645 6558 6473 6390 6309 6230 6153 6078 6005 7116 7017 6920 6825 6732 6641 6552 6465 6380 6297 6216 6137 6060 5985 5912 7025 6926 6829 6734 6641 6550 6461 6374 6289 6206 6125 6046 5969 5894 5821 6936 6837 6740 6645 6552 6461 6372 6285 6200 6117 6036 5957 5880 5805 5732 6849 6750 6653 6558 6465 6374 6285 6198 6113 6030 5949 5870 5793 5718 5645 6764 6665 6568 6473 6380 6289 6200 6113 6028 5945 5864 5785 5708 5633 5560 6681 6582 6485 6390 6297 6206 6117 6030 5945 5862 5781 5702 5625 5550 5477 6600 6501 6404 6309 6216 6125 6036 5949 5864 5781 5700 5621 5544 5469 5396 6521 6422 6325 6230 6137 6046 5957 5870 5785 5702 5621 5542 5465 5390 5317 6444 6345 6248 6153 6060 5969 5880 5793 5708 5625 5544 5465 5388 5313 5240 6369 6270 6173 6078 5985 5894 5805 5718 5633 5550 5469 5390 5313 5238 5165 6296 6197 6100 6005 5912 5821 5732 5645 5560 5477 5396 5317 5240 5165 5092 6225 6126 6029 5934 5841 5750 5661 5574 5489 5406 5325 5246 5169 5094 5021 6156 6057 5960 5865 5772 5681 5592 5505 5420 5337 5256 5177 5100 5025 4952 6089 5990 5893 5798 5705 5614 5525 5438 5353 5270 5189 5110 5033 4958 4885 6024 5925 5828 5733 5640 5549 5460 5373 5288 5205 5124 5045 4968 4893 4820 5961 5862 5765 5670 5577 5486 5397 5310 5225 5142 5061 4982 4905 4830 4757 5900 5801 5704 5609 5516 5425 5336 5249 5164 5081 5000 4921 4844 4769 4696 5841 5742 5645 5550 5457 5366 5277 5190 5105 5022 4941 4862 4785 4710 4637 5784 5685 5588 5493 5400 5309 5220 5133 5048 4965 4884 4805 4728 4653 4580 5729 5630 5533 5438 5345 5254 5165 5078 4993 4910 4829 4750 4673 4598 4525 5676 5577 5480 5385 5292 5201 5112 5025 4940 4857 4776 4697 4620 4545 4472 5625 5526 5429 5334 5241 5150 5061 4974 4889 4806 4725 4646 4569 4494 4421 5576 5477 5380 5285 5192 5101 5012 4925 4840 4757 4676 4597 4520 4445 4372 5529 5430 5333 5238 5145 5054 4965 4878 4793 4710 4629 4550 4473 4398 4325 5484 5385 5288 5193 5100 5009 4920 4833 4748 4665 4584 4505 4428 4353 4280 5441 5342 5245 5150 5057 4966 4877 4790 4705 4622 4541 4462 4385 4310 4237 5400 5301 5204 5109 5016 4925 4836 4749 4664 4581 4500 4421 4344 4269 4196 5361 5262 5165 5070 4977 4886 4797 4710 4625 4542 4461 4382 4305 4230 4157 5324 5225 5128 5033 4940 4849 4760 4673 4588 4505 4424 4345 4268 4193 4120 5289 5190 5093 4998 4905 4814 4725 4638 4553 4470 4389 4310 4233 4158 4085 5256 5157 5060 4965 4872 4781 4692 4605 4520 4437 4356 4277 4200 4125 4052 5225 5126 5029 4934 4841 4750 4661 4574 4489 4406 4325 4246 4169 4094 4021 5196 5097 5000 4905 4812 4721 4632 4545 4460 4377 4296 4217 4140 4065 3992 5169 5070 4973 4878 4785 4694 4605 4518 4433 4350 4269 4190 4113 4038 3965 5144 5045 4948 4853 4760 4669 4580 4493 4408 4325 4244 4165 4088 4013 3940 5121 5022 4925 4830 4737 4646 4557 4470 4385 4302 4221 4142 4065 3990 3917 5100 5001 4904 4809 4716 4625 4536 4449 4364 4281 4200 4121 4044 3969 3896 5081 4982 4885 4790 4697 4606 4517 4430 4345 4262 4181 4102 4025 3950 3877 5064 4965 4868 4773 4680 4589 4500 4413 4328 4245 4164 4085 4008 3933 3860 5049 4950 4853 4758 4665 4574 4485 4398 4313 4230 4149 4070 3993 3918 3845 5036 4937 4840 4745 4652 4561 4472 4385 4300 4217 4136 4057 3980 3905 3832 5025 4926 4829 4734 4641 4550 4461 4374 4289 4206 4125 4046 3969 3894 3821 5016 4917 4820 4725 4632 4541 4452 4365 4280 4197 4116 4037 3960 3885 3812 5009 4910 4813 4718 4625 4534 4445 4358 4273 4190 4109 4030 3953 3878 3805 5004 4905 4808 4713 4620 4529 4440 4353 4268 4185 4104 4025 3948 3873 3800 5001 4902 4805 4710 4617 4526 4437 4350 4265 4182 4101 4022 3945 3870 3797 5000 4901 4804 4709 4616 4525 4436 4349 4264 4181 4100 4021 3944 3869 3796 5001 4902 4805 4710 4617 4526 4437 4350 4265 4182 4101 4022 3945 3870 3797 5004 4905 4808 4713 4620 4529 4440 4353 4268 4185 4104 4025 3948 3873 3800 5009 4910 4813 4718 4625 4534 4445 4358 4273 4190 4109 4030 3953 3878 3805 5016 4917 4820 4725 4632 4541 4452 4365 4280 4197 4116 4037 3960 3885 3812 ...

16   Pass
r = 100; n_correct = 4187857; assert(isequal(dots_in_sphere(r),n_correct));

RR(:,:,1) = Columns 1 through 15 30000 29801 29604 29409 29216 29025 28836 28649 28464 28281 28100 27921 27744 27569 27396 29801 29602 29405 29210 29017 28826 28637 28450 28265 28082 27901 27722 27545 27370 27197 29604 29405 29208 29013 28820 28629 28440 28253 28068 27885 27704 27525 27348 27173 27000 29409 29210 29013 28818 28625 28434 28245 28058 27873 27690 27509 27330 27153 26978 26805 29216 29017 28820 28625 28432 28241 28052 27865 27680 27497 27316 27137 26960 26785 26612 29025 28826 28629 28434 28241 28050 27861 27674 27489 27306 27125 26946 26769 26594 26421 28836 28637 28440 28245 28052 27861 27672 27485 27300 27117 26936 26757 26580 26405 26232 28649 28450 28253 28058 27865 27674 27485 27298 27113 26930 26749 26570 26393 26218 26045 28464 28265 28068 27873 27680 27489 27300 27113 26928 26745 26564 26385 26208 26033 25860 28281 28082 27885 27690 27497 27306 27117 26930 26745 26562 26381 26202 26025 25850 25677 28100 27901 27704 27509 27316 27125 26936 26749 26564 26381 26200 26021 25844 25669 25496 27921 27722 27525 27330 27137 26946 26757 26570 26385 26202 26021 25842 25665 25490 25317 27744 27545 27348 27153 26960 26769 26580 26393 26208 26025 25844 25665 25488 25313 25140 27569 27370 27173 26978 26785 26594 26405 26218 26033 25850 25669 25490 25313 25138 24965 27396 27197 27000 26805 26612 26421 26232 26045 25860 25677 25496 25317 25140 24965 24792 27225 27026 26829 26634 26441 26250 26061 25874 25689 25506 25325 25146 24969 24794 24621 27056 26857 26660 26465 26272 26081 25892 25705 25520 25337 25156 24977 24800 24625 24452 26889 26690 26493 26298 26105 25914 25725 25538 25353 25170 24989 24810 24633 24458 24285 26724 26525 26328 26133 25940 25749 25560 25373 25188 25005 24824 24645 24468 24293 24120 26561 26362 26165 25970 25777 25586 25397 25210 25025 24842 24661 24482 24305 24130 23957 26400 26201 26004 25809 25616 25425 25236 25049 24864 24681 24500 24321 24144 23969 23796 26241 26042 25845 25650 25457 25266 25077 24890 24705 24522 24341 24162 23985 23810 23637 26084 25885 25688 25493 25300 25109 24920 24733 24548 24365 24184 24005 23828 23653 23480 25929 25730 25533 25338 25145 24954 24765 24578 24393 24210 24029 23850 23673 23498 23325 25776 25577 25380 25185 24992 24801 24612 24425 24240 24057 23876 23697 23520 23345 23172 25625 25426 25229 25034 24841 24650 24461 24274 24089 23906 23725 23546 23369 23194 23021 25476 25277 25080 24885 24692 24501 24312 24125 23940 23757 23576 23397 23220 23045 22872 25329 25130 24933 24738 24545 24354 24165 23978 23793 23610 23429 23250 23073 22898 22725 25184 24985 24788 24593 24400 24209 24020 23833 23648 23465 23284 23105 22928 22753 22580 25041 24842 24645 24450 24257 24066 23877 23690 23505 23322 23141 22962 22785 22610 22437 24900 24701 24504 24309 24116 23925 23736 23549 23364 23181 23000 22821 22644 22469 22296 24761 24562 24365 24170 23977 23786 23597 23410 23225 23042 22861 22682 22505 22330 22157 24624 24425 24228 24033 23840 23649 23460 23273 23088 22905 22724 22545 22368 22193 22020 24489 24290 24093 23898 23705 23514 23325 23138 22953 22770 22589 22410 22233 22058 21885 24356 24157 23960 23765 23572 23381 23192 23005 22820 22637 22456 22277 22100 21925 21752 24225 24026 23829 23634 23441 23250 23061 22874 22689 22506 22325 22146 21969 21794 21621 24096 23897 23700 23505 23312 23121 22932 22745 22560 22377 22196 22017 21840 21665 21492 23969 23770 23573 23378 23185 22994 22805 22618 22433 22250 22069 21890 21713 21538 21365 23844 23645 23448 23253 23060 22869 22680 22493 22308 22125 21944 21765 21588 21413 21240 23721 23522 23325 23130 22937 22746 22557 22370 22185 22002 21821 21642 21465 21290 21117 23600 23401 23204 23009 22816 22625 22436 22249 22064 21881 21700 21521 21344 21169 20996 23481 23282 23085 22890 22697 22506 22317 22130 21945 21762 21581 21402 21225 21050 20877 23364 23165 22968 22773 22580 22389 22200 22013 21828 21645 21464 21285 21108 20933 20760 23249 23050 22853 22658 22465 22274 22085 21898 21713 21530 21349 21170 20993 20818 20645 23136 22937 22740 22545 22352 22161 21972 21785 21600 21417 21236 21057 20880 20705 20532 23025 22826 22629 22434 22241 22050 21861 21674 21489 21306 21125 20946 20769 20594 20421 22916 22717 22520 22325 22132 21941 21752 21565 21380 21197 21016 20837 20660 20485 20312 22809 22610 22413 22218 22025 21834 21645 21458 21273 21090 20909 20730 20553 20378 20205 22704 22505 22308 22113 21920 21729 21540 21353 21168 20985 20804 20625 20448 20273 20100 22601 22402 22205 22010 21817 21626 21437 21250 21065 20882 20701 20522 20345 20170 19997 22500 22301 22104 21909 21716 21525 21336 21149 20964 20781 20600 20421 20244 20069 19896 22401 22202 22005 21810 21617 21426 21237 21050 20865 20682 20501 20322 20145 19970 19797 22304 22105 21908 21713 21520 21329 21140 20953 20768 20585 20404 20225 20048 19873 19700 22209 22010 21813 21618 21425 21234 21045 20858 20673 20490 20309 20130 19953 19778 19605 22116 21917 21720 21525 21332 21141 20952 20765 20580 20397 20216 20037 19860 19685 19512 2...

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!