clear; clear all; clc; format long; E=200000; mu=0.33; dispbcelements=[409 0 0 0;410 0 0 0;411 0 0 0;412 0 0 0; 413 0 0 0; 414 0 0 0; 415 0 0 0; 416 0 0 0; 417 0 0 0; 418 0 0 0]; fbcelements=[419 0 0 10000;420 0 0 10000;421 0 0 10000;422 0 0 10000; 423 0 0 10000;424 0 0 10000; 425 0 0 10000; 426 0 0 10000; 427 0 0 10000; 428 0 0 10000]; fbczeroelements(1:408,1:4)=0; for i=1:408 fbczeroelements(i,1)=i; end fbcelements=[fbczeroelements;fbcelements]; xyzofelements=[1 0.000000e+000 0.000000e+000 0.000000e+000; 1 1.333333e+000 0.000000e+000 0.000000e+000; 1 0.000000e+000 0.000000e+000 2.040816e+000; 2 1.333333e+000 0.000000e+000 0.000000e+000; 2 2.666667e+000 0.000000e+000 0.000000e+000; 2 4.000000e+000 0.000000e+000 2.040816e+000; 3 2.666667e+000 0.000000e+000 0.000000e+000; 3 4.000000e+000 0.000000e+000 0.000000e+000; 3 4.000000e+000 0.000000e+000 2.040816e+000; 4 4.000000e+000 0.000000e+000 1.000000e+002; 4 2.666667e+000 0.000000e+000 1.000000e+002; 4 4.000000e+000 0.000000e+000 9.795918e+001; 5 2.666667e+000 0.000000e+000 1.000000e+002; 5 1.333333e+000 0.000000e+000 1.000000e+002; 5 4.000000e+000 0.000000e+000 9.795918e+001; 6 1.333333e+000 0.000000e+000 1.000000e+002; 6 0.000000e+000 0.000000e+000 1.000000e+002; 6 0.000000e+000 0.000000e+000 9.795918e+001; 7 0.000000e+000 0.000000e+000 9.795918e+001; 7 0.000000e+000 0.000000e+000 9.591837e+001; 7 1.333333e+000 0.000000e+000 1.000000e+002; 8 0.000000e+000 0.000000e+000 9.591837e+001; 8 0.000000e+000 0.000000e+000 9.387755e+001; 8 4.000000e+000 0.000000e+000 9.591837e+001; 9 0.000000e+000 0.000000e+000 9.387755e+001; 9 0.000000e+000 0.000000e+000 9.183673e+001; 9 4.000000e+000 0.000000e+000 9.183673e+001; 10 0.000000e+000 0.000000e+000 9.183673e+001; 10 0.000000e+000 0.000000e+000 8.979592e+001; 10 4.000000e+000 0.000000e+000 8.979592e+001; 11 0.000000e+000 0.000000e+000 8.979592e+001; 11 0.000000e+000 0.000000e+000 8.775510e+001; 11 4.000000e+000 0.000000e+000 8.775510e+001; 12 0.000000e+000 0.000000e+000 8.775510e+001; 12 0.000000e+000 0.000000e+000 8.571429e+001; 12 4.000000e+000 0.000000e+000 8.571429e+001; 13 0.000000e+000 0.000000e+000 8.571429e+001; 13 0.000000e+000 0.000000e+000 8.367347e+001; 13 4.000000e+000 0.000000e+000 8.367347e+001; 14 0.000000e+000 0.000000e+000 8.367347e+001; 14 0.000000e+000 0.000000e+000 8.163265e+001; 14 4.000000e+000 0.000000e+000 8.163265e+001; 15 0.000000e+000 0.000000e+000 8.163265e+001; 15 0.000000e+000 0.000000e+000 7.959184e+001; 15 4.000000e+000 0.000000e+000 7.959184e+001; 16 0.000000e+000 0.000000e+000 7.959184e+001; 16 0.000000e+000 0.000000e+000 7.755102e+001; 16 4.000000e+000 0.000000e+000 7.755102e+001; 17 0.000000e+000 0.000000e+000 7.755102e+001; 17 0.000000e+000 0.000000e+000 7.551020e+001; 17 4.000000e+000 0.000000e+000 7.551020e+001; 18 0.000000e+000 0.000000e+000 7.551020e+001; 18 0.000000e+000 0.000000e+000 7.346939e+001; 18 4.000000e+000 0.000000e+000 7.346939e+001; 19 0.000000e+000 0.000000e+000 7.346939e+001; 19 0.000000e+000 0.000000e+000 7.142857e+001; 19 4.000000e+000 0.000000e+000 7.142857e+001; 20 0.000000e+000 0.000000e+000 7.142857e+001; 20 0.000000e+000 0.000000e+000 6.938776e+001; 20 4.000000e+000 0.000000e+000 6.938776e+001; 21 0.000000e+000 0.000000e+000 6.938776e+001; 21 0.000000e+000 0.000000e+000 6.734694e+001; 21 4.000000e+000 0.000000e+000 6.734694e+001; 22 0.000000e+000 0.000000e+000 6.734694e+001; 22 0.000000e+000 0.000000e+000 6.530612e+001; 22 4.000000e+000 0.000000e+000 6.530612e+001; 23 0.000000e+000 0.000000e+000 6.530612e+001; 23 0.000000e+000 0.000000e+000 6.326530e+001; 23 4.000000e+000 0.000000e+000 6.326530e+001; 24 0.000000e+000 0.000000e+000 6.326530e+001; 24 0.000000e+000 0.000000e+000 6.122449e+001; 24 4.000000e+000 0.000000e+000 6.122449e+001; 25 0.000000e+000 0.000000e+000 6.122449e+001; 25 0.000000e+000 0.000000e+000 5.918367e+001; 25 4.000000e+000 0.000000e+000 5.918367e+001; 26 0.000000e+000 0.000000e+000 5.918367e+001; 26 0.000000e+000 0.000000e+000 5.714286e+001; 26 4.000000e+000 0.000000e+000 5.714286e+001; 27 0.000000e+000 0.000000e+000 5.714286e+001; 27 0.000000e+000 0.000000e+000 5.510204e+001; 27 4.000000e+000 0.000000e+000 5.510204e+001; 28 0.000000e+000 0.000000e+000 5.510204e+001; 28 0.000000e+000 0.000000e+000 5.306123e+001; 28 4.000000e+000 0.000000e+000 5.306123e+001; 29 0.000000e+000 0.000000e+000 5.306123e+001; 29 0.000000e+000 0.000000e+000 5.102041e+001; 29 4.000000e+000 0.000000e+000 5.102041e+001; 30 0.000000e+000 0.000000e+000 5.102041e+001; 30 0.000000e+000 0.000000e+000 4.897959e+001; 30 4.000000e+000 0.000000e+000 4.897959e+001; 31 0.000000e+000 0.000000e+000 4.897959e+001; 31 0.000000e+000 0.000000e+000 4.693877e+001; 31 4.000000e+000 0.000000e+000 4.693877e+001; 32 0.000000e+000 0.000000e+000 4.693877e+001; 32 0.000000e+000 0.000000e+000 4.489796e+001; 32 4.000000e+000 0.000000e+000 4.489796e+001; 33 0.000000e+000 0.000000e+000 4.489796e+001; 33 0.000000e+000 0.000000e+000 4.285714e+001; 33 4.000000e+000 0.000000e+000 4.489796e+001; 34 0.000000e+000 0.000000e+000 4.285714e+001; 34 0.000000e+000 0.000000e+000 4.081633e+001; 34 4.000000e+000 0.000000e+000 4.285714e+001; 35 0.000000e+000 0.000000e+000 4.081633e+001; 35 0.000000e+000 0.000000e+000 3.877551e+001; 35 4.000000e+000 0.000000e+000 4.081633e+001; 36 0.000000e+000 0.000000e+000 3.877551e+001; 36 0.000000e+000 0.000000e+000 3.673470e+001; 36 4.000000e+000 0.000000e+000 3.673470e+001; 37 0.000000e+000 0.000000e+000 3.673470e+001; 37 0.000000e+000 0.000000e+000 3.469388e+001; 37 4.000000e+000 0.000000e+000 3.469388e+001; 38 0.000000e+000 0.000000e+000 3.469388e+001; 38 0.000000e+000 0.000000e+000 3.265306e+001; 38 4.000000e+000 0.000000e+000 3.265306e+001; 39 0.000000e+000 0.000000e+000 3.265306e+001; 39 0.000000e+000 0.000000e+000 3.061225e+001; 39 4.000000e+000 0.000000e+000 3.061225e+001; 40 0.000000e+000 0.000000e+000 3.061225e+001; 40 0.000000e+000 0.000000e+000 2.857143e+001; 40 4.000000e+000 0.000000e+000 2.857143e+001; 41 0.000000e+000 0.000000e+000 2.857143e+001; 41 0.000000e+000 0.000000e+000 2.653061e+001; 41 4.000000e+000 0.000000e+000 2.653061e+001; 42 0.000000e+000 0.000000e+000 2.653061e+001; 42 0.000000e+000 0.000000e+000 2.448980e+001; 42 4.000000e+000 0.000000e+000 2.448980e+001; 43 0.000000e+000 0.000000e+000 2.448980e+001; 43 0.000000e+000 0.000000e+000 2.244898e+001; 43 4.000000e+000 0.000000e+000 2.244898e+001; 44 0.000000e+000 0.000000e+000 2.244898e+001; 44 0.000000e+000 0.000000e+000 2.040816e+001; 44 4.000000e+000 0.000000e+000 2.244898e+001; 45 0.000000e+000 0.000000e+000 2.040816e+001; 45 0.000000e+000 0.000000e+000 1.836735e+001; 45 4.000000e+000 0.000000e+000 2.040816e+001; 46 0.000000e+000 0.000000e+000 1.836735e+001; 46 0.000000e+000 0.000000e+000 1.632653e+001; 46 4.000000e+000 0.000000e+000 1.632653e+001; 47 0.000000e+000 0.000000e+000 1.632653e+001; 47 0.000000e+000 0.000000e+000 1.428571e+001; 47 4.000000e+000 0.000000e+000 1.428571e+001; 48 0.000000e+000 0.000000e+000 1.428571e+001; 48 0.000000e+000 0.000000e+000 1.224490e+001; 48 4.000000e+000 0.000000e+000 1.224490e+001; 49 0.000000e+000 0.000000e+000 1.224490e+001; 49 0.000000e+000 0.000000e+000 1.020408e+001; 49 4.000000e+000 0.000000e+000 1.224490e+001; 50 0.000000e+000 0.000000e+000 1.020408e+001; 50 0.000000e+000 0.000000e+000 8.163265e+000; 50 4.000000e+000 0.000000e+000 1.020408e+001; 51 0.000000e+000 0.000000e+000 8.163265e+000; 51 0.000000e+000 0.000000e+000 6.122449e+000; 51 4.000000e+000 0.000000e+000 6.122449e+000; 52 0.000000e+000 0.000000e+000 6.122449e+000; 52 0.000000e+000 0.000000e+000 4.081633e+000; 52 4.000000e+000 0.000000e+000 6.122449e+000; 53 0.000000e+000 0.000000e+000 4.081633e+000; 53 0.000000e+000 0.000000e+000 2.040816e+000; 53 1.333333e+000 0.000000e+000 0.000000e+000; 54 4.000000e+000 0.000000e+000 2.040816e+000; 54 4.000000e+000 0.000000e+000 4.081633e+000; 54 0.000000e+000 0.000000e+000 4.081633e+000; 55 4.000000e+000 0.000000e+000 4.081633e+000; 55 4.000000e+000 0.000000e+000 6.122449e+000; 55 0.000000e+000 0.000000e+000 4.081633e+000; 56 4.000000e+000 0.000000e+000 6.122449e+000; 56 4.000000e+000 0.000000e+000 8.163265e+000; 56 0.000000e+000 0.000000e+000 8.163265e+000; 57 4.000000e+000 0.000000e+000 8.163265e+000; 57 4.000000e+000 0.000000e+000 1.020408e+001; 57 0.000000e+000 0.000000e+000 8.163265e+000; 58 4.000000e+000 0.000000e+000 1.020408e+001; 58 4.000000e+000 0.000000e+000 1.224490e+001; 58 0.000000e+000 0.000000e+000 1.020408e+001; 59 4.000000e+000 0.000000e+000 1.224490e+001; 59 4.000000e+000 0.000000e+000 1.428571e+001; 59 0.000000e+000 0.000000e+000 1.428571e+001; 60 4.000000e+000 0.000000e+000 1.428571e+001; 60 4.000000e+000 0.000000e+000 1.632653e+001; 60 0.000000e+000 0.000000e+000 1.632653e+001; 61 4.000000e+000 0.000000e+000 1.632653e+001; 61 4.000000e+000 0.000000e+000 1.836735e+001; 61 0.000000e+000 0.000000e+000 1.836735e+001; 62 4.000000e+000 0.000000e+000 1.836735e+001; 62 4.000000e+000 0.000000e+000 2.040816e+001; 62 0.000000e+000 0.000000e+000 1.836735e+001; 63 4.000000e+000 0.000000e+000 2.040816e+001; 63 4.000000e+000 0.000000e+000 2.244898e+001; 63 0.000000e+000 0.000000e+000 2.040816e+001; 64 4.000000e+000 0.000000e+000 2.244898e+001; 64 4.000000e+000 0.000000e+000 2.448980e+001; 64 0.000000e+000 0.000000e+000 2.448980e+001; 65 4.000000e+000 0.000000e+000 2.448980e+001; 65 4.000000e+000 0.000000e+000 2.653061e+001; 65 0.000000e+000 0.000000e+000 2.653061e+001; 66 4.000000e+000 0.000000e+000 2.653061e+001; 66 4.000000e+000 0.000000e+000 2.857143e+001; 66 0.000000e+000 0.000000e+000 2.857143e+001; 67 4.000000e+000 0.000000e+000 2.857143e+001; 67 4.000000e+000 0.000000e+000 3.061225e+001; 67 0.000000e+000 0.000000e+000 3.061225e+001; 68 4.000000e+000 0.000000e+000 3.061225e+001; 68 4.000000e+000 0.000000e+000 3.265306e+001; 68 0.000000e+000 0.000000e+000 3.265306e+001; 69 4.000000e+000 0.000000e+000 3.265306e+001; 69 4.000000e+000 0.000000e+000 3.469388e+001; 69 0.000000e+000 0.000000e+000 3.469388e+001; 70 4.000000e+000 0.000000e+000 3.469388e+001; 70 4.000000e+000 0.000000e+000 3.673470e+001; 70 0.000000e+000 0.000000e+000 3.673470e+001; 71 4.000000e+000 0.000000e+000 3.673470e+001; 71 4.000000e+000 0.000000e+000 3.877551e+001; 71 0.000000e+000 0.000000e+000 3.877551e+001; 72 4.000000e+000 0.000000e+000 3.877551e+001; 72 4.000000e+000 0.000000e+000 4.081633e+001; 72 0.000000e+000 0.000000e+000 3.877551e+001; 73 4.000000e+000 0.000000e+000 4.081633e+001; 73 4.000000e+000 0.000000e+000 4.285714e+001; 73 0.000000e+000 0.000000e+000 4.081633e+001; 74 4.000000e+000 0.000000e+000 4.285714e+001; 74 4.000000e+000 0.000000e+000 4.489796e+001; 74 0.000000e+000 0.000000e+000 4.285714e+001; 75 4.000000e+000 0.000000e+000 4.489796e+001; 75 4.000000e+000 0.000000e+000 4.693877e+001; 75 0.000000e+000 0.000000e+000 4.693877e+001; 76 4.000000e+000 0.000000e+000 4.693877e+001; 76 4.000000e+000 0.000000e+000 4.897959e+001; 76 0.000000e+000 0.000000e+000 4.897959e+001; 77 4.000000e+000 0.000000e+000 4.897959e+001; 77 4.000000e+000 0.000000e+000 5.102041e+001; 77 0.000000e+000 0.000000e+000 5.102041e+001; 78 4.000000e+000 0.000000e+000 5.102041e+001; 78 4.000000e+000 0.000000e+000 5.306123e+001; 78 0.000000e+000 0.000000e+000 5.306123e+001; 79 4.000000e+000 0.000000e+000 5.306123e+001; 79 4.000000e+000 0.000000e+000 5.510204e+001; 79 0.000000e+000 0.000000e+000 5.510204e+001; 80 4.000000e+000 0.000000e+000 5.510204e+001; 80 4.000000e+000 0.000000e+000 5.714286e+001; 80 0.000000e+000 0.000000e+000 5.714286e+001; 81 4.000000e+000 0.000000e+000 5.714286e+001; 81 4.000000e+000 0.000000e+000 5.918367e+001; 81 0.000000e+000 0.000000e+000 5.918367e+001; 82 4.000000e+000 0.000000e+000 5.918367e+001; 82 4.000000e+000 0.000000e+000 6.122449e+001; 82 0.000000e+000 0.000000e+000 6.122449e+001; 83 4.000000e+000 0.000000e+000 6.122449e+001; 83 4.000000e+000 0.000000e+000 6.326530e+001; 83 0.000000e+000 0.000000e+000 6.326530e+001; 84 4.000000e+000 0.000000e+000 6.326530e+001; 84 4.000000e+000 0.000000e+000 6.530612e+001; 84 0.000000e+000 0.000000e+000 6.530612e+001; 85 4.000000e+000 0.000000e+000 6.530612e+001; 85 4.000000e+000 0.000000e+000 6.734694e+001; 85 0.000000e+000 0.000000e+000 6.734694e+001; 86 4.000000e+000 0.000000e+000 6.734694e+001; 86 4.000000e+000 0.000000e+000 6.938776e+001; 86 0.000000e+000 0.000000e+000 6.938776e+001; 87 4.000000e+000 0.000000e+000 6.938776e+001; 87 4.000000e+000 0.000000e+000 7.142857e+001; 87 0.000000e+000 0.000000e+000 7.142857e+001; 88 4.000000e+000 0.000000e+000 7.142857e+001; 88 4.000000e+000 0.000000e+000 7.346939e+001; 88 0.000000e+000 0.000000e+000 7.346939e+001; 89 4.000000e+000 0.000000e+000 7.346939e+001; 89 4.000000e+000 0.000000e+000 7.551020e+001; 89 0.000000e+000 0.000000e+000 7.551020e+001; 90 4.000000e+000 0.000000e+000 7.551020e+001; 90 4.000000e+000 0.000000e+000 7.755102e+001; 90 0.000000e+000 0.000000e+000 7.755102e+001; 91 4.000000e+000 0.000000e+000 7.755102e+001; 91 4.000000e+000 0.000000e+000 7.959184e+001; 91 0.000000e+000 0.000000e+000 7.959184e+001; 92 4.000000e+000 0.000000e+000 7.959184e+001; 92 4.000000e+000 0.000000e+000 8.163265e+001; 92 0.000000e+000 0.000000e+000 8.163265e+001; 93 4.000000e+000 0.000000e+000 8.163265e+001; 93 4.000000e+000 0.000000e+000 8.367347e+001; 93 0.000000e+000 0.000000e+000 8.367347e+001; 94 4.000000e+000 0.000000e+000 8.367347e+001; 94 4.000000e+000 0.000000e+000 8.571429e+001; 94 0.000000e+000 0.000000e+000 8.571429e+001; 95 4.000000e+000 0.000000e+000 8.571429e+001; 95 4.000000e+000 0.000000e+000 8.775510e+001; 95 0.000000e+000 0.000000e+000 8.775510e+001; 96 4.000000e+000 0.000000e+000 8.775510e+001; 96 4.000000e+000 0.000000e+000 8.979592e+001; 96 0.000000e+000 0.000000e+000 8.979592e+001; 97 4.000000e+000 0.000000e+000 8.979592e+001; 97 4.000000e+000 0.000000e+000 9.183673e+001; 97 0.000000e+000 0.000000e+000 9.183673e+001; 98 4.000000e+000 0.000000e+000 9.183673e+001; 98 4.000000e+000 0.000000e+000 9.387755e+001; 98 0.000000e+000 0.000000e+000 9.387755e+001; 99 4.000000e+000 0.000000e+000 9.387755e+001; 99 4.000000e+000 0.000000e+000 9.591837e+001; 99 0.000000e+000 0.000000e+000 9.387755e+001; 100 4.000000e+000 0.000000e+000 9.591837e+001; 100 4.000000e+000 0.000000e+000 9.795918e+001; 100 0.000000e+000 0.000000e+000 9.591837e+001; 101 1.333333e+000 0.000000e+000 0.000000e+000; 101 4.000000e+000 0.000000e+000 2.040816e+000; 101 0.000000e+000 0.000000e+000 4.081633e+000; 102 4.000000e+000 0.000000e+000 9.795918e+001; 102 1.333333e+000 0.000000e+000 1.000000e+002; 102 0.000000e+000 0.000000e+000 9.591837e+001; 103 4.000000e+000 0.000000e+000 0.000000e+000; 103 4.000000e+000 1.333333e+000 0.000000e+000; 103 4.000000e+000 0.000000e+000 2.040816e+000; 104 4.000000e+000 1.333333e+000 0.000000e+000; 104 4.000000e+000 2.666667e+000 0.000000e+000; 104 4.000000e+000 4.000000e+000 2.040816e+000; 105 4.000000e+000 2.666667e+000 0.000000e+000; 105 4.000000e+000 4.000000e+000 0.000000e+000; 105 4.000000e+000 4.000000e+000 2.040816e+000; 106 4.000000e+000 4.000000e+000 1.000000e+002; 106 4.000000e+000 2.666667e+000 1.000000e+002; 106 4.000000e+000 4.000000e+000 9.795918e+001; 107 4.000000e+000 2.666667e+000 1.000000e+002; 107 4.000000e+000 1.333333e+000 1.000000e+002; 107 4.000000e+000 4.000000e+000 9.795918e+001; 108 4.000000e+000 1.333333e+000 1.000000e+002; 108 4.000000e+000 0.000000e+000 1.000000e+002; 108 4.000000e+000 0.000000e+000 9.795918e+001; 109 4.000000e+000 0.000000e+000 9.795918e+001; 109 4.000000e+000 0.000000e+000 9.591837e+001; 109 4.000000e+000 1.333333e+000 1.000000e+002; 110 4.000000e+000 0.000000e+000 9.591837e+001; 110 4.000000e+000 0.000000e+000 9.387755e+001; 110 4.000000e+000 4.000000e+000 9.591837e+001; 111 4.000000e+000 0.000000e+000 9.387755e+001; 111 4.000000e+000 0.000000e+000 9.183673e+001; 111 4.000000e+000 4.000000e+000 9.183673e+001; 112 4.000000e+000 0.000000e+000 9.183673e+001; 112 4.000000e+000 0.000000e+000 8.979592e+001; 112 4.000000e+000 4.000000e+000 8.979592e+001; 113 4.000000e+000 0.000000e+000 8.979592e+001; 113 4.000000e+000 0.000000e+000 8.775510e+001; 113 4.000000e+000 4.000000e+000 8.775510e+001; 114 4.000000e+000 0.000000e+000 8.775510e+001; 114 4.000000e+000 0.000000e+000 8.571429e+001; 114 4.000000e+000 4.000000e+000 8.571429e+001; 115 4.000000e+000 0.000000e+000 8.571429e+001; 115 4.000000e+000 0.000000e+000 8.367347e+001; 115 4.000000e+000 4.000000e+000 8.367347e+001; 116 4.000000e+000 0.000000e+000 8.367347e+001; 116 4.000000e+000 0.000000e+000 8.163265e+001; 116 4.000000e+000 4.000000e+000 8.163265e+001; 117 4.000000e+000 0.000000e+000 8.163265e+001; 117 4.000000e+000 0.000000e+000 7.959184e+001; 117 4.000000e+000 4.000000e+000 7.959184e+001; 118 4.000000e+000 0.000000e+000 7.959184e+001; 118 4.000000e+000 0.000000e+000 7.755102e+001; 118 4.000000e+000 4.000000e+000 7.755102e+001; 119 4.000000e+000 0.000000e+000 7.755102e+001; 119 4.000000e+000 0.000000e+000 7.551020e+001; 119 4.000000e+000 4.000000e+000 7.551020e+001; 120 4.000000e+000 0.000000e+000 7.551020e+001; 120 4.000000e+000 0.000000e+000 7.346939e+001; 120 4.000000e+000 4.000000e+000 7.346939e+001; 121 4.000000e+000 0.000000e+000 7.346939e+001; 121 4.000000e+000 0.000000e+000 7.142857e+001; 121 4.000000e+000 4.000000e+000 7.142857e+001; 122 4.000000e+000 0.000000e+000 7.142857e+001; 122 4.000000e+000 0.000000e+000 6.938776e+001; 122 4.000000e+000 4.000000e+000 6.938776e+001; 123 4.000000e+000 0.000000e+000 6.938776e+001; 123 4.000000e+000 0.000000e+000 6.734694e+001; 123 4.000000e+000 4.000000e+000 6.734694e+001; 124 4.000000e+000 0.000000e+000 6.734694e+001; 124 4.000000e+000 0.000000e+000 6.530612e+001; 124 4.000000e+000 4.000000e+000 6.530612e+001; 125 4.000000e+000 0.000000e+000 6.530612e+001; 125 4.000000e+000 0.000000e+000 6.326530e+001; 125 4.000000e+000 4.000000e+000 6.326530e+001; 126 4.000000e+000 0.000000e+000 6.326530e+001; 126 4.000000e+000 0.000000e+000 6.122449e+001; 126 4.000000e+000 4.000000e+000 6.122449e+001; 127 4.000000e+000 0.000000e+000 6.122449e+001; 127 4.000000e+000 0.000000e+000 5.918367e+001; 127 4.000000e+000 4.000000e+000 5.918367e+001; 128 4.000000e+000 0.000000e+000 5.918367e+001; 128 4.000000e+000 0.000000e+000 5.714286e+001; 128 4.000000e+000 4.000000e+000 5.714286e+001; 129 4.000000e+000 0.000000e+000 5.714286e+001; 129 4.000000e+000 0.000000e+000 5.510204e+001; 129 4.000000e+000 4.000000e+000 5.510204e+001; 130 4.000000e+000 0.000000e+000 5.510204e+001; 130 4.000000e+000 0.000000e+000 5.306123e+001; 130 4.000000e+000 4.000000e+000 5.306123e+001; 131 4.000000e+000 0.000000e+000 5.306123e+001; 131 4.000000e+000 0.000000e+000 5.102041e+001; 131 4.000000e+000 4.000000e+000 5.102041e+001; 132 4.000000e+000 0.000000e+000 5.102041e+001; 132 4.000000e+000 0.000000e+000 4.897959e+001; 132 4.000000e+000 4.000000e+000 4.897959e+001; 133 4.000000e+000 0.000000e+000 4.897959e+001; 133 4.000000e+000 0.000000e+000 4.693877e+001; 133 4.000000e+000 4.000000e+000 4.693877e+001; 134 4.000000e+000 0.000000e+000 4.693877e+001; 134 4.000000e+000 0.000000e+000 4.489796e+001; 134 4.000000e+000 4.000000e+000 4.489796e+001; 135 4.000000e+000 0.000000e+000 4.489796e+001; 135 4.000000e+000 0.000000e+000 4.285714e+001; 135 4.000000e+000 4.000000e+000 4.489796e+001; 136 4.000000e+000 0.000000e+000 4.285714e+001; 136 4.000000e+000 0.000000e+000 4.081633e+001; 136 4.000000e+000 4.000000e+000 4.285714e+001; 137 4.000000e+000 0.000000e+000 4.081633e+001; 137 4.000000e+000 0.000000e+000 3.877551e+001; 137 4.000000e+000 4.000000e+000 4.081633e+001; 138 4.000000e+000 0.000000e+000 3.877551e+001; 138 4.000000e+000 0.000000e+000 3.673470e+001; 138 4.000000e+000 4.000000e+000 3.673470e+001; 139 4.000000e+000 0.000000e+000 3.673470e+001; 139 4.000000e+000 0.000000e+000 3.469388e+001; 139 4.000000e+000 4.000000e+000 3.469388e+001; 140 4.000000e+000 0.000000e+000 3.469388e+001; 140 4.000000e+000 0.000000e+000 3.265306e+001; 140 4.000000e+000 4.000000e+000 3.265306e+001; 141 4.000000e+000 0.000000e+000 3.265306e+001; 141 4.000000e+000 0.000000e+000 3.061225e+001; 141 4.000000e+000 4.000000e+000 3.061225e+001; 142 4.000000e+000 0.000000e+000 3.061225e+001; 142 4.000000e+000 0.000000e+000 2.857143e+001; 142 4.000000e+000 4.000000e+000 2.857143e+001; 143 4.000000e+000 0.000000e+000 2.857143e+001; 143 4.000000e+000 0.000000e+000 2.653061e+001; 143 4.000000e+000 4.000000e+000 2.653061e+001; 144 4.000000e+000 0.000000e+000 2.653061e+001; 144 4.000000e+000 0.000000e+000 2.448980e+001; 144 4.000000e+000 4.000000e+000 2.448980e+001; 145 4.000000e+000 0.000000e+000 2.448980e+001; 145 4.000000e+000 0.000000e+000 2.244898e+001; 145 4.000000e+000 4.000000e+000 2.244898e+001; 146 4.000000e+000 0.000000e+000 2.244898e+001; 146 4.000000e+000 0.000000e+000 2.040816e+001; 146 4.000000e+000 4.000000e+000 2.244898e+001; 147 4.000000e+000 0.000000e+000 2.040816e+001; 147 4.000000e+000 0.000000e+000 1.836735e+001; 147 4.000000e+000 4.000000e+000 2.040816e+001; 148 4.000000e+000 0.000000e+000 1.836735e+001; 148 4.000000e+000 0.000000e+000 1.632653e+001; 148 4.000000e+000 4.000000e+000 1.632653e+001; 149 4.000000e+000 0.000000e+000 1.632653e+001; 149 4.000000e+000 0.000000e+000 1.428571e+001; 149 4.000000e+000 4.000000e+000 1.428571e+001; 150 4.000000e+000 0.000000e+000 1.428571e+001; 150 4.000000e+000 0.000000e+000 1.224490e+001; 150 4.000000e+000 4.000000e+000 1.224490e+001; 151 4.000000e+000 0.000000e+000 1.224490e+001; 151 4.000000e+000 0.000000e+000 1.020408e+001; 151 4.000000e+000 4.000000e+000 1.224490e+001; 152 4.000000e+000 0.000000e+000 1.020408e+001; 152 4.000000e+000 0.000000e+000 8.163265e+000; 152 4.000000e+000 4.000000e+000 1.020408e+001; 153 4.000000e+000 0.000000e+000 8.163265e+000; 153 4.000000e+000 0.000000e+000 6.122449e+000; 153 4.000000e+000 4.000000e+000 6.122449e+000; 154 4.000000e+000 0.000000e+000 6.122449e+000; 154 4.000000e+000 0.000000e+000 4.081633e+000; 154 4.000000e+000 4.000000e+000 6.122449e+000; 155 4.000000e+000 0.000000e+000 4.081633e+000; 155 4.000000e+000 0.000000e+000 2.040816e+000; 155 4.000000e+000 1.333333e+000 0.000000e+000; 156 4.000000e+000 4.000000e+000 2.040816e+000; 156 4.000000e+000 4.000000e+000 4.081633e+000; 156 4.000000e+000 0.000000e+000 4.081633e+000; 157 4.000000e+000 4.000000e+000 4.081633e+000; 157 4.000000e+000 4.000000e+000 6.122449e+000; 157 4.000000e+000 0.000000e+000 4.081633e+000; 158 4.000000e+000 4.000000e+000 6.122449e+000; 158 4.000000e+000 4.000000e+000 8.163265e+000; 158 4.000000e+000 0.000000e+000 8.163265e+000; 159 4.000000e+000 4.000000e+000 8.163265e+000; 159 4.000000e+000 4.000000e+000 1.020408e+001; 159 4.000000e+000 0.000000e+000 8.163265e+000; 160 4.000000e+000 4.000000e+000 1.020408e+001; 160 4.000000e+000 4.000000e+000 1.224490e+001; 160 4.000000e+000 0.000000e+000 1.020408e+001; 161 4.000000e+000 4.000000e+000 1.224490e+001; 161 4.000000e+000 4.000000e+000 1.428571e+001; 161 4.000000e+000 0.000000e+000 1.428571e+001; 162 4.000000e+000 4.000000e+000 1.428571e+001; 162 4.000000e+000 4.000000e+000 1.632653e+001; 162 4.000000e+000 0.000000e+000 1.632653e+001; 163 4.000000e+000 4.000000e+000 1.632653e+001; 163 4.000000e+000 4.000000e+000 1.836735e+001; 163 4.000000e+000 0.000000e+000 1.836735e+001; 164 4.000000e+000 4.000000e+000 1.836735e+001; 164 4.000000e+000 4.000000e+000 2.040816e+001; 164 4.000000e+000 0.000000e+000 1.836735e+001; 165 4.000000e+000 4.000000e+000 2.040816e+001; 165 4.000000e+000 4.000000e+000 2.244898e+001; 165 4.000000e+000 0.000000e+000 2.040816e+001; 166 4.000000e+000 4.000000e+000 2.244898e+001; 166 4.000000e+000 4.000000e+000 2.448980e+001; 166 4.000000e+000 0.000000e+000 2.448980e+001; 167 4.000000e+000 4.000000e+000 2.448980e+001; 167 4.000000e+000 4.000000e+000 2.653061e+001; 167 4.000000e+000 0.000000e+000 2.653061e+001; 168 4.000000e+000 4.000000e+000 2.653061e+001; 168 4.000000e+000 4.000000e+000 2.857143e+001; 168 4.000000e+000 0.000000e+000 2.857143e+001; 169 4.000000e+000 4.000000e+000 2.857143e+001; 169 4.000000e+000 4.000000e+000 3.061225e+001; 169 4.000000e+000 0.000000e+000 3.061225e+001; 170 4.000000e+000 4.000000e+000 3.061225e+001; 170 4.000000e+000 4.000000e+000 3.265306e+001; 170 4.000000e+000 0.000000e+000 3.265306e+001; 171 4.000000e+000 4.000000e+000 3.265306e+001; 171 4.000000e+000 4.000000e+000 3.469388e+001; 171 4.000000e+000 0.000000e+000 3.469388e+001; 172 4.000000e+000 4.000000e+000 3.469388e+001; 172 4.000000e+000 4.000000e+000 3.673470e+001; 172 4.000000e+000 0.000000e+000 3.673470e+001; 173 4.000000e+000 4.000000e+000 3.673470e+001; 173 4.000000e+000 4.000000e+000 3.877551e+001; 173 4.000000e+000 0.000000e+000 3.877551e+001; 174 4.000000e+000 4.000000e+000 3.877551e+001; 174 4.000000e+000 4.000000e+000 4.081633e+001; 174 4.000000e+000 0.000000e+000 3.877551e+001; 175 4.000000e+000 4.000000e+000 4.081633e+001; 175 4.000000e+000 4.000000e+000 4.285714e+001; 175 4.000000e+000 0.000000e+000 4.081633e+001; 176 4.000000e+000 4.000000e+000 4.285714e+001; 176 4.000000e+000 4.000000e+000 4.489796e+001; 176 4.000000e+000 0.000000e+000 4.285714e+001; 177 4.000000e+000 4.000000e+000 4.489796e+001; 177 4.000000e+000 4.000000e+000 4.693877e+001; 177 4.000000e+000 0.000000e+000 4.693877e+001; 178 4.000000e+000 4.000000e+000 4.693877e+001; 178 4.000000e+000 4.000000e+000 4.897959e+001; 178 4.000000e+000 0.000000e+000 4.897959e+001; 179 4.000000e+000 4.000000e+000 4.897959e+001; 179 4.000000e+000 4.000000e+000 5.102041e+001; 179 4.000000e+000 0.000000e+000 5.102041e+001; 180 4.000000e+000 4.000000e+000 5.102041e+001; 180 4.000000e+000 4.000000e+000 5.306123e+001; 180 4.000000e+000 0.000000e+000 5.306123e+001; 181 4.000000e+000 4.000000e+000 5.306123e+001; 181 4.000000e+000 4.000000e+000 5.510204e+001; 181 4.000000e+000 0.000000e+000 5.510204e+001; 182 4.000000e+000 4.000000e+000 5.510204e+001; 182 4.000000e+000 4.000000e+000 5.714286e+001; 182 4.000000e+000 0.000000e+000 5.714286e+001; 183 4.000000e+000 4.000000e+000 5.714286e+001; 183 4.000000e+000 4.000000e+000 5.918367e+001; 183 4.000000e+000 0.000000e+000 5.918367e+001; 184 4.000000e+000 4.000000e+000 5.918367e+001; 184 4.000000e+000 4.000000e+000 6.122449e+001; 184 4.000000e+000 0.000000e+000 6.122449e+001; 185 4.000000e+000 4.000000e+000 6.122449e+001; 185 4.000000e+000 4.000000e+000 6.326530e+001; 185 4.000000e+000 0.000000e+000 6.326530e+001; 186 4.000000e+000 4.000000e+000 6.326530e+001; 186 4.000000e+000 4.000000e+000 6.530612e+001; 186 4.000000e+000 0.000000e+000 6.530612e+001; 187 4.000000e+000 4.000000e+000 6.530612e+001; 187 4.000000e+000 4.000000e+000 6.734694e+001; 187 4.000000e+000 0.000000e+000 6.734694e+001; 188 4.000000e+000 4.000000e+000 6.734694e+001; 188 4.000000e+000 4.000000e+000 6.938776e+001; 188 4.000000e+000 0.000000e+000 6.938776e+001; 189 4.000000e+000 4.000000e+000 6.938776e+001; 189 4.000000e+000 4.000000e+000 7.142857e+001; 189 4.000000e+000 0.000000e+000 7.142857e+001; 190 4.000000e+000 4.000000e+000 7.142857e+001; 190 4.000000e+000 4.000000e+000 7.346939e+001; 190 4.000000e+000 0.000000e+000 7.346939e+001; 191 4.000000e+000 4.000000e+000 7.346939e+001; 191 4.000000e+000 4.000000e+000 7.551020e+001; 191 4.000000e+000 0.000000e+000 7.551020e+001; 192 4.000000e+000 4.000000e+000 7.551020e+001; 192 4.000000e+000 4.000000e+000 7.755102e+001; 192 4.000000e+000 0.000000e+000 7.755102e+001; 193 4.000000e+000 4.000000e+000 7.755102e+001; 193 4.000000e+000 4.000000e+000 7.959184e+001; 193 4.000000e+000 0.000000e+000 7.959184e+001; 194 4.000000e+000 4.000000e+000 7.959184e+001; 194 4.000000e+000 4.000000e+000 8.163265e+001; 194 4.000000e+000 0.000000e+000 8.163265e+001; 195 4.000000e+000 4.000000e+000 8.163265e+001; 195 4.000000e+000 4.000000e+000 8.367347e+001; 195 4.000000e+000 0.000000e+000 8.367347e+001; 196 4.000000e+000 4.000000e+000 8.367347e+001; 196 4.000000e+000 4.000000e+000 8.571429e+001; 196 4.000000e+000 0.000000e+000 8.571429e+001; 197 4.000000e+000 4.000000e+000 8.571429e+001; 197 4.000000e+000 4.000000e+000 8.775510e+001; 197 4.000000e+000 0.000000e+000 8.775510e+001; 198 4.000000e+000 4.000000e+000 8.775510e+001; 198 4.000000e+000 4.000000e+000 8.979592e+001; 198 4.000000e+000 0.000000e+000 8.979592e+001; 199 4.000000e+000 4.000000e+000 8.979592e+001; 199 4.000000e+000 4.000000e+000 9.183673e+001; 199 4.000000e+000 0.000000e+000 9.183673e+001; 200 4.000000e+000 4.000000e+000 9.183673e+001; 200 4.000000e+000 4.000000e+000 9.387755e+001; 200 4.000000e+000 0.000000e+000 9.387755e+001; 201 4.000000e+000 4.000000e+000 9.387755e+001; 201 4.000000e+000 4.000000e+000 9.591837e+001; 201 4.000000e+000 0.000000e+000 9.387755e+001; 202 4.000000e+000 4.000000e+000 9.591837e+001; 202 4.000000e+000 4.000000e+000 9.795918e+001; 202 4.000000e+000 0.000000e+000 9.591837e+001; 203 4.000000e+000 1.333333e+000 0.000000e+000; 203 4.000000e+000 4.000000e+000 2.040816e+000; 203 4.000000e+000 0.000000e+000 4.081633e+000; 204 4.000000e+000 4.000000e+000 9.795918e+001; 204 4.000000e+000 1.333333e+000 1.000000e+002; 204 4.000000e+000 0.000000e+000 9.591837e+001; 205 4.000000e+000 4.000000e+000 0.000000e+000; 205 2.666667e+000 4.000000e+000 0.000000e+000; 205 4.000000e+000 4.000000e+000 2.040816e+000; 206 2.666667e+000 4.000000e+000 0.000000e+000; 206 1.333333e+000 4.000000e+000 0.000000e+000; 206 0.000000e+000 4.000000e+000 2.040816e+000; 207 1.333333e+000 4.000000e+000 0.000000e+000; 207 0.000000e+000 4.000000e+000 0.000000e+000; 207 0.000000e+000 4.000000e+000 2.040816e+000; 208 0.000000e+000 4.000000e+000 1.000000e+002; 208 1.333333e+000 4.000000e+000 1.000000e+002; 208 0.000000e+000 4.000000e+000 9.795918e+001; 209 1.333333e+000 4.000000e+000 1.000000e+002; 209 2.666667e+000 4.000000e+000 1.000000e+002; 209 0.000000e+000 4.000000e+000 9.795918e+001; 210 2.666667e+000 4.000000e+000 1.000000e+002; 210 4.000000e+000 4.000000e+000 1.000000e+002; 210 4.000000e+000 4.000000e+000 9.795918e+001; 211 4.000000e+000 4.000000e+000 9.795918e+001; 211 4.000000e+000 4.000000e+000 9.591837e+001; 211 2.666667e+000 4.000000e+000 1.000000e+002; 212 4.000000e+000 4.000000e+000 9.591837e+001; 212 4.000000e+000 4.000000e+000 9.387755e+001; 212 0.000000e+000 4.000000e+000 9.591837e+001; 213 4.000000e+000 4.000000e+000 9.387755e+001; 213 4.000000e+000 4.000000e+000 9.183673e+001; 213 0.000000e+000 4.000000e+000 9.183673e+001; 214 4.000000e+000 4.000000e+000 9.183673e+001; 214 4.000000e+000 4.000000e+000 8.979592e+001; 214 0.000000e+000 4.000000e+000 8.979592e+001; 215 4.000000e+000 4.000000e+000 8.979592e+001; 215 4.000000e+000 4.000000e+000 8.775510e+001; 215 0.000000e+000 4.000000e+000 8.775510e+001; 216 4.000000e+000 4.000000e+000 8.775510e+001; 216 4.000000e+000 4.000000e+000 8.571429e+001; 216 0.000000e+000 4.000000e+000 8.571429e+001; 217 4.000000e+000 4.000000e+000 8.571429e+001; 217 4.000000e+000 4.000000e+000 8.367347e+001; 217 0.000000e+000 4.000000e+000 8.367347e+001; 218 4.000000e+000 4.000000e+000 8.367347e+001; 218 4.000000e+000 4.000000e+000 8.163265e+001; 218 0.000000e+000 4.000000e+000 8.163265e+001; 219 4.000000e+000 4.000000e+000 8.163265e+001; 219 4.000000e+000 4.000000e+000 7.959184e+001; 219 0.000000e+000 4.000000e+000 7.959184e+001; 220 4.000000e+000 4.000000e+000 7.959184e+001; 220 4.000000e+000 4.000000e+000 7.755102e+001; 220 0.000000e+000 4.000000e+000 7.755102e+001; 221 4.000000e+000 4.000000e+000 7.755102e+001; 221 4.000000e+000 4.000000e+000 7.551020e+001; 221 0.000000e+000 4.000000e+000 7.551020e+001; 222 4.000000e+000 4.000000e+000 7.551020e+001; 222 4.000000e+000 4.000000e+000 7.346939e+001; 222 0.000000e+000 4.000000e+000 7.346939e+001; 223 4.000000e+000 4.000000e+000 7.346939e+001; 223 4.000000e+000 4.000000e+000 7.142857e+001; 223 0.000000e+000 4.000000e+000 7.142857e+001; 224 4.000000e+000 4.000000e+000 7.142857e+001; 224 4.000000e+000 4.000000e+000 6.938776e+001; 224 0.000000e+000 4.000000e+000 6.938776e+001; 225 4.000000e+000 4.000000e+000 6.938776e+001; 225 4.000000e+000 4.000000e+000 6.734694e+001; 225 0.000000e+000 4.000000e+000 6.734694e+001; 226 4.000000e+000 4.000000e+000 6.734694e+001; 226 4.000000e+000 4.000000e+000 6.530612e+001; 226 0.000000e+000 4.000000e+000 6.530612e+001; 227 4.000000e+000 4.000000e+000 6.530612e+001; 227 4.000000e+000 4.000000e+000 6.326530e+001; 227 0.000000e+000 4.000000e+000 6.326530e+001; 228 4.000000e+000 4.000000e+000 6.326530e+001; 228 4.000000e+000 4.000000e+000 6.122449e+001; 228 0.000000e+000 4.000000e+000 6.122449e+001; 229 4.000000e+000 4.000000e+000 6.122449e+001; 229 4.000000e+000 4.000000e+000 5.918367e+001; 229 0.000000e+000 4.000000e+000 5.918367e+001; 230 4.000000e+000 4.000000e+000 5.918367e+001; 230 4.000000e+000 4.000000e+000 5.714286e+001; 230 0.000000e+000 4.000000e+000 5.714286e+001; 231 4.000000e+000 4.000000e+000 5.714286e+001; 231 4.000000e+000 4.000000e+000 5.510204e+001; 231 0.000000e+000 4.000000e+000 5.510204e+001; 232 4.000000e+000 4.000000e+000 5.510204e+001; 232 4.000000e+000 4.000000e+000 5.306123e+001; 232 0.000000e+000 4.000000e+000 5.306123e+001; 233 4.000000e+000 4.000000e+000 5.306123e+001; 233 4.000000e+000 4.000000e+000 5.102041e+001; 233 0.000000e+000 4.000000e+000 5.102041e+001; 234 4.000000e+000 4.000000e+000 5.102041e+001; 234 4.000000e+000 4.000000e+000 4.897959e+001; 234 0.000000e+000 4.000000e+000 4.897959e+001; 235 4.000000e+000 4.000000e+000 4.897959e+001; 235 4.000000e+000 4.000000e+000 4.693877e+001; 235 0.000000e+000 4.000000e+000 4.693877e+001; 236 4.000000e+000 4.000000e+000 4.693877e+001; 236 4.000000e+000 4.000000e+000 4.489796e+001; 236 0.000000e+000 4.000000e+000 4.489796e+001; 237 4.000000e+000 4.000000e+000 4.489796e+001; 237 4.000000e+000 4.000000e+000 4.285714e+001; 237 0.000000e+000 4.000000e+000 4.489796e+001; 238 4.000000e+000 4.000000e+000 4.285714e+001; 238 4.000000e+000 4.000000e+000 4.081633e+001; 238 0.000000e+000 4.000000e+000 4.285714e+001; 239 4.000000e+000 4.000000e+000 4.081633e+001; 239 4.000000e+000 4.000000e+000 3.877551e+001; 239 0.000000e+000 4.000000e+000 4.081633e+001; 240 4.000000e+000 4.000000e+000 3.877551e+001; 240 4.000000e+000 4.000000e+000 3.673470e+001; 240 0.000000e+000 4.000000e+000 3.673470e+001; 241 4.000000e+000 4.000000e+000 3.673470e+001; 241 4.000000e+000 4.000000e+000 3.469388e+001; 241 0.000000e+000 4.000000e+000 3.469388e+001; 242 4.000000e+000 4.000000e+000 3.469388e+001; 242 4.000000e+000 4.000000e+000 3.265306e+001; 242 0.000000e+000 4.000000e+000 3.265306e+001; 243 4.000000e+000 4.000000e+000 3.265306e+001; 243 4.000000e+000 4.000000e+000 3.061225e+001; 243 0.000000e+000 4.000000e+000 3.061225e+001; 244 4.000000e+000 4.000000e+000 3.061225e+001; 244 4.000000e+000 4.000000e+000 2.857143e+001; 244 0.000000e+000 4.000000e+000 2.857143e+001; 245 4.000000e+000 4.000000e+000 2.857143e+001; 245 4.000000e+000 4.000000e+000 2.653061e+001; 245 0.000000e+000 4.000000e+000 2.653061e+001; 246 4.000000e+000 4.000000e+000 2.653061e+001; 246 4.000000e+000 4.000000e+000 2.448980e+001; 246 0.000000e+000 4.000000e+000 2.448980e+001; 247 4.000000e+000 4.000000e+000 2.448980e+001; 247 4.000000e+000 4.000000e+000 2.244898e+001; 247 0.000000e+000 4.000000e+000 2.244898e+001; 248 4.000000e+000 4.000000e+000 2.244898e+001; 248 4.000000e+000 4.000000e+000 2.040816e+001; 248 0.000000e+000 4.000000e+000 2.244898e+001; 249 4.000000e+000 4.000000e+000 2.040816e+001; 249 4.000000e+000 4.000000e+000 1.836735e+001; 249 0.000000e+000 4.000000e+000 2.040816e+001; 250 4.000000e+000 4.000000e+000 1.836735e+001; 250 4.000000e+000 4.000000e+000 1.632653e+001; 250 0.000000e+000 4.000000e+000 1.632653e+001; 251 4.000000e+000 4.000000e+000 1.632653e+001; 251 4.000000e+000 4.000000e+000 1.428571e+001; 251 0.000000e+000 4.000000e+000 1.428571e+001; 252 4.000000e+000 4.000000e+000 1.428571e+001; 252 4.000000e+000 4.000000e+000 1.224490e+001; 252 0.000000e+000 4.000000e+000 1.224490e+001; 253 4.000000e+000 4.000000e+000 1.224490e+001; 253 4.000000e+000 4.000000e+000 1.020408e+001; 253 0.000000e+000 4.000000e+000 1.224490e+001; 254 4.000000e+000 4.000000e+000 1.020408e+001; 254 4.000000e+000 4.000000e+000 8.163265e+000; 254 0.000000e+000 4.000000e+000 1.020408e+001; 255 4.000000e+000 4.000000e+000 8.163265e+000; 255 4.000000e+000 4.000000e+000 6.122449e+000; 255 0.000000e+000 4.000000e+000 6.122449e+000; 256 4.000000e+000 4.000000e+000 6.122449e+000; 256 4.000000e+000 4.000000e+000 4.081633e+000; 256 0.000000e+000 4.000000e+000 6.122449e+000; 257 4.000000e+000 4.000000e+000 4.081633e+000; 257 4.000000e+000 4.000000e+000 2.040816e+000; 257 2.666667e+000 4.000000e+000 0.000000e+000; 258 0.000000e+000 4.000000e+000 2.040816e+000; 258 0.000000e+000 4.000000e+000 4.081633e+000; 258 4.000000e+000 4.000000e+000 4.081633e+000; 259 0.000000e+000 4.000000e+000 4.081633e+000; 259 0.000000e+000 4.000000e+000 6.122449e+000; 259 4.000000e+000 4.000000e+000 4.081633e+000; 260 0.000000e+000 4.000000e+000 6.122449e+000; 260 0.000000e+000 4.000000e+000 8.163265e+000; 260 4.000000e+000 4.000000e+000 8.163265e+000; 261 0.000000e+000 4.000000e+000 8.163265e+000; 261 0.000000e+000 4.000000e+000 1.020408e+001; 261 4.000000e+000 4.000000e+000 8.163265e+000; 262 0.000000e+000 4.000000e+000 1.020408e+001; 262 0.000000e+000 4.000000e+000 1.224490e+001; 262 4.000000e+000 4.000000e+000 1.020408e+001; 263 0.000000e+000 4.000000e+000 1.224490e+001; 263 0.000000e+000 4.000000e+000 1.428571e+001; 263 4.000000e+000 4.000000e+000 1.428571e+001; 264 0.000000e+000 4.000000e+000 1.428571e+001; 264 0.000000e+000 4.000000e+000 1.632653e+001; 264 4.000000e+000 4.000000e+000 1.632653e+001; 265 0.000000e+000 4.000000e+000 1.632653e+001; 265 0.000000e+000 4.000000e+000 1.836735e+001; 265 4.000000e+000 4.000000e+000 1.836735e+001; 266 0.000000e+000 4.000000e+000 1.836735e+001; 266 0.000000e+000 4.000000e+000 2.040816e+001; 266 4.000000e+000 4.000000e+000 1.836735e+001; 267 0.000000e+000 4.000000e+000 2.040816e+001; 267 0.000000e+000 4.000000e+000 2.244898e+001; 267 4.000000e+000 4.000000e+000 2.040816e+001; 268 0.000000e+000 4.000000e+000 2.244898e+001; 268 0.000000e+000 4.000000e+000 2.448980e+001; 268 4.000000e+000 4.000000e+000 2.448980e+001; 269 0.000000e+000 4.000000e+000 2.448980e+001; 269 0.000000e+000 4.000000e+000 2.653061e+001; 269 4.000000e+000 4.000000e+000 2.653061e+001; 270 0.000000e+000 4.000000e+000 2.653061e+001; 270 0.000000e+000 4.000000e+000 2.857143e+001; 270 4.000000e+000 4.000000e+000 2.857143e+001; 271 0.000000e+000 4.000000e+000 2.857143e+001; 271 0.000000e+000 4.000000e+000 3.061225e+001; 271 4.000000e+000 4.000000e+000 3.061225e+001; 272 0.000000e+000 4.000000e+000 3.061225e+001; 272 0.000000e+000 4.000000e+000 3.265306e+001; 272 4.000000e+000 4.000000e+000 3.265306e+001; 273 0.000000e+000 4.000000e+000 3.265306e+001; 273 0.000000e+000 4.000000e+000 3.469388e+001; 273 4.000000e+000 4.000000e+000 3.469388e+001; 274 0.000000e+000 4.000000e+000 3.469388e+001; 274 0.000000e+000 4.000000e+000 3.673470e+001; 274 4.000000e+000 4.000000e+000 3.673470e+001; 275 0.000000e+000 4.000000e+000 3.673470e+001; 275 0.000000e+000 4.000000e+000 3.877551e+001; 275 4.000000e+000 4.000000e+000 3.877551e+001; 276 0.000000e+000 4.000000e+000 3.877551e+001; 276 0.000000e+000 4.000000e+000 4.081633e+001; 276 4.000000e+000 4.000000e+000 3.877551e+001; 277 0.000000e+000 4.000000e+000 4.081633e+001; 277 0.000000e+000 4.000000e+000 4.285714e+001; 277 4.000000e+000 4.000000e+000 4.081633e+001; 278 0.000000e+000 4.000000e+000 4.285714e+001; 278 0.000000e+000 4.000000e+000 4.489796e+001; 278 4.000000e+000 4.000000e+000 4.285714e+001; 279 0.000000e+000 4.000000e+000 4.489796e+001; 279 0.000000e+000 4.000000e+000 4.693877e+001; 279 4.000000e+000 4.000000e+000 4.693877e+001; 280 0.000000e+000 4.000000e+000 4.693877e+001; 280 0.000000e+000 4.000000e+000 4.897959e+001; 280 4.000000e+000 4.000000e+000 4.897959e+001; 281 0.000000e+000 4.000000e+000 4.897959e+001; 281 0.000000e+000 4.000000e+000 5.102041e+001; 281 4.000000e+000 4.000000e+000 5.102041e+001; 282 0.000000e+000 4.000000e+000 5.102041e+001; 282 0.000000e+000 4.000000e+000 5.306123e+001; 282 4.000000e+000 4.000000e+000 5.306123e+001; 283 0.000000e+000 4.000000e+000 5.306123e+001; 283 0.000000e+000 4.000000e+000 5.510204e+001; 283 4.000000e+000 4.000000e+000 5.510204e+001; 284 0.000000e+000 4.000000e+000 5.510204e+001; 284 0.000000e+000 4.000000e+000 5.714286e+001; 284 4.000000e+000 4.000000e+000 5.714286e+001; 285 0.000000e+000 4.000000e+000 5.714286e+001; 285 0.000000e+000 4.000000e+000 5.918367e+001; 285 4.000000e+000 4.000000e+000 5.918367e+001; 286 0.000000e+000 4.000000e+000 5.918367e+001; 286 0.000000e+000 4.000000e+000 6.122449e+001; 286 4.000000e+000 4.000000e+000 6.122449e+001; 287 0.000000e+000 4.000000e+000 6.122449e+001; 287 0.000000e+000 4.000000e+000 6.326530e+001; 287 4.000000e+000 4.000000e+000 6.326530e+001; 288 0.000000e+000 4.000000e+000 6.326530e+001; 288 0.000000e+000 4.000000e+000 6.530612e+001; 288 4.000000e+000 4.000000e+000 6.530612e+001; 289 0.000000e+000 4.000000e+000 6.530612e+001; 289 0.000000e+000 4.000000e+000 6.734694e+001; 289 4.000000e+000 4.000000e+000 6.734694e+001; 290 0.000000e+000 4.000000e+000 6.734694e+001; 290 0.000000e+000 4.000000e+000 6.938776e+001; 290 4.000000e+000 4.000000e+000 6.938776e+001; 291 0.000000e+000 4.000000e+000 6.938776e+001; 291 0.000000e+000 4.000000e+000 7.142857e+001; 291 4.000000e+000 4.000000e+000 7.142857e+001; 292 0.000000e+000 4.000000e+000 7.142857e+001; 292 0.000000e+000 4.000000e+000 7.346939e+001; 292 4.000000e+000 4.000000e+000 7.346939e+001; 293 0.000000e+000 4.000000e+000 7.346939e+001; 293 0.000000e+000 4.000000e+000 7.551020e+001; 293 4.000000e+000 4.000000e+000 7.551020e+001; 294 0.000000e+000 4.000000e+000 7.551020e+001; 294 0.000000e+000 4.000000e+000 7.755102e+001; 294 4.000000e+000 4.000000e+000 7.755102e+001; 295 0.000000e+000 4.000000e+000 7.755102e+001; 295 0.000000e+000 4.000000e+000 7.959184e+001; 295 4.000000e+000 4.000000e+000 7.959184e+001; 296 0.000000e+000 4.000000e+000 7.959184e+001; 296 0.000000e+000 4.000000e+000 8.163265e+001; 296 4.000000e+000 4.000000e+000 8.163265e+001; 297 0.000000e+000 4.000000e+000 8.163265e+001; 297 0.000000e+000 4.000000e+000 8.367347e+001; 297 4.000000e+000 4.000000e+000 8.367347e+001; 298 0.000000e+000 4.000000e+000 8.367347e+001; 298 0.000000e+000 4.000000e+000 8.571429e+001; 298 4.000000e+000 4.000000e+000 8.571429e+001; 299 0.000000e+000 4.000000e+000 8.571429e+001; 299 0.000000e+000 4.000000e+000 8.775510e+001; 299 4.000000e+000 4.000000e+000 8.775510e+001; 300 0.000000e+000 4.000000e+000 8.775510e+001; 300 0.000000e+000 4.000000e+000 8.979592e+001; 300 4.000000e+000 4.000000e+000 8.979592e+001; 301 0.000000e+000 4.000000e+000 8.979592e+001; 301 0.000000e+000 4.000000e+000 9.183673e+001; 301 4.000000e+000 4.000000e+000 9.183673e+001; 302 0.000000e+000 4.000000e+000 9.183673e+001; 302 0.000000e+000 4.000000e+000 9.387755e+001; 302 4.000000e+000 4.000000e+000 9.387755e+001; 303 0.000000e+000 4.000000e+000 9.387755e+001; 303 0.000000e+000 4.000000e+000 9.591837e+001; 303 4.000000e+000 4.000000e+000 9.387755e+001; 304 0.000000e+000 4.000000e+000 9.591837e+001; 304 0.000000e+000 4.000000e+000 9.795918e+001; 304 4.000000e+000 4.000000e+000 9.591837e+001; 305 2.666667e+000 4.000000e+000 0.000000e+000; 305 0.000000e+000 4.000000e+000 2.040816e+000; 305 4.000000e+000 4.000000e+000 4.081633e+000; 306 0.000000e+000 4.000000e+000 9.795918e+001; 306 2.666667e+000 4.000000e+000 1.000000e+002; 306 4.000000e+000 4.000000e+000 9.591837e+001; 307 0.000000e+000 4.000000e+000 0.000000e+000; 307 0.000000e+000 2.666667e+000 0.000000e+000; 307 0.000000e+000 4.000000e+000 2.040816e+000; 308 0.000000e+000 2.666667e+000 0.000000e+000; 308 0.000000e+000 1.333333e+000 0.000000e+000; 308 0.000000e+000 0.000000e+000 2.040816e+000; 309 0.000000e+000 1.333333e+000 0.000000e+000; 309 0.000000e+000 0.000000e+000 0.000000e+000; 309 0.000000e+000 0.000000e+000 2.040816e+000; 310 0.000000e+000 0.000000e+000 1.000000e+002; 310 0.000000e+000 1.333333e+000 1.000000e+002; 310 0.000000e+000 0.000000e+000 9.795918e+001; 311 0.000000e+000 1.333333e+000 1.000000e+002; 311 0.000000e+000 2.666667e+000 1.000000e+002; 311 0.000000e+000 0.000000e+000 9.795918e+001; 312 0.000000e+000 2.666667e+000 1.000000e+002; 312 0.000000e+000 4.000000e+000 1.000000e+002; 312 0.000000e+000 4.000000e+000 9.795918e+001; 313 0.000000e+000 0.000000e+000 2.040816e+000; 313 0.000000e+000 0.000000e+000 4.081633e+000; 313 0.000000e+000 4.000000e+000 2.040816e+000; 314 0.000000e+000 0.000000e+000 4.081633e+000; 314 0.000000e+000 0.000000e+000 6.122449e+000; 314 0.000000e+000 4.000000e+000 4.081633e+000; 315 0.000000e+000 0.000000e+000 6.122449e+000; 315 0.000000e+000 0.000000e+000 8.163265e+000; 315 0.000000e+000 4.000000e+000 8.163265e+000; 316 0.000000e+000 0.000000e+000 8.163265e+000; 316 0.000000e+000 0.000000e+000 1.020408e+001; 316 0.000000e+000 4.000000e+000 1.020408e+001; 317 0.000000e+000 0.000000e+000 1.020408e+001; 317 0.000000e+000 0.000000e+000 1.224490e+001; 317 0.000000e+000 4.000000e+000 1.020408e+001; 318 0.000000e+000 0.000000e+000 1.224490e+001; 318 0.000000e+000 0.000000e+000 1.428571e+001; 318 0.000000e+000 4.000000e+000 1.428571e+001; 319 0.000000e+000 0.000000e+000 1.428571e+001; 319 0.000000e+000 0.000000e+000 1.632653e+001; 319 0.000000e+000 4.000000e+000 1.632653e+001; 320 0.000000e+000 0.000000e+000 1.632653e+001; 320 0.000000e+000 0.000000e+000 1.836735e+001; 320 0.000000e+000 4.000000e+000 1.836735e+001; 321 0.000000e+000 0.000000e+000 1.836735e+001; 321 0.000000e+000 0.000000e+000 2.040816e+001; 321 0.000000e+000 4.000000e+000 1.836735e+001; 322 0.000000e+000 0.000000e+000 2.040816e+001; 322 0.000000e+000 0.000000e+000 2.244898e+001; 322 0.000000e+000 4.000000e+000 2.040816e+001; 323 0.000000e+000 0.000000e+000 2.244898e+001; 323 0.000000e+000 0.000000e+000 2.448980e+001; 323 0.000000e+000 4.000000e+000 2.448980e+001; 324 0.000000e+000 0.000000e+000 2.448980e+001; 324 0.000000e+000 0.000000e+000 2.653061e+001; 324 0.000000e+000 4.000000e+000 2.653061e+001; 325 0.000000e+000 0.000000e+000 2.653061e+001; 325 0.000000e+000 0.000000e+000 2.857143e+001; 325 0.000000e+000 4.000000e+000 2.857143e+001; 326 0.000000e+000 0.000000e+000 2.857143e+001; 326 0.000000e+000 0.000000e+000 3.061225e+001; 326 0.000000e+000 4.000000e+000 3.061225e+001; 327 0.000000e+000 0.000000e+000 3.061225e+001; 327 0.000000e+000 0.000000e+000 3.265306e+001; 327 0.000000e+000 4.000000e+000 3.265306e+001; 328 0.000000e+000 0.000000e+000 3.265306e+001; 328 0.000000e+000 0.000000e+000 3.469388e+001; 328 0.000000e+000 4.000000e+000 3.469388e+001; 329 0.000000e+000 0.000000e+000 3.469388e+001; 329 0.000000e+000 0.000000e+000 3.673470e+001; 329 0.000000e+000 4.000000e+000 3.673470e+001; 330 0.000000e+000 0.000000e+000 3.673470e+001; 330 0.000000e+000 0.000000e+000 3.877551e+001; 330 0.000000e+000 4.000000e+000 3.877551e+001; 331 0.000000e+000 0.000000e+000 3.877551e+001; 331 0.000000e+000 0.000000e+000 4.081633e+001; 331 0.000000e+000 4.000000e+000 4.081633e+001; 332 0.000000e+000 0.000000e+000 4.081633e+001; 332 0.000000e+000 0.000000e+000 4.285714e+001; 332 0.000000e+000 4.000000e+000 4.081633e+001; 333 0.000000e+000 0.000000e+000 4.285714e+001; 333 0.000000e+000 0.000000e+000 4.489796e+001; 333 0.000000e+000 4.000000e+000 4.285714e+001; 334 0.000000e+000 0.000000e+000 4.489796e+001; 334 0.000000e+000 0.000000e+000 4.693877e+001; 334 0.000000e+000 4.000000e+000 4.693877e+001; 335 0.000000e+000 0.000000e+000 4.693877e+001; 335 0.000000e+000 0.000000e+000 4.897959e+001; 335 0.000000e+000 4.000000e+000 4.897959e+001; 336 0.000000e+000 0.000000e+000 4.897959e+001; 336 0.000000e+000 0.000000e+000 5.102041e+001; 336 0.000000e+000 4.000000e+000 5.102041e+001; 337 0.000000e+000 0.000000e+000 5.102041e+001; 337 0.000000e+000 0.000000e+000 5.306123e+001; 337 0.000000e+000 4.000000e+000 5.306123e+001; 338 0.000000e+000 0.000000e+000 5.306123e+001; 338 0.000000e+000 0.000000e+000 5.510204e+001; 338 0.000000e+000 4.000000e+000 5.510204e+001; 339 0.000000e+000 0.000000e+000 5.510204e+001; 339 0.000000e+000 0.000000e+000 5.714286e+001; 339 0.000000e+000 4.000000e+000 5.714286e+001; 340 0.000000e+000 0.000000e+000 5.714286e+001; 340 0.000000e+000 0.000000e+000 5.918367e+001; 340 0.000000e+000 4.000000e+000 5.918367e+001; 341 0.000000e+000 0.000000e+000 5.918367e+001; 341 0.000000e+000 0.000000e+000 6.122449e+001; 341 0.000000e+000 4.000000e+000 6.122449e+001; 342 0.000000e+000 0.000000e+000 6.122449e+001; 342 0.000000e+000 0.000000e+000 6.326530e+001; 342 0.000000e+000 4.000000e+000 6.326530e+001; 343 0.000000e+000 0.000000e+000 6.326530e+001; 343 0.000000e+000 0.000000e+000 6.530612e+001; 343 0.000000e+000 4.000000e+000 6.530612e+001; 344 0.000000e+000 0.000000e+000 6.530612e+001; 344 0.000000e+000 0.000000e+000 6.734694e+001; 344 0.000000e+000 4.000000e+000 6.734694e+001; 345 0.000000e+000 0.000000e+000 6.734694e+001; 345 0.000000e+000 0.000000e+000 6.938776e+001; 345 0.000000e+000 4.000000e+000 6.938776e+001; 346 0.000000e+000 0.000000e+000 6.938776e+001; 346 0.000000e+000 0.000000e+000 7.142857e+001; 346 0.000000e+000 4.000000e+000 7.142857e+001; 347 0.000000e+000 0.000000e+000 7.142857e+001; 347 0.000000e+000 0.000000e+000 7.346939e+001; 347 0.000000e+000 4.000000e+000 7.346939e+001; 348 0.000000e+000 0.000000e+000 7.346939e+001; 348 0.000000e+000 0.000000e+000 7.551020e+001; 348 0.000000e+000 4.000000e+000 7.551020e+001; 349 0.000000e+000 0.000000e+000 7.551020e+001; 349 0.000000e+000 0.000000e+000 7.755102e+001; 349 0.000000e+000 4.000000e+000 7.755102e+001; 350 0.000000e+000 0.000000e+000 7.755102e+001; 350 0.000000e+000 0.000000e+000 7.959184e+001; 350 0.000000e+000 4.000000e+000 7.959184e+001; 351 0.000000e+000 0.000000e+000 7.959184e+001; 351 0.000000e+000 0.000000e+000 8.163265e+001; 351 0.000000e+000 4.000000e+000 8.163265e+001; 352 0.000000e+000 0.000000e+000 8.163265e+001; 352 0.000000e+000 0.000000e+000 8.367347e+001; 352 0.000000e+000 4.000000e+000 8.367347e+001; 353 0.000000e+000 0.000000e+000 8.367347e+001; 353 0.000000e+000 0.000000e+000 8.571429e+001; 353 0.000000e+000 4.000000e+000 8.571429e+001; 354 0.000000e+000 0.000000e+000 8.571429e+001; 354 0.000000e+000 0.000000e+000 8.775510e+001; 354 0.000000e+000 4.000000e+000 8.775510e+001; 355 0.000000e+000 0.000000e+000 8.775510e+001; 355 0.000000e+000 0.000000e+000 8.979592e+001; 355 0.000000e+000 4.000000e+000 8.979592e+001; 356 0.000000e+000 0.000000e+000 8.979592e+001; 356 0.000000e+000 0.000000e+000 9.183673e+001; 356 0.000000e+000 4.000000e+000 9.183673e+001; 357 0.000000e+000 0.000000e+000 9.183673e+001; 357 0.000000e+000 0.000000e+000 9.387755e+001; 357 0.000000e+000 4.000000e+000 9.387755e+001; 358 0.000000e+000 0.000000e+000 9.387755e+001; 358 0.000000e+000 0.000000e+000 9.591837e+001; 358 0.000000e+000 4.000000e+000 9.591837e+001; 359 0.000000e+000 0.000000e+000 9.591837e+001; 359 0.000000e+000 0.000000e+000 9.795918e+001; 359 0.000000e+000 4.000000e+000 9.795918e+001; 360 0.000000e+000 4.000000e+000 9.795918e+001; 360 0.000000e+000 4.000000e+000 9.591837e+001; 360 0.000000e+000 0.000000e+000 9.591837e+001; 361 0.000000e+000 4.000000e+000 9.591837e+001; 361 0.000000e+000 4.000000e+000 9.387755e+001; 361 0.000000e+000 0.000000e+000 9.387755e+001; 362 0.000000e+000 4.000000e+000 9.387755e+001; 362 0.000000e+000 4.000000e+000 9.183673e+001; 362 0.000000e+000 0.000000e+000 9.183673e+001; 363 0.000000e+000 4.000000e+000 9.183673e+001; 363 0.000000e+000 4.000000e+000 8.979592e+001; 363 0.000000e+000 0.000000e+000 8.979592e+001; 364 0.000000e+000 4.000000e+000 8.979592e+001; 364 0.000000e+000 4.000000e+000 8.775510e+001; 364 0.000000e+000 0.000000e+000 8.775510e+001; 365 0.000000e+000 4.000000e+000 8.775510e+001; 365 0.000000e+000 4.000000e+000 8.571429e+001; 365 0.000000e+000 0.000000e+000 8.571429e+001; 366 0.000000e+000 4.000000e+000 8.571429e+001; 366 0.000000e+000 4.000000e+000 8.367347e+001; 366 0.000000e+000 0.000000e+000 8.367347e+001; 367 0.000000e+000 4.000000e+000 8.367347e+001; 367 0.000000e+000 4.000000e+000 8.163265e+001; 367 0.000000e+000 0.000000e+000 8.163265e+001; 368 0.000000e+000 4.000000e+000 8.163265e+001; 368 0.000000e+000 4.000000e+000 7.959184e+001; 368 0.000000e+000 0.000000e+000 7.959184e+001; 369 0.000000e+000 4.000000e+000 7.959184e+001; 369 0.000000e+000 4.000000e+000 7.755102e+001; 369 0.000000e+000 0.000000e+000 7.755102e+001; 370 0.000000e+000 4.000000e+000 7.755102e+001; 370 0.000000e+000 4.000000e+000 7.551020e+001; 370 0.000000e+000 0.000000e+000 7.551020e+001; 371 0.000000e+000 4.000000e+000 7.551020e+001; 371 0.000000e+000 4.000000e+000 7.346939e+001; 371 0.000000e+000 0.000000e+000 7.346939e+001; 372 0.000000e+000 4.000000e+000 7.346939e+001; 372 0.000000e+000 4.000000e+000 7.142857e+001; 372 0.000000e+000 0.000000e+000 7.142857e+001; 373 0.000000e+000 4.000000e+000 7.142857e+001; 373 0.000000e+000 4.000000e+000 6.938776e+001; 373 0.000000e+000 0.000000e+000 6.938776e+001; 374 0.000000e+000 4.000000e+000 6.938776e+001; 374 0.000000e+000 4.000000e+000 6.734694e+001; 374 0.000000e+000 0.000000e+000 6.734694e+001; 375 0.000000e+000 4.000000e+000 6.734694e+001; 375 0.000000e+000 4.000000e+000 6.530612e+001; 375 0.000000e+000 0.000000e+000 6.530612e+001; 376 0.000000e+000 4.000000e+000 6.530612e+001; 376 0.000000e+000 4.000000e+000 6.326530e+001; 376 0.000000e+000 0.000000e+000 6.326530e+001; 377 0.000000e+000 4.000000e+000 6.326530e+001; 377 0.000000e+000 4.000000e+000 6.122449e+001; 377 0.000000e+000 0.000000e+000 6.122449e+001; 378 0.000000e+000 4.000000e+000 6.122449e+001; 378 0.000000e+000 4.000000e+000 5.918367e+001; 378 0.000000e+000 0.000000e+000 5.918367e+001; 379 0.000000e+000 4.000000e+000 5.918367e+001; 379 0.000000e+000 4.000000e+000 5.714286e+001; 379 0.000000e+000 0.000000e+000 5.714286e+001; 380 0.000000e+000 4.000000e+000 5.714286e+001; 380 0.000000e+000 4.000000e+000 5.510204e+001; 380 0.000000e+000 0.000000e+000 5.510204e+001; 381 0.000000e+000 4.000000e+000 5.510204e+001; 381 0.000000e+000 4.000000e+000 5.306123e+001; 381 0.000000e+000 0.000000e+000 5.306123e+001; 382 0.000000e+000 4.000000e+000 5.306123e+001; 382 0.000000e+000 4.000000e+000 5.102041e+001; 382 0.000000e+000 0.000000e+000 5.102041e+001; 383 0.000000e+000 4.000000e+000 5.102041e+001; 383 0.000000e+000 4.000000e+000 4.897959e+001; 383 0.000000e+000 0.000000e+000 4.897959e+001; 384 0.000000e+000 4.000000e+000 4.897959e+001; 384 0.000000e+000 4.000000e+000 4.693877e+001; 384 0.000000e+000 0.000000e+000 4.693877e+001; 385 0.000000e+000 4.000000e+000 4.693877e+001; 385 0.000000e+000 4.000000e+000 4.489796e+001; 385 0.000000e+000 0.000000e+000 4.489796e+001; 386 0.000000e+000 4.000000e+000 4.489796e+001; 386 0.000000e+000 4.000000e+000 4.285714e+001; 386 0.000000e+000 0.000000e+000 4.489796e+001; 387 0.000000e+000 4.000000e+000 4.285714e+001; 387 0.000000e+000 4.000000e+000 4.081633e+001; 387 0.000000e+000 0.000000e+000 4.285714e+001; 388 0.000000e+000 4.000000e+000 4.081633e+001; 388 0.000000e+000 4.000000e+000 3.877551e+001; 388 0.000000e+000 0.000000e+000 3.877551e+001; 389 0.000000e+000 4.000000e+000 3.877551e+001; 389 0.000000e+000 4.000000e+000 3.673470e+001; 389 0.000000e+000 0.000000e+000 3.673470e+001; 390 0.000000e+000 4.000000e+000 3.673470e+001; 390 0.000000e+000 4.000000e+000 3.469388e+001; 390 0.000000e+000 0.000000e+000 3.469388e+001; 391 0.000000e+000 4.000000e+000 3.469388e+001; 391 0.000000e+000 4.000000e+000 3.265306e+001; 391 0.000000e+000 0.000000e+000 3.265306e+001; 392 0.000000e+000 4.000000e+000 3.265306e+001; 392 0.000000e+000 4.000000e+000 3.061225e+001; 392 0.000000e+000 0.000000e+000 3.061225e+001; 393 0.000000e+000 4.000000e+000 3.061225e+001; 393 0.000000e+000 4.000000e+000 2.857143e+001; 393 0.000000e+000 0.000000e+000 2.857143e+001; 394 0.000000e+000 4.000000e+000 2.857143e+001; 394 0.000000e+000 4.000000e+000 2.653061e+001; 394 0.000000e+000 0.000000e+000 2.653061e+001; 395 0.000000e+000 4.000000e+000 2.653061e+001; 395 0.000000e+000 4.000000e+000 2.448980e+001; 395 0.000000e+000 0.000000e+000 2.448980e+001; 396 0.000000e+000 4.000000e+000 2.448980e+001; 396 0.000000e+000 4.000000e+000 2.244898e+001; 396 0.000000e+000 0.000000e+000 2.244898e+001; 397 0.000000e+000 4.000000e+000 2.244898e+001; 397 0.000000e+000 4.000000e+000 2.040816e+001; 397 0.000000e+000 0.000000e+000 2.244898e+001; 398 0.000000e+000 4.000000e+000 2.040816e+001; 398 0.000000e+000 4.000000e+000 1.836735e+001; 398 0.000000e+000 0.000000e+000 2.040816e+001; 399 0.000000e+000 4.000000e+000 1.836735e+001; 399 0.000000e+000 4.000000e+000 1.632653e+001; 399 0.000000e+000 0.000000e+000 1.632653e+001; 400 0.000000e+000 4.000000e+000 1.632653e+001; 400 0.000000e+000 4.000000e+000 1.428571e+001; 400 0.000000e+000 0.000000e+000 1.428571e+001; 401 0.000000e+000 4.000000e+000 1.428571e+001; 401 0.000000e+000 4.000000e+000 1.224490e+001; 401 0.000000e+000 0.000000e+000 1.224490e+001; 402 0.000000e+000 4.000000e+000 1.224490e+001; 402 0.000000e+000 4.000000e+000 1.020408e+001; 402 0.000000e+000 0.000000e+000 1.224490e+001; 403 0.000000e+000 4.000000e+000 1.020408e+001; 403 0.000000e+000 4.000000e+000 8.163265e+000; 403 0.000000e+000 0.000000e+000 8.163265e+000; 404 0.000000e+000 4.000000e+000 8.163265e+000; 404 0.000000e+000 4.000000e+000 6.122449e+000; 404 0.000000e+000 0.000000e+000 6.122449e+000; 405 0.000000e+000 4.000000e+000 6.122449e+000; 405 0.000000e+000 4.000000e+000 4.081633e+000; 405 0.000000e+000 0.000000e+000 6.122449e+000; 406 0.000000e+000 4.000000e+000 4.081633e+000; 406 0.000000e+000 4.000000e+000 2.040816e+000; 406 0.000000e+000 0.000000e+000 4.081633e+000; 407 0.000000e+000 4.000000e+000 2.040816e+000; 407 0.000000e+000 2.666667e+000 0.000000e+000; 407 0.000000e+000 0.000000e+000 2.040816e+000; 408 0.000000e+000 0.000000e+000 9.795918e+001; 408 0.000000e+000 2.666667e+000 1.000000e+002; 408 0.000000e+000 4.000000e+000 9.795918e+001; 409 4.000000e+000 0.000000e+000 0.000000e+000; 409 2.666667e+000 0.000000e+000 0.000000e+000; 409 4.000000e+000 1.333333e+000 0.000000e+000; 410 2.666667e+000 0.000000e+000 0.000000e+000; 410 1.333333e+000 0.000000e+000 0.000000e+000; 410 0.000000e+000 1.333333e+000 0.000000e+000; 411 1.333333e+000 0.000000e+000 0.000000e+000; 411 0.000000e+000 0.000000e+000 0.000000e+000; 411 0.000000e+000 1.333333e+000 0.000000e+000; 412 4.000000e+000 4.000000e+000 0.000000e+000; 412 4.000000e+000 2.666667e+000 0.000000e+000; 412 2.666667e+000 4.000000e+000 0.000000e+000; 413 4.000000e+000 2.666667e+000 0.000000e+000; 413 4.000000e+000 1.333333e+000 0.000000e+000; 413 2.666667e+000 0.000000e+000 0.000000e+000; 414 0.000000e+000 4.000000e+000 0.000000e+000; 414 1.333333e+000 4.000000e+000 0.000000e+000; 414 0.000000e+000 2.666667e+000 0.000000e+000; 415 1.333333e+000 4.000000e+000 0.000000e+000; 415 2.666667e+000 4.000000e+000 0.000000e+000; 415 4.000000e+000 2.666667e+000 0.000000e+000; 416 0.000000e+000 1.333333e+000 0.000000e+000; 416 0.000000e+000 2.666667e+000 0.000000e+000; 416 1.333333e+000 4.000000e+000 0.000000e+000; 417 2.666667e+000 0.000000e+000 0.000000e+000; 417 0.000000e+000 1.333333e+000 0.000000e+000; 417 1.333333e+000 4.000000e+000 0.000000e+000; 418 4.000000e+000 2.666667e+000 0.000000e+000; 418 2.666667e+000 0.000000e+000 0.000000e+000; 418 1.333333e+000 4.000000e+000 0.000000e+000; 419 0.000000e+000 0.000000e+000 1.000000e+002; 419 1.333333e+000 0.000000e+000 1.000000e+002; 419 0.000000e+000 1.333333e+000 1.000000e+002; 420 1.333333e+000 0.000000e+000 1.000000e+002; 420 2.666667e+000 0.000000e+000 1.000000e+002; 420 4.000000e+000 1.333333e+000 1.000000e+002; 421 2.666667e+000 0.000000e+000 1.000000e+002; 421 4.000000e+000 0.000000e+000 1.000000e+002; 421 4.000000e+000 1.333333e+000 1.000000e+002; 422 4.000000e+000 1.333333e+000 1.000000e+002; 422 4.000000e+000 2.666667e+000 1.000000e+002; 422 2.666667e+000 4.000000e+000 1.000000e+002; 423 4.000000e+000 2.666667e+000 1.000000e+002; 423 4.000000e+000 4.000000e+000 1.000000e+002; 423 2.666667e+000 4.000000e+000 1.000000e+002; 424 2.666667e+000 4.000000e+000 1.000000e+002; 424 1.333333e+000 4.000000e+000 1.000000e+002; 424 0.000000e+000 2.666667e+000 1.000000e+002; 425 1.333333e+000 4.000000e+000 1.000000e+002; 425 0.000000e+000 4.000000e+000 1.000000e+002; 425 0.000000e+000 2.666667e+000 1.000000e+002; 426 0.000000e+000 2.666667e+000 1.000000e+002; 426 0.000000e+000 1.333333e+000 1.000000e+002; 426 1.333333e+000 0.000000e+000 1.000000e+002; 427 1.333333e+000 0.000000e+000 1.000000e+002; 427 4.000000e+000 1.333333e+000 1.000000e+002; 427 2.666667e+000 4.000000e+000 1.000000e+002; 428 2.666667e+000 4.000000e+000 1.000000e+002; 428 0.000000e+000 2.666667e+000 1.000000e+002; 428 1.333333e+000 0.000000e+000 1.000000e+002]; G=E/(2*(1+mu)); C=1/(16*pi*G*(1-mu)); C1=3-4*mu; C2=1/(8*pi*(1-mu)); C3=1-2*mu; n=2; totalsize=size(xyzofelements); nelement=totalsize(1,1)/3; sizedispbcelements=size(dispbcelements); ndispbcelements=sizedispbcelements(1,1); sizefbcelements=size(fbcelements); nfbcelements=sizefbcelements(1,1); nx(1:nelement)=0; ny(1:nelement)=0; nz(1:nelement)=0; xbar(1:nelement)=0; ybar(1:nelement)=0; zbar(1:nelement)=0; J(1:nelement)=0; for i=1:nelement x1=xyzofelements(3*i-2,2); y1=xyzofelements(3*i-2,3); z1=xyzofelements(3*i-2,4); x2=xyzofelements(3*i-1,2); y2=xyzofelements(3*i-1,3); z2=xyzofelements(3*i-1,4); x3=xyzofelements(3*i,2); y3=xyzofelements(3*i,3); z3=xyzofelements(3*i,4); d=(((y2-y1)*(z3-z1)-(z2-z1)*(y3-y1))^2+((z2-z1)*(x3-x1)-(x2-x1)*(z3-z1))^2+((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1))^2)^0.5; nx(i)=((y2-y1)*(z3-z1)-(z2-z1)*(y3-y1))/d; ny(i)=((z2-z1)*(x3-x1)-(x2-x1)*(z3-z1))/d; nz(i)=((x2-x1)*(y3-y1)-(y2-y1)*(x3-x1))/d; alpha=sqrt((x1-x2)^2+(y1-y2)^2+(z1-z2)^2); beta=sqrt((x2-x3)^2+(y2-y3)^2+(z2-z3)^2); gamma=sqrt((x3-x1)^2+(y3-y1)^2+(z3-z1)^2); sigma=(alpha+beta+gamma)/2; J(i)=2*sqrt(sigma*(sigma-alpha)*(sigma-beta)*(sigma-gamma)); if (abs(nz(i))>=1/sqrt(3)) xbar(i)=(x2-x1)*0.25+(x3-x1)*0.5+x1; ybar(i)=(y2-y1)*0.25+(y3-y1)*0.5+y1; zbar(i)=-(nz(i))^(-1)*(nx(i)*(xbar(i)-x1)+ny(i)*(ybar(i)-y1))+z1; elseif (abs(nz(i))<1/sqrt(3) && abs(ny(i))>=1/sqrt(3)) xbar(i)=(x2-x1)*0.25+(x3-x1)*0.5+x1; zbar(i)=(z2-z1)*0.25+(z3-z1)*0.5+z1; ybar(i)=-(ny(i))^(-1)*(nx(i)*(xbar(i)-x1)+nz(i)*(zbar(i)-z1))+y1; elseif (abs(nz(i))<1/sqrt(3) && abs(ny(i))<1/sqrt(3)) ybar(i)=(y2-y1)*0.25+(y3-y1)*0.5+y1; zbar(i)=(z2-z1)*0.25+(z3-z1)*0.5+z1; xbar(i)=-(nx(i))^(-1)*(ny(i)*(ybar(i)-y1)+nz(i)*(zbar(i)-z1))+x1; end end t=[1/4+1/(4*sqrt(3)) 1/4+1/(4*sqrt(3)) 1/4-1/(4*sqrt(3)) 1/4-1/(4*sqrt(3)) 3/4+1/(4*sqrt(3)) 3/4+1/(4*sqrt(3)) 3/4-1/(4*sqrt(3)) 3/4-1/(4*sqrt(3)) 3/4+1/(4*sqrt(3)) 3/4+1/(4*sqrt(3)) 3/4-1/(4*sqrt(3)) 3/4-1/(4*sqrt(3)) 1/4+1/(4*sqrt(3)) 1/4+1/(4*sqrt(3)) 1/4-1/(4*sqrt(3)) 1/4-1/(4*sqrt(3))]; v=[1/4+1/(4*sqrt(3)) 1/4-1/(4*sqrt(3)) 1/4+1/(4*sqrt(3)) 1/4-1/(4*sqrt(3)) 1/4+1/(4*sqrt(3)) 1/4-1/(4*sqrt(3)) 1/4+1/(4*sqrt(3)) 1/4-1/(4*sqrt(3)) 3/4+1/(4*sqrt(3)) 3/4-1/(4*sqrt(3)) 3/4+1/(4*sqrt(3)) 3/4-1/(4*sqrt(3)) 3/4+1/(4*sqrt(3)) 3/4-1/(4*sqrt(3)) 3/4+1/(4*sqrt(3)) 3/4-1/(4*sqrt(3))]; u=t.*(1-v); K=zeros(3*i); F(1:3*i,1)=0; for i=1:nelement for j=1:nfbcelements if (abs(nz(fbcelements(j,1)))>=1/sqrt(3)) x1=xyzofelements(3*fbcelements(j,1)-2,2); y1=xyzofelements(3*fbcelements(j,1)-2,3); z1=xyzofelements(3*fbcelements(j,1)-2,4); x2=xyzofelements(3*fbcelements(j,1)-1,2); y2=xyzofelements(3*fbcelements(j,1)-1,3); z2=xyzofelements(3*fbcelements(j,1)-1,4); x3=xyzofelements(3*fbcelements(j,1),2); y3=xyzofelements(3*fbcelements(j,1),3); z3=xyzofelements(3*fbcelements(j,1),4); x=(x2-x1)*u+(x3-x1)*v+x1; y=(y2-y1)*u+(y3-y1)*v+y1; z=-(nz(fbcelements(j,1)))^(-1)*(nx(fbcelements(j,1))*(x-x1)+ny(fbcelements(j,1))*(y-y1))+z1; elseif (abs(nz(fbcelements(j,1)))<1/sqrt(3) && abs(ny(fbcelements(j,1)))>=1/sqrt(3)) x1=xyzofelements(3*fbcelements(j,1)-2,2); y1=xyzofelements(3*fbcelements(j,1)-2,3); z1=xyzofelements(3*fbcelements(j,1)-2,4); x2=xyzofelements(3*fbcelements(j,1)-1,2); y2=xyzofelements(3*fbcelements(j,1)-1,3); z2=xyzofelements(3*fbcelements(j,1)-1,4); x3=xyzofelements(3*fbcelements(j,1),2); y3=xyzofelements(3*fbcelements(j,1),3); z3=xyzofelements(3*fbcelements(j,1),4); x=(x2-x1)*u+(x3-x1)*v+x1; z=(z2-z1)*u+(z3-z1)*v+z1; y=-(ny(fbcelements(j,1)))^(-1)*(nx(fbcelements(j,1))*(x-x1)+nz(fbcelements(j,1))*(z-z1))+y1; elseif (abs(nz(fbcelements(j,1)))<1/sqrt(3) && abs(ny(fbcelements(j,1)))<1/sqrt(3)) x1=xyzofelements(3*fbcelements(j,1)-2,2); y1=xyzofelements(3*fbcelements(j,1)-2,3); z1=xyzofelements(3*fbcelements(j,1)-2,4); x2=xyzofelements(3*fbcelements(j,1)-1,2); y2=xyzofelements(3*fbcelements(j,1)-1,3); z2=xyzofelements(3*fbcelements(j,1)-1,4); x3=xyzofelements(3*fbcelements(j,1),2); y3=xyzofelements(3*fbcelements(j,1),3); z3=xyzofelements(3*fbcelements(j,1),4); y=(y2-y1)*u+(y3-y1)*v+y1; z=(z2-z1)*u+(z3-z1)*v+z1; x=-(nx(fbcelements(j,1)))^(-1)*(ny(fbcelements(j,1))*(y-y1)+nz(fbcelements(j,1))*(z-z1))+x1; end ITxxE=0; ITxyE=0; ITxzE=0; ITyxE=0; ITyyE=0; ITyzE=0; ITzxE=0; ITzyE=0; ITzzE=0; IUxxE=0; IUxyE=0; IUxzE=0; IUyxE=0; IUyyE=0; IUyzE=0; IUzxE=0; IUzyE=0; IUzzE=0; for k=1:16 r=sqrt((x(k)-xbar(i))^2+(y(k)-ybar(i))^2+(z(k)-zbar(i))^2); drbydx=(x(k)-xbar(i))/r; drbydy=(y(k)-ybar(i))/r; drbydz=(z(k)-zbar(i))/r; costheta=(1/r)*((x(k)-xbar(i))*nx(fbcelements(j,1))+(y(k)-ybar(i))*ny(fbcelements(j,1))+(z(k)-zbar(i))*nz(fbcelements(j,1))); Txx=(-C2/r^2)*((C3+3*drbydx^2)*costheta); Tyy=(-C2/r^2)*((C3+3*drbydy^2)*costheta); Tzz=(-C2/r^2)*((C3+3*drbydz^2)*costheta); Txy=(-C2/r^2)*(3*drbydx*drbydy*costheta-C3*(ny(fbcelements(j,1))*drbydx-nx(fbcelements(j,1))*drbydy)); Tyx=(-C2/r^2)*(3*drbydy*drbydx*costheta-C3*(nx(fbcelements(j,1))*drbydy-ny(fbcelements(j,1))*drbydx)); Tyz=(-C2/r^2)*(3*drbydy*drbydz*costheta-C3*(nz(fbcelements(j,1))*drbydy-ny(fbcelements(j,1))*drbydz)); Tzy=(-C2/r^2)*(3*drbydz*drbydy*costheta-C3*(ny(fbcelements(j,1))*drbydz-nz(fbcelements(j,1))*drbydy)); Txz=(-C2/r^2)*(3*drbydx*drbydz*costheta-C3*(nz(fbcelements(j,1))*drbydx-nx(fbcelements(j,1))*drbydz)); Tzx=(-C2/r^2)*(3*drbydz*drbydx*costheta-C3*(nx(fbcelements(j,1))*drbydz-nz(fbcelements(j,1))*drbydx)); ITxxE=ITxxE+Txx*(1-v(k))*J(fbcelements(j,1)); ITxyE=ITxyE+Txy*(1-v(k))*J(fbcelements(j,1)); ITxzE=ITxzE+Txz*(1-v(k))*J(fbcelements(j,1)); ITyxE=ITyxE+Tyx*(1-v(k))*J(fbcelements(j,1)); ITyyE=ITyyE+Tyy*(1-v(k))*J(fbcelements(j,1)); ITyzE=ITyzE+Tyz*(1-v(k))*J(fbcelements(j,1)); ITzxE=ITzxE+Tzx*(1-v(k))*J(fbcelements(j,1)); ITzyE=ITzyE+Tzy*(1-v(k))*J(fbcelements(j,1)); ITzzE=ITzzE+Tzz*(1-v(k))*J(fbcelements(j,1)); Uxx=(C/r)*(C1+drbydx^2); Uyy=(C/r)*(C1+drbydy^2); Uzz=(C/r)*(C1+drbydz^2); Uxy=(C/r)*drbydx*drbydy; Uyx=Uxy; Uyz=(C/r)*drbydy*drbydz; Uzy=Uyz; Uxz=(C/r)*drbydx*drbydz; Uzx=Uxz; IUxxE=IUxxE+Uxx*(1-v(k))*J(fbcelements(j,1)); IUxyE=IUxyE+Uxy*(1-v(k))*J(fbcelements(j,1)); IUxzE=IUxzE+Uxz*(1-v(k))*J(fbcelements(j,1)); IUyxE=IUyxE+Uyx*(1-v(k))*J(fbcelements(j,1)); IUyyE=IUyyE+Uyy*(1-v(k))*J(fbcelements(j,1)); IUyzE=IUyzE+Uyz*(1-v(k))*J(fbcelements(j,1)); IUzxE=IUzxE+Uzx*(1-v(k))*J(fbcelements(j,1)); IUzyE=IUzyE+Uzy*(1-v(k))*J(fbcelements(j,1)); IUzzE=IUzzE+Uzz*(1-v(k))*J(fbcelements(j,1)); end ITxxE=(1/16)*ITxxE; ITxyE=(1/16)*ITxyE; ITxzE=(1/16)*ITxzE; ITyxE=(1/16)*ITyxE; ITyyE=(1/16)*ITyyE; ITyzE=(1/16)*ITyzE; ITzxE=(1/16)*ITzxE; ITzyE=(1/16)*ITzyE; ITzzE=(1/16)*ITzzE; IUxxE=(1/16)*IUxxE; IUxyE=(1/16)*IUxyE; IUxzE=(1/16)*IUxzE; IUyxE=(1/16)*IUyxE; IUyyE=(1/16)*IUyyE; IUyzE=(1/16)*IUyzE; IUzxE=(1/16)*IUzxE; IUzyE=(1/16)*IUzyE; IUzzE=(1/16)*IUzzE; K(3*i-2,3*fbcelements(j,1)-2)=ITxxE; K(3*i-2,3*fbcelements(j,1)-1)=ITxyE; K(3*i-2,3*fbcelements(j,1))=ITxzE; K(3*i-1,3*fbcelements(j,1)-2)=ITyxE; K(3*i-1,3*fbcelements(j,1)-1)=ITyyE; K(3*i-1,3*fbcelements(j,1))=ITyzE; K(3*i,3*fbcelements(j,1)-2)=ITzxE; K(3*i,3*fbcelements(j,1)-1)=ITzyE; K(3*i,3*fbcelements(j,1))=ITzzE; F(3*i-2)=F(3*i-2)+fbcelements(j,2)*IUxxE+fbcelements(j,3)*IUxyE+fbcelements(j,4)*IUxzE; F(3*i-1)=F(3*i-1)+fbcelements(j,2)*IUyxE+fbcelements(j,3)*IUyyE+fbcelements(j,4)*IUyzE; F(3*i)=F(3*i)+fbcelements(j,2)*IUzxE+fbcelements(j,3)*IUzyE+fbcelements(j,4)*IUzzE; if i==fbcelements(j,1) K(3*fbcelements(j,1)-2,3*fbcelements(j,1)-2)=K(3*fbcelements(j,1)-2,3*fbcelements(j,1)-2)+0.5; K(3*fbcelements(j,1)-1,3*fbcelements(j,1)-1)=K(3*fbcelements(j,1)-1,3*fbcelements(j,1)-1)+0.5; K(3*fbcelements(j,1),3*fbcelements(j,1))=K(3*fbcelements(j,1),3*fbcelements(j,1))+0.5; end end for j=1:ndispbcelements if (abs(nz(dispbcelements(j,1)))>=1/sqrt(3)) x1=xyzofelements(3*dispbcelements(j,1)-2,2); y1=xyzofelements(3*dispbcelements(j,1)-2,3); z1=xyzofelements(3*dispbcelements(j,1)-2,4); x2=xyzofelements(3*dispbcelements(j,1)-1,2); y2=xyzofelements(3*dispbcelements(j,1)-1,3); z2=xyzofelements(3*dispbcelements(j,1)-1,4); x3=xyzofelements(3*dispbcelements(j,1),2); y3=xyzofelements(3*dispbcelements(j,1),3); z3=xyzofelements(3*dispbcelements(j,1),4); x=(x2-x1)*u+(x3-x1)*v+x1; y=(y2-y1)*u+(y3-y1)*v+y1; z=-(nz(dispbcelements(j,1)))^(-1)*(nx(dispbcelements(j,1))*(x-x1)+ny(dispbcelements(j,1))*(y-y1))+z1; elseif (abs(nz(dispbcelements(j,1)))<1/sqrt(3) && abs(ny(dispbcelements(j,1)))>=1/sqrt(3)) x1=xyzofelements(3*dispbcelements(j,1)-2,2); y1=xyzofelements(3*dispbcelements(j,1)-2,3); z1=xyzofelements(3*dispbcelements(j,1)-2,4); x2=xyzofelements(3*dispbcelements(j,1)-1,2); y2=xyzofelements(3*dispbcelements(j,1)-1,3); z2=xyzofelements(3*dispbcelements(j,1)-1,4); x3=xyzofelements(3*dispbcelements(j,1),2); y3=xyzofelements(3*dispbcelements(j,1),3); z3=xyzofelements(3*dispbcelements(j,1),4); x=(x2-x1)*u+(x3-x1)*v+x1; z=(z2-z1)*u+(z3-z1)*v+z1; y=-(ny(dispbcelements(j,1)))^(-1)*(nx(dispbcelements(j,1))*(x-x1)+nz(dispbcelements(j,1))*(z-z1))+y1; elseif (abs(nz(dispbcelements(j,1)))<1/sqrt(3) && abs(ny(dispbcelements(j,1)))<1/sqrt(3)) x1=xyzofelements(3*dispbcelements(j,1)-2,2); y1=xyzofelements(3*dispbcelements(j,1)-2,3); z1=xyzofelements(3*dispbcelements(j,1)-2,4); x2=xyzofelements(3*dispbcelements(j,1)-1,2); y2=xyzofelements(3*dispbcelements(j,1)-1,3); z2=xyzofelements(3*dispbcelements(j,1)-1,4); x3=xyzofelements(3*dispbcelements(j,1),2); y3=xyzofelements(3*dispbcelements(j,1),3); z3=xyzofelements(3*dispbcelements(j,1),4); y=(y2-y1)*u+(y3-y1)*v+y1; z=(z2-z1)*u+(z3-z1)*v+z1; x=-(nx(dispbcelements(j,1)))^(-1)*(ny(dispbcelements(j,1))*(y-y1)+nz(dispbcelements(j,1))*(z-z1))+x1; end ITxxE=0; ITxyE=0; ITxzE=0; ITyxE=0; ITyyE=0; ITyzE=0; ITzxE=0; ITzyE=0; ITzzE=0; IUxxE=0; IUxyE=0; IUxzE=0; IUyxE=0; IUyyE=0; IUyzE=0; IUzxE=0; IUzyE=0; IUzzE=0; for k=1:16 r=sqrt((x(k)-xbar(i))^2+(y(k)-ybar(i))^2+(z(k)-zbar(i))^2); drbydx=(x(k)-xbar(i))/r; drbydy=(y(k)-ybar(i))/r; drbydz=(z(k)-zbar(i))/r; costheta=(1/r)*((x(k)-xbar(i))*nx(dispbcelements(j,1))+(y(k)-ybar(i))*ny(dispbcelements(j,1))+(z(k)-zbar(i))*nz(dispbcelements(j,1))); Txx=(-C2/r^2)*((C3+3*drbydx^2)*costheta); Tyy=(-C2/r^2)*((C3+3*drbydy^2)*costheta); Tzz=(-C2/r^2)*((C3+3*drbydz^2)*costheta); Txy=(-C2/r^2)*(3*drbydx*drbydy*costheta-C3*(ny(dispbcelements(j,1))*drbydx-nx(dispbcelements(j,1))*drbydy)); Tyx=(-C2/r^2)*(3*drbydy*drbydx*costheta-C3*(nx(dispbcelements(j,1))*drbydy-ny(dispbcelements(j,1))*drbydx)); Tyz=(-C2/r^2)*(3*drbydy*drbydz*costheta-C3*(nz(dispbcelements(j,1))*drbydy-ny(dispbcelements(j,1))*drbydz)); Tzy=(-C2/r^2)*(3*drbydz*drbydy*costheta-C3*(ny(dispbcelements(j,1))*drbydz-nz(dispbcelements(j,1))*drbydy)); Txz=(-C2/r^2)*(3*drbydx*drbydz*costheta-C3*(nz(dispbcelements(j,1))*drbydx-nx(dispbcelements(j,1))*drbydz)); Tzx=(-C2/r^2)*(3*drbydz*drbydx*costheta-C3*(nx(dispbcelements(j,1))*drbydz-nz(dispbcelements(j,1))*drbydx)); ITxxE=ITxxE+Txx*(1-v(k))*J(dispbcelements(j,1)); ITxyE=ITxyE+Txy*(1-v(k))*J(dispbcelements(j,1)); ITxzE=ITxzE+Txz*(1-v(k))*J(dispbcelements(j,1)); ITyxE=ITyxE+Tyx*(1-v(k))*J(dispbcelements(j,1)); ITyyE=ITyyE+Tyy*(1-v(k))*J(dispbcelements(j,1)); ITyzE=ITyzE+Tyz*(1-v(k))*J(dispbcelements(j,1)); ITzxE=ITzxE+Tzx*(1-v(k))*J(dispbcelements(j,1)); ITzyE=ITzyE+Tzy*(1-v(k))*J(dispbcelements(j,1)); ITzzE=ITzzE+Tzz*(1-v(k))*J(dispbcelements(j,1)); Uxx=(C/r)*(C1+drbydx^2); Uyy=(C/r)*(C1+drbydy^2); Uzz=(C/r)*(C1+drbydz^2); Uxy=(C/r)*drbydx*drbydy; Uyx=Uxy; Uyz=(C/r)*drbydy*drbydz; Uzy=Uyz; Uxz=(C/r)*drbydx*drbydz; Uzx=Uxz; IUxxE=IUxxE+Uxx*(1-v(k))*J(dispbcelements(j,1)); IUxyE=IUxyE+Uxy*(1-v(k))*J(dispbcelements(j,1)); IUxzE=IUxzE+Uxz*(1-v(k))*J(dispbcelements(j,1)); IUyxE=IUyxE+Uyx*(1-v(k))*J(dispbcelements(j,1)); IUyyE=IUyyE+Uyy*(1-v(k))*J(dispbcelements(j,1)); IUyzE=IUyzE+Uyz*(1-v(k))*J(dispbcelements(j,1)); IUzxE=IUzxE+Uzx*(1-v(k))*J(dispbcelements(j,1)); IUzyE=IUzyE+Uzy*(1-v(k))*J(dispbcelements(j,1)); IUzzE=IUzzE+Uzz*(1-v(k))*J(dispbcelements(j,1)); end ITxxE=(1/16)*ITxxE; ITxyE=(1/16)*ITxyE; ITxzE=(1/16)*ITxzE; ITyxE=(1/16)*ITyxE; ITyyE=(1/16)*ITyyE; ITyzE=(1/16)*ITyzE; ITzxE=(1/16)*ITzxE; ITzyE=(1/16)*ITzyE; ITzzE=(1/16)*ITzzE; IUxxE=(1/16)*IUxxE; IUxyE=(1/16)*IUxyE; IUxzE=(1/16)*IUxzE; IUyxE=(1/16)*IUyxE; IUyyE=(1/16)*IUyyE; IUyzE=(1/16)*IUyzE; IUzxE=(1/16)*IUzxE; IUzyE=(1/16)*IUzyE; IUzzE=(1/16)*IUzzE; K(3*i-2,3*dispbcelements(j,1)-2)=IUxxE; K(3*i-2,3*dispbcelements(j,1)-1)=IUxyE; K(3*i-2,3*dispbcelements(j,1))=IUxzE; K(3*i-1,3*dispbcelements(j,1)-2)=IUyxE; K(3*i-1,3*dispbcelements(j,1)-1)=IUyyE; K(3*i-1,3*dispbcelements(j,1))=IUyzE; K(3*i,3*dispbcelements(j,1)-2)=IUzxE; K(3*i,3*dispbcelements(j,1)-1)=IUzyE; K(3*i,3*dispbcelements(j,1))=IUzzE; F(3*i-2)=F(3*i-2)+dispbcelements(j,2)*ITxxE+dispbcelements(j,3)*ITxyE+dispbcelements(j,4)*ITxzE; F(3*i-1)=F(3*i-1)+dispbcelements(j,2)*ITyxE+dispbcelements(j,3)*ITyyE+dispbcelements(j,4)*ITyzE; F(3*i)=F(3*i)+dispbcelements(j,2)*ITzxE+dispbcelements(j,3)*ITzyE+dispbcelements(j,4)*ITzzE; if i==dispbcelements(j,1) F(3*i-2)=F(3*i-2)+0.5*dispbcelements(j,2); F(3*i-1)=F(3*i-1)+0.5*dispbcelements(j,3); F(3*i)=F(3*i)+0.5*dispbcelements(j,4); end end end %Solve U=K\F; % Display unknown dofs for each element U