Skip to content
Snippets Groups Projects
Metavoxel GPU Validation 3.ipynb 395 KiB
Newer Older
Amira Abdel-Rahman's avatar
Amira Abdel-Rahman committed
2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 2017 2018 2019 2020 2021 2022 2023 2024 2025 2026 2027 2028 2029 2030 2031 2032 2033 2034 2035 2036 2037 2038 2039 2040 2041 2042 2043 2044 2045 2046 2047 2048 2049 2050 2051 2052 2053 2054 2055 2056 2057 2058 2059 2060 2061 2062 2063 2064 2065 2066 2067 2068 2069 2070 2071 2072 2073 2074 2075 2076 2077 2078 2079 2080 2081 2082 2083 2084 2085 2086 2087 2088 2089 2090 2091 2092 2093 2094 2095 2096 2097 2098 2099 2100 2101 2102 2103 2104 2105 2106 2107 2108 2109 2110 2111 2112 2113 2114 2115 2116 2117 2118 2119 2120 2121 2122 2123 2124 2125 2126 2127 2128 2129 2130 2131 2132 2133 2134 2135 2136 2137 2138 2139 2140 2141 2142 2143 2144 2145 2146 2147 2148 2149 2150 2151 2152 2153 2154 2155 2156 2157 2158 2159 2160 2161 2162 2163 2164 2165 2166 2167 2168 2169 2170 2171 2172 2173 2174 2175 2176 2177 2178 2179 2180 2181 2182 2183 2184 2185 2186 2187 2188 2189 2190 2191 2192 2193 2194 2195 2196 2197 2198 2199 2200 2201 2202 2203 2204 2205 2206 2207 2208 2209 2210 2211 2212 2213 2214 2215 2216 2217 2218 2219 2220 2221 2222 2223 2224 2225 2226 2227 2228 2229 2230 2231 2232 2233 2234 2235 2236 2237 2238 2239 2240 2241 2242 2243 2244 2245 2246 2247 2248 2249 2250 2251 2252 2253 2254 2255 2256 2257 2258 2259 2260 2261 2262 2263 2264 2265 2266 2267 2268 2269 2270 2271 2272 2273 2274 2275 2276 2277 2278 2279 2280 2281 2282 2283 2284 2285 2286 2287 2288 2289 2290 2291 2292 2293 2294 2295 2296 2297 2298 2299 2300 2301 2302 2303 2304 2305 2306 2307 2308 2309 2310 2311 2312 2313 2314 2315 2316 2317 2318 2319 2320 2321 2322 2323 2324 2325 2326 2327 2328 2329 2330 2331 2332 2333 2334 2335 2336 2337 2338 2339 2340 2341 2342 2343 2344 2345 2346 2347 2348 2349 2350 2351 2352 2353 2354 2355 2356 2357 2358 2359 2360 2361 2362 2363 2364 2365 2366 2367 2368 2369 2370 2371 2372 2373 2374 2375 2376 2377 2378 2379 2380 2381 2382 2383 2384 2385 2386 2387 2388 2389 2390 2391 2392 2393 2394 2395 2396 2397 2398 2399 2400 2401 2402 2403 2404 2405 2406 2407 2408 2409 2410 2411 2412 2413 2414 2415 2416 2417 2418 2419 2420 2421 2422 2423 2424 2425 2426 2427 2428 2429 2430 2431 2432 2433 2434 2435 2436 2437 2438 2439 2440 2441 2442 2443 2444 2445 2446 2447 2448 2449 2450 2451 2452 2453 2454 2455 2456 2457 2458 2459 2460 2461 2462 2463 2464 2465 2466 2467 2468 2469 2470 2471 2472 2473 2474 2475 2476 2477 2478 2479 2480 2481 2482 2483 2484 2485 2486 2487 2488 2489 2490 2491 2492 2493 2494 2495 2496 2497 2498 2499 2500 2501 2502 2503 2504 2505 2506 2507 2508 2509 2510 2511 2512 2513 2514 2515 2516 2517 2518 2519 2520 2521 2522 2523 2524 2525 2526 2527 2528 2529 2530 2531 2532 2533 2534 2535 2536 2537 2538 2539 2540 2541 2542 2543 2544 2545 2546 2547 2548 2549 2550 2551 2552 2553 2554 2555 2556 2557 2558 2559 2560 2561 2562 2563 2564 2565 2566 2567 2568 2569 2570 2571 2572 2573 2574 2575 2576 2577 2578 2579 2580 2581 2582 2583 2584 2585 2586 2587 2588 2589 2590 2591 2592 2593 2594 2595 2596 2597 2598 2599 2600 2601 2602 2603 2604 2605 2606 2607 2608 2609 2610 2611 2612 2613 2614 2615 2616 2617 2618 2619 2620 2621 2622 2623 2624 2625 2626 2627 2628 2629 2630 2631 2632 2633 2634 2635 2636 2637 2638 2639 2640 2641 2642 2643 2644 2645 2646 2647 2648 2649 2650 2651 2652 2653 2654 2655 2656 2657 2658 2659 2660 2661 2662 2663 2664 2665 2666 2667 2668 2669 2670 2671 2672 2673 2674 2675 2676 2677 2678 2679 2680 2681 2682 2683 2684 2685 2686 2687 2688 2689 2690 2691 2692 2693 2694 2695 2696 2697 2698 2699 2700 2701 2702 2703 2704 2705 2706 2707 2708 2709 2710 2711 2712 2713 2714 2715 2716 2717 2718 2719 2720 2721 2722 2723 2724 2725 2726 2727 2728 2729 2730 2731 2732 2733 2734 2735 2736 2737 2738 2739 2740 2741 2742 2743 2744 2745 2746 2747 2748 2749 2750 2751 2752 2753 2754 2755 2756 2757 2758 2759 2760 2761 2762 2763 2764 2765 2766 2767 2768 2769 2770 2771 2772 2773 2774 2775 2776 2777 2778 2779 2780 2781 2782 2783 2784 2785 2786 2787 2788 2789 2790 2791 2792 2793 2794 2795 2796 2797 2798 2799 2800 2801 2802 2803 2804 2805 2806 2807 2808 2809 2810 2811 2812 2813 2814 2815 2816 2817 2818 2819 2820 2821 2822 2823 2824 2825 2826 2827 2828 2829 2830 2831 2832 2833 2834 2835 2836 2837 2838 2839 2840 2841 2842 2843 2844 2845 2846 2847 2848 2849 2850 2851 2852 2853 2854 2855 2856 2857 2858 2859 2860 2861 2862 2863 2864 2865 2866 2867 2868 2869 2870 2871 2872 2873 2874 2875 2876 2877 2878 2879 2880 2881 2882 2883 2884 2885 2886 2887 2888 2889 2890 2891 2892 2893 2894 2895 2896 2897 2898 2899 2900 2901 2902 2903 2904 2905 2906 2907 2908 2909 2910 2911 2912 2913 2914 2915 2916 2917 2918 2919 2920 2921 2922 2923 2924 2925 2926 2927 2928 2929 2930 2931 2932 2933 2934 2935 2936 2937 2938 2939 2940 2941 2942 2943 2944 2945 2946 2947 2948 2949 2950 2951 2952 2953 2954 2955 2956 2957 2958 2959 2960 2961 2962 2963 2964 2965 2966 2967 2968 2969 2970 2971 2972 2973 2974 2975 2976 2977 2978 2979 2980 2981 2982 2983 2984 2985 2986 2987 2988 2989 2990 2991 2992 2993 2994 2995 2996 2997 2998 2999 3000
    "            A = edge[\"area\"]\n",
    "            G=1.0#todo shear_modulus\n",
    "            ixx = 1.0#todo section ixx\n",
    "            iyy = 1.0#todo section.iyy#\n",
    "            l0=length\n",
    "            j=1.0;#todo check\n",
    "            l02 = l0 * l0\n",
    "            l03 = l0 * l0 * l0\n",
    "            \n",
    "            ################################\n",
    "            mass=1\n",
    "            nu=0.35\n",
    "            W = 75\n",
    "            L = W/sqrt(2)\n",
    "            L=length\n",
    "            n_min = 1\n",
    "            n_max = 7\n",
    "            # Cross Section inputs, must be floats\n",
    "            E = 2000e9  # MPa\n",
    "            G = E * 1 / 3  # MPa\n",
    "            h = 2.38/100.0  # mm\n",
    "            b = 2.38/100.0  # mm\n",
    "            rho = 7.85e-9 / 3  # kg/mm^3\n",
    "            S = h * b\n",
    "            Sy = (S * (6 + 12 * nu + 6 * nu^2)/ (7 + 12 * nu + 4 * nu^2))\n",
    "            # For solid rectangular cross section (width=b, depth=d & ( b < d )):\n",
    "            Q = 1 / 3 - 0.2244 / (min(h / b, b / h) + 0.1607)\n",
    "            J = Q * min(h * b^3, b * h^3)\n",
    "            \n",
    "            \n",
    "            ##if voxels\n",
    "            #nu=0\n",
    "            #L=l\n",
    "            #a1 = E*L # EA/L : Units of N/m\n",
    "            #a2 = E * L*L*L / (12.0*(1+nu)) # GJ/L : Units of N-m\n",
    "            #b1 = E*L # 12EI/L^3 : Units of N/m\n",
    "            #b2 = E*L*L/2.0 # 6EI/L^2 : Units of N (or N-m/m: torque related to linear distance)\n",
    "            #b3 = E*L*L*L/6.0 # 2EI/L : Units of N-m\n",
    "            \n",
    "            I= b*h^3/12\n",
    "#             J=b*h*(b*b+h*h)/12\n",
    "            a1=E*b*h/L\n",
    "            a2=G*J/L\n",
    "            b1=12*E*I/(L^3)\n",
    "            b2=6*E*I/(L^2)\n",
    "            b3=2*E*I/(L)\n",
    "            \n",
    "            # Cm = rho * area * length /6.0\n",
    "            # Ck= E * area / length \n",
    "\n",
    "            # m = [[2 1];[1 2]]\n",
    "            # k = [[1 -1];[-1 1]]\n",
    "            \n",
    "\n",
    "            k = [[E*A/l0  0  0  0  0  0  -E*A/l0  0  0  0  0  0];[0  12*E*ixx/l03  0  0  0  6*E*ixx/l02  0  -12*E*ixx/l03  0  0  0  6*E*ixx/l02];[0  0  12*E*iyy/l03  0  -6*E*iyy/l02  0  0  0  -12*E*iyy/l03  0  -6*E*iyy/l02  0];[0  0  0  G*j/l0  0  0  0  0  0  -G*j/l0  0  0];[0  0  -6*E*iyy/l02  0  4*E*iyy/l0  0  0  0  6*E*iyy/l02  0  2*E*iyy/l0  0];[0  6*E*ixx/l02  0  0  0  4*E*ixx/l0  0  -6*E*ixx/l02  0  0  0  2*E*ixx/l0];[-E*A/l0  0  0  0  0  0  E*A/l0  0  0  0  0  0];[0  -12*E*ixx/l03  0  0  0  -6*E*ixx/l02  0  12*E*ixx/l03  0  0  0  -6*E*ixx/l02];[0  0  -12*E*iyy/l03  0  6*E*iyy/l02  0  0  0  12*E*iyy/l03  0  6*E*iyy/l02  0];[0  0  0  -G*j/l0  0  0  0  0  0  G*j/l0  0  0];[0  0  -6*E*iyy/l02  0  2*E*iyy/l0  0  0  0  6*E*iyy/l02  0  4*E*iyy/l0  0];[0  6*E*ixx/l02  0  0  0  2*E*ixx/l0  0  -6*E*ixx/l02  0  0  0  4*E*ixx/l0]]\n",
    "            k= [[ a1  0  0  0   0   0  -a1  0   0   0   0   0  ];\n",
    "                [  0 b1  0  0   0   b2  0  -b1  0   0   0   b2 ];\n",
    "                [  0  0  b1 0  -b2  0   0   0  -b1  0  -b2  0  ];\n",
    "                [  0  0  0  a2  0   0   0   0   0  -a2  0   0  ];\n",
    "                [  0  0  0  0  2b3  0   0   0   b2  0   b3  0  ];\n",
    "                [  0  0  0  0   0  2b3  0  -b2  0   0   0   b3 ];\n",
    "                [  0  0  0  0   0   0   a1  0   0   0   0   0  ];\n",
    "                [  0  0  0  0   0   0   0   b1  0   0   0  -b2 ];\n",
    "                [  0  0  0  0   0   0   0   0   b1  0   b2  0  ];\n",
    "                [  0  0  0  0   0   0   0   0   0   a2  0   0  ];\n",
    "                [  0  0  0  0   0   0   0   0   0    0 2b3  0  ];\n",
    "                [  0  0  0  0   0   0   0   0   0    0  0  2b3 ]]\n",
    "            # find rotated mass and stifness matrices\n",
    "            tau = rotation_matrix(element_vector, x_axis,y_axis,z_axis)\n",
    "\n",
    "            # m_r=transpose(tau)*m*tau\n",
    "            k_r=transpose(tau)*k*tau\n",
    "\n",
    "            # change from element to global coordinate\n",
    "            index= dofs.+1\n",
    "\n",
    "            B=zeros((12,ndofs))\n",
    "            for i in 1:12\n",
    "                  B[i,index[i]]=1.0\n",
    "            end\n",
    "\n",
    "\n",
    "            # M_rG= transpose(B)*m_r*B\n",
    "            K_rG= transpose(B)*k_r*B\n",
    "\n",
    "            # M += Cm .* M_rG\n",
    "            # K += Ck .* K_rG\n",
    "            K +=  K_rG\n",
    "\n",
    "        end\n",
    "        \n",
    "        \n",
    "        # construct the force vector\n",
    "        F=zeros(ndofs)\n",
    "        remove_indices=[];\n",
    "        for node in nodes\n",
    "            #insert!(F,i, value);\n",
    "            #F=vcat(F,value)\n",
    "            \n",
    "            \n",
    "            i=parse(Int,node[\"id\"][2:end])\n",
    "            f=node[\"force\"]\n",
    "            \n",
    "            # println(f)\n",
    "            F[(i)*6+1]=f[\"x\"]/1.0\n",
    "            F[(i)*6+2]=f[\"y\"]/1.0\n",
    "            F[(i)*6+3]=f[\"z\"]/1.0\n",
    "            F[(i)*6+4]=0\n",
    "            F[(i)*6+5]=0\n",
    "            F[(i)*6+6]=0\n",
    "            Load+=f[\"y\"]/1.0\n",
    "            if (F[(i)*6+2]!=0)\n",
    "                append!(topNodesIndices,i)\n",
    "            end\n",
    "            \n",
    "            dofs = convert(Array{Int}, node[\"degrees_of_freedom\"]).+1\n",
    "            restrained_dofs=node[\"restrained_degrees_of_freedom\"]\n",
    "            for (index, value) in enumerate(dofs)\n",
    "                if restrained_dofs[index]\n",
    "                    append!( remove_indices, value)\n",
    "                end\n",
    "            end\n",
    "            \n",
    "        end\n",
    "\n",
    "        #println(remove_indices)\n",
    "        #print(K)\n",
    "        #print(F)\n",
    "        \n",
    "\n",
    "        #M = M[setdiff(1:end, remove_indices), :]\n",
    "        K = K[setdiff(1:end, remove_indices), :]\n",
    "\n",
    "        #M = M[:,setdiff(1:end, remove_indices)]\n",
    "        K = K[:,setdiff(1:end, remove_indices)]\n",
    "\n",
    "        F = F[setdiff(1:end, remove_indices)]\n",
    "        U=zeros(ndofs)\n",
    "        \n",
    "        return M,K,F,U,remove_indices\n",
    "    end\n",
    "\n",
    "    \n",
    "    function updateDisplacement(setup, X)\n",
    "        nodes= setup[\"nodes\"]\n",
    "        i=0\n",
    "        for node in nodes\n",
    "            \n",
    "#             if !node[\"restrained_degrees_of_freedom\"][2]\n",
    "                #i=parse(Int,node[\"id\"][2:end])\n",
    "                node[\"displacement\"][\"x\"]=X[(i)*6+1]/15\n",
    "                node[\"displacement\"][\"y\"]=X[(i)*6+2]/15\n",
    "                node[\"displacement\"][\"z\"]=X[(i)*6+3]/15\n",
    "                node[\"angle\"][\"x\"]=X[(i)*6+4]\n",
    "                node[\"angle\"][\"y\"]=X[(i)*6+5]\n",
    "                node[\"angle\"][\"z\"]=X[(i)*6+6]\n",
    "                append!(displacementFEA,[Vector3(X[(i)*6+1],X[(i)*6+2],X[(i)*6+3])])\n",
    "                i=i+1\n",
    "#             else\n",
    "#                 append!(displacementFEA,[Vector3(0,0,0)])\n",
    "#             end\n",
    "        end\n",
    "    end\n",
    "    \n",
    "    #######################################################\n",
    "\n",
    "    function get_stresses(setup)\n",
    "        nodes      = setup[\"nodes\"]\n",
    "        edges      = setup[\"edges\"]\n",
    "        ndofs      = length(nodes)*6\n",
    "\n",
    "        x_axis     = [1 0 0]\n",
    "        y_axis     = [0 1 0]\n",
    "        z_axis     = [0 0 1]\n",
    "\n",
    "        # find the stresses in each member\n",
    "        stresses=zeros(length(edges))\n",
    "        max11=-10e6\n",
    "        min11=10e6\n",
    "        for edge in edges\n",
    "            #degrees_of_freedom = properties[\"degrees_of_freedom\"]\n",
    "\n",
    "            element=parse(Int,edge[\"id\"][2:end])\n",
    "\n",
    "            # find the nodes that the lements connects\n",
    "            fromNode = nodes[edge[\"source\"]+1]\n",
    "            toNode = nodes[edge[\"target\"]+1]\n",
    "\n",
    "            # the coordinates for each node\n",
    "            fromPoint = [fromNode[\"position\"][\"x\"]/100.0 fromNode[\"position\"][\"y\"]/100.0 fromNode[\"position\"][\"z\"]/100.0]\n",
    "            toPoint = [toNode[\"position\"][\"x\"]/100.0 toNode[\"position\"][\"y\"]/100.0 toNode[\"position\"][\"z\"]/100.0]\n",
    "\n",
    "            # find the degrees of freedom for each node\n",
    "            dofs = convert(Array{Int}, fromNode[\"degrees_of_freedom\"])\n",
    "            dofs=vcat(dofs,convert(Array{Int}, toNode[\"degrees_of_freedom\"]))\n",
    "\n",
    "            element_vector=toPoint-fromPoint\n",
    "\n",
    "\n",
    "            # find rotated mass and stifness matrices\n",
    "            tau = rotation_matrix(element_vector, x_axis,y_axis,z_axis)\n",
    "\n",
    "            # i1=parse(Int,fromNode[\"id\"][2:end])\n",
    "            # i2=parse(Int,toNode[\"id\"][2:end])\n",
    "\n",
    "            # global_displacements=[X[(i1)*6+1] X[(i1)*6+2] X[(i1)*6+3] X[(i1)*6+4] X[(i1)*6+5] X[(i1)*6+6] X[(i2)*6+1] X[(i2)*6+2] X[(i2)*6+3] X[(i2)*6+4] X[(i2)*6+5] X[(i2)*6+6]] # todo change\n",
    "            global_displacements=[fromNode[\"displacement\"][\"x\"] fromNode[\"displacement\"][\"y\"] fromNode[\"displacement\"][\"z\"] fromNode[\"angle\"][\"x\"] fromNode[\"angle\"][\"y\"] fromNode[\"angle\"][\"z\"] toNode[\"displacement\"][\"x\"] toNode[\"displacement\"][\"y\"] toNode[\"displacement\"][\"z\"] toNode[\"angle\"][\"x\"] toNode[\"angle\"][\"y\"] toNode[\"angle\"][\"z\"]] # todo change\n",
    "\n",
    "            # nodal displacement\n",
    "\n",
    "            q=tau*transpose(global_displacements)\n",
    "            # println(q)\n",
    "            # calculate the strain and stresses\n",
    "            strain =(q[7]-q[1])/norm(element_vector)\n",
    "            E = edge[\"stiffness\"]# youngs modulus\n",
    "            stress=E.*strain\n",
    "            edge[\"stress\"]=stress\n",
    "            if stress>max11\n",
    "                max11=stress\n",
    "            end\n",
    "            if stress<min11\n",
    "                min11=stress\n",
    "            end\n",
    "            # println(element)\n",
    "            # println(stress)\n",
    "        end\n",
    "\n",
    "\n",
    "\n",
    "        setup[\"viz\"][\"minStress\"]=min11\n",
    "        setup[\"viz\"][\"maxStress\"]=max11\n",
    "        return stresses\n",
    "    end\n",
    "    \n",
    "    function initialize(setup)\n",
    "        nodes      = setup[\"nodes\"]\n",
    "        ndofs      = length(nodes)*6\n",
    "        \n",
    "        i=0\n",
    "        for node in nodes\n",
    "            dg=[]\n",
    "            for ii in 0:5\n",
    "                append!(dg,i+ii) \n",
    "            end\n",
    "            i+=6\n",
    "            node[\"degrees_of_freedom\"]=dg\n",
    "        end\n",
    "    end\n",
    "\n",
    "    #######################################################\n",
    "    function solveFea(setup)\n",
    "        // # determine the global matrices\n",
    "        initialize(setup)\n",
    "        M,K,F,U,ind=get_matrices(setup)\n",
    "        \n",
    "        #println(M)\n",
    "        #println(K)\n",
    "        #println(F)\n",
    "\n",
    "        #evals=eigvals(K,M)\n",
    "        #evecs=eigvecs(K,M)\n",
    "        #frequencies=sqrt.(evals)\n",
    "        X=inv(K)*F\n",
    "        U[setdiff(1:end, ind)]=X\n",
    "\n",
    "        updateDisplacement(setup, U)\n",
    "\n",
    "        # determine the stresses in each element\n",
    "        stresses=get_stresses(setup)\n",
    "    end\n",
    "    #######################################################\n",
    "    displacementFEA=[]\n",
    "    Load=0\n",
    "    topNodesIndices=[]\n",
    "    solveFea(setup)\n",
    "    return displacementFEA,Load,topNodesIndices\n",
    "end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "recommendedTimeStep (generic function with 1 method)"
      ]
     },
     "execution_count": 17,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "function recommendedTimeStep()\n",
    "    \n",
    "    \n",
    "#     #find the largest natural frequency (sqrt(k/m)) that anything in the simulation will experience, then multiply by 2*pi and invert to get the optimally largest timestep that should retain stability\n",
    "#     MaxFreq2 = 0.0; #maximum frequency in the simulation in rad/sec\n",
    "    \n",
    "#     for e in 1:length(edges) #for each edge\n",
    "#         CVX_Link* pL = (*it);\n",
    "#         #axial\n",
    "#         m1 = mass(),  m2 = mass()\n",
    "#         thisMaxFreq2 = axialStiffness()/(m1<m2?m1:m2)\n",
    "#         if (thisMaxFreq2 > MaxFreq2) \n",
    "#             MaxFreq2 = thisMaxFreq2;\n",
    "#         end\n",
    "\n",
    "#         #rotational will always be less than or equal\n",
    "#     end\n",
    "        \n",
    "\n",
    "#     if (MaxFreq2 <= 0.0) #didn't find anything (i.e no links) check for individual voxelss\n",
    "#         for n in 1:length(nodes) #for each node\n",
    "#             thisMaxFreq2 = youngsModulus*nomSize/mass;\n",
    "#             #thisMaxFreq2 = (*it)->mat->youngsModulus()*(*it)->mat->nomSize/(*it)->mat->mass();\n",
    "#             if (thisMaxFreq2 > MaxFreq2) \n",
    "#                     MaxFreq2 = thisMaxFreq2;\n",
    "#             end\n",
    "#         end\n",
    "#     end\n",
    "\n",
    "#     if (MaxFreq2 <= 0.0) \n",
    "#         return 0.0\n",
    "#     else \n",
    "#         return 1.0/(6.283185*sqrt(MaxFreq2)) #the optimal timestep is to advance one radian of the highest natural frequency\n",
    "#     end\n",
    "    MaxFreq2=E*size/mass\n",
    "    return 1/(6.283185*sqrt(MaxFreq2))\n",
    "end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "getSetup (generic function with 1 method)"
      ]
     },
     "execution_count": 20,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "function getSetup(latticeSize)\n",
    "    setup = Dict()\n",
    "    name=string(\"../json/setupTestUni$latticeSize\",\".json\")\n",
    "#     open(\"../json/setupValid2.json\", \"r\") do f\n",
    "#     open(\"../json/setupTest.json\", \"r\") do f\n",
    "    # open(\"../json/trialJulia.json\", \"r\") do f\n",
    "#     open(\"../json/setupTestUni4.json\", \"r\") do f\n",
    "    # open(\"../json/setupChiral.json\", \"r\") do f\n",
    "#     open(\"../json/setupTestCubeUni10.json\", \"r\") do f\n",
    "    open(name, \"r\") do f\n",
    "#         global setup\n",
    "        dicttxt = String(read(f))  # file information to string\n",
    "        setup=JSON.parse(dicttxt)  # parse and transform data\n",
    "    end\n",
    "\n",
    "    setup=setup[\"setup\"]\n",
    "    return setup\n",
    "end"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 21,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "4-element Array{Any,1}:\n",
       " 4\n",
       " 7\n",
       " 8\n",
       " 9"
      ]
     },
     "execution_count": 21,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "latticeSize=1\n",
    "setup=getSetup(latticeSize)\n",
    "displacementFEA,Load,topNodesIndices=fea(setup)\n",
    "topNodesIndices"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 22,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "5-element Array{Float64,1}:\n",
       " 0.0\n",
       " 0.0\n",
       " 0.0\n",
       " 0.0\n",
       " 0.0"
      ]
     },
     "execution_count": 22,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "DDisplacements=[[],[],[],[],[]]\n",
    "DDisplacementsFEA=[[],[],[],[],[]]\n",
    "EsFEA=[0.0,0,0,0,0]\n",
    "Es=[0.0,0,0,0,0]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 46,
   "metadata": {},
   "outputs": [],
   "source": [
    "latticeSize=5\n",
    "setup=getSetup(latticeSize)\n",
    "displacementFEA,Load,topNodesIndices=fea(setup);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 48,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "dt: 7.294855212986817e-7\n",
      "first timestep took 0.0006283 seconds\n",
      "ran latticeSize 5 with 540 voxels and 1800 edges for 2000 time steps took 0.8313075 seconds\n",
      "EsFEA:[4.516432560172212, 7.389712166701619, 9.998185539814536, 12.509610413282042, 14.980612218071958]\n",
      "Es:[-1.7207658256817755, -1.7187653726043761, -1.8909709475438676, -2.125645742684135, -2.322015017271702]\n",
      "FEA displacement= -8.696319301969073e-10,converged displacement= 3.0929663544073403e-9\n"
     ]
    },
    {
     "data": {
      "image/svg+xml": [
       "<?xml version=\"1.0\" encoding=\"utf-8\"?>\n",
       "<svg xmlns=\"http://www.w3.org/2000/svg\" xmlns:xlink=\"http://www.w3.org/1999/xlink\" width=\"600\" height=\"400\" viewBox=\"0 0 2400 1600\">\n",
       "<defs>\n",
       "  <clipPath id=\"clip0800\">\n",
       "    <rect x=\"0\" y=\"0\" width=\"2400\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<path clip-path=\"url(#clip0800)\" d=\"\n",
       "M0 1600 L2400 1600 L2400 0 L0 0  Z\n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip0801\">\n",
       "    <rect x=\"480\" y=\"0\" width=\"1681\" height=\"1600\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<path clip-path=\"url(#clip0800)\" d=\"\n",
       "M407.037 1425.62 L2352.76 1425.62 L2352.76 121.675 L407.037 121.675  Z\n",
       "  \" fill=\"#ffffff\" fill-rule=\"evenodd\" fill-opacity=\"1\"/>\n",
       "<defs>\n",
       "  <clipPath id=\"clip0802\">\n",
       "    <rect x=\"407\" y=\"121\" width=\"1947\" height=\"1305\"/>\n",
       "  </clipPath>\n",
       "</defs>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  461.187,1425.62 461.187,121.675 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  920.083,1425.62 920.083,121.675 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1378.98,1425.62 1378.98,121.675 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  1837.87,1425.62 1837.87,121.675 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  2296.77,1425.62 2296.77,121.675 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  407.037,1415.51 2352.76,1415.51 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  407.037,1209.98 2352.76,1209.98 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  407.037,1004.45 2352.76,1004.45 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  407.037,798.913 2352.76,798.913 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  407.037,593.381 2352.76,593.381 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  407.037,387.849 2352.76,387.849 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#000000; stroke-width:2; stroke-opacity:0.1; fill:none\" points=\"\n",
       "  407.037,182.317 2352.76,182.317 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,1425.62 2352.76,1425.62 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,1425.62 407.037,121.675 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  461.187,1425.62 461.187,1409.97 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  920.083,1425.62 920.083,1409.97 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1378.98,1425.62 1378.98,1409.97 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  1837.87,1425.62 1837.87,1409.97 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  2296.77,1425.62 2296.77,1409.97 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,1415.51 430.386,1415.51 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,1209.98 430.386,1209.98 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,1004.45 430.386,1004.45 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,798.913 430.386,798.913 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,593.381 430.386,593.381 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,387.849 430.386,387.849 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0800)\" style=\"stroke:#000000; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  407.037,182.317 430.386,182.317 \n",
       "  \"/>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 461.187, 1479.62)\" x=\"461.187\" y=\"1479.62\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 920.083, 1479.62)\" x=\"920.083\" y=\"1479.62\">500</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1378.98, 1479.62)\" x=\"1378.98\" y=\"1479.62\">1000</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 1837.87, 1479.62)\" x=\"1837.87\" y=\"1479.62\">1500</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:middle;\" transform=\"rotate(0, 2296.77, 1479.62)\" x=\"2296.77\" y=\"1479.62\">2000</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 201.946, 1439.24)\" x=\"201.946\" y=\"1439.24\">-</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 230.056, 1439.24)\" x=\"230.056\" y=\"1439.24\">1×10</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 338.454, 1411.83)\" x=\"338.454\" y=\"1411.83\">-</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 361.293, 1411.83)\" x=\"361.293\" y=\"1411.83\">7</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:end;\" transform=\"rotate(0, 383.037, 1227.48)\" x=\"383.037\" y=\"1227.48\">0</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 230.056, 1028.17)\" x=\"230.056\" y=\"1028.17\">1×10</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 338.454, 1000.76)\" x=\"338.454\" y=\"1000.76\">-</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 361.293, 1000.76)\" x=\"361.293\" y=\"1000.76\">7</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 230.056, 822.641)\" x=\"230.056\" y=\"822.641\">2×10</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 338.454, 795.231)\" x=\"338.454\" y=\"795.231\">-</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 361.293, 795.231)\" x=\"361.293\" y=\"795.231\">7</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 230.056, 617.109)\" x=\"230.056\" y=\"617.109\">3×10</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 338.454, 589.699)\" x=\"338.454\" y=\"589.699\">-</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 361.293, 589.699)\" x=\"361.293\" y=\"589.699\">7</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 230.056, 411.577)\" x=\"230.056\" y=\"411.577\">4×10</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 338.454, 384.167)\" x=\"338.454\" y=\"384.167\">-</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 361.293, 384.167)\" x=\"361.293\" y=\"384.167\">7</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:48px; text-anchor:start;\" transform=\"rotate(0, 230.056, 206.045)\" x=\"230.056\" y=\"206.045\">5×10</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 338.454, 178.634)\" x=\"338.454\" y=\"178.634\">-</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:38px; text-anchor:start;\" transform=\"rotate(0, 361.293, 178.634)\" x=\"361.293\" y=\"178.634\">7</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:84px; text-anchor:middle;\" transform=\"rotate(0, 1379.9, 73.2)\" x=\"1379.9\" y=\"73.2\">5 Voxel Convergence Study</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:66px; text-anchor:middle;\" transform=\"rotate(0, 1379.9, 1559.48)\" x=\"1379.9\" y=\"1559.48\">timestep</text>\n",
       "</g>\n",
       "<g clip-path=\"url(#clip0800)\">\n",
       "<text style=\"fill:#000000; fill-opacity:1; font-family:Arial,Helvetica Neue,Helvetica,sans-serif; font-size:66px; text-anchor:middle;\" transform=\"rotate(-90, 89.2861, 773.647)\" x=\"89.2861\" y=\"773.647\">displacement</text>\n",
       "</g>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#009af9; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  462.105,1209.98 463.023,1206.46 463.941,1200.48 464.858,1192.66 465.776,1183.43 466.694,1173.09 467.612,1161.87 468.529,1149.96 469.447,1137.5 470.365,1124.62 \n",
       "  471.283,1111.41 472.201,1097.95 473.118,1084.3 474.036,1070.52 474.954,1056.65 475.872,1042.73 476.79,1028.78 477.707,1014.84 478.625,1000.92 479.543,987.061 \n",
       "  480.461,973.266 481.379,959.559 482.296,945.958 483.214,932.481 484.132,919.147 485.05,905.972 485.968,892.976 486.885,880.175 487.803,867.589 488.721,855.233 \n",
       "  489.639,843.125 490.556,831.282 491.474,819.717 492.392,808.444 493.31,797.476 494.228,786.823 495.145,776.493 496.063,766.492 496.981,756.824 497.899,747.491 \n",
       "  498.817,738.491 499.734,729.822 500.652,721.477 501.57,713.447 502.488,705.721 503.406,698.287 504.323,691.128 505.241,684.226 506.159,677.562 507.077,671.115 \n",
       "  507.995,664.861 508.912,658.777 509.83,652.838 510.748,647.016 511.666,641.286 512.583,635.62 513.501,629.992 514.419,624.375 515.337,618.743 516.255,613.07 \n",
       "  517.172,607.331 518.09,601.502 519.008,595.562 519.926,589.488 520.844,583.263 521.761,576.868 522.679,570.287 523.597,563.508 524.515,556.518 525.433,549.308 \n",
       "  526.35,541.871 527.268,534.201 528.186,526.297 529.104,518.157 530.022,509.783 530.939,501.179 531.857,492.351 532.775,483.307 533.693,474.057 534.61,464.612 \n",
       "  535.528,454.987 536.446,445.196 537.364,435.256 538.282,425.187 539.199,415.006 540.117,404.735 541.035,394.396 541.953,384.01 542.871,373.602 543.788,363.196 \n",
       "  544.706,352.814 545.624,342.483 546.542,332.226 547.46,322.069 548.377,312.036 549.295,302.152 550.213,292.441 551.131,282.925 552.049,273.629 552.966,264.575 \n",
       "  553.884,255.784 554.802,247.276 555.72,239.073 556.637,231.191 557.555,223.65 558.473,216.465 559.391,209.652 560.309,203.225 561.226,197.196 562.144,191.578 \n",
       "  563.062,186.381 563.98,181.612 564.898,177.281 565.815,173.393 566.733,169.954 567.651,166.966 568.569,164.433 569.487,162.355 570.404,160.733 571.322,159.564 \n",
       "  572.24,158.848 573.158,158.579 574.076,158.754 574.993,159.367 575.911,160.412 576.829,161.88 577.747,163.764 578.664,166.055 579.582,168.742 580.5,171.816 \n",
       "  581.418,175.265 582.336,179.078 583.253,183.241 584.171,187.744 585.089,192.573 586.007,197.714 586.925,203.154 587.842,208.879 588.76,214.876 589.678,221.131 \n",
       "  590.596,227.628 591.514,234.356 592.431,241.298 593.349,248.442 594.267,255.774 595.185,263.279 596.103,270.946 597.02,278.759 597.938,286.707 598.856,294.777 \n",
       "  599.774,302.956 600.691,311.233 601.609,319.595 602.527,328.032 603.445,336.533 604.363,345.086 605.28,353.682 606.198,362.311 607.116,370.965 608.034,379.633 \n",
       "  608.952,388.307 609.869,396.981 610.787,405.645 611.705,414.293 612.623,422.917 613.541,431.513 614.458,440.072 615.376,448.591 616.294,457.064 617.212,465.485 \n",
       "  618.13,473.85 619.047,482.155 619.965,490.397 620.883,498.571 621.801,506.674 622.718,514.704 623.636,522.656 624.554,530.53 625.472,538.323 626.39,546.032 \n",
       "  627.307,553.655 628.225,561.192 629.143,568.64 630.061,575.998 630.979,583.264 631.896,590.438 632.814,597.519 633.732,604.505 634.65,611.396 635.568,618.191 \n",
       "  636.485,624.889 637.403,631.489 638.321,637.991 639.239,644.395 640.157,650.699 641.074,656.903 641.992,663.006 642.91,669.008 643.828,674.908 644.745,680.705 \n",
       "  645.663,686.399 646.581,691.989 647.499,697.475 648.417,702.855 649.334,708.129 650.252,713.295 651.17,718.354 652.088,723.304 653.006,728.144 653.923,732.874 \n",
       "  654.841,737.491 655.759,741.996 656.677,746.387 657.595,750.662 658.512,754.822 659.43,758.863 660.348,762.786 661.266,766.589 662.184,770.271 663.101,773.83 \n",
       "  664.019,777.266 664.937,780.576 665.855,783.76 666.772,786.816 667.69,789.743 668.608,792.54 669.526,795.205 670.444,797.737 671.361,800.135 672.279,802.398 \n",
       "  673.197,804.525 674.115,806.514 675.033,808.364 675.95,810.075 676.868,811.645 677.786,813.073 678.704,814.359 679.622,815.502 680.539,816.5 681.457,817.354 \n",
       "  682.375,818.063 683.293,818.626 684.211,819.043 685.128,819.313 686.046,819.436 686.964,819.412 687.882,819.242 688.799,818.924 689.717,818.459 690.635,817.848 \n",
       "  691.553,817.09 692.471,816.186 693.388,815.137 694.306,813.943 695.224,812.605 696.142,811.124 697.06,809.501 697.977,807.736 698.895,805.832 699.813,803.788 \n",
       "  700.731,801.607 701.649,799.291 702.566,796.84 703.484,794.257 704.402,791.543 705.32,788.7 706.238,785.731 707.155,782.636 708.073,779.42 708.991,776.083 \n",
       "  709.909,772.629 710.826,769.059 711.744,765.377 712.662,761.586 713.58,757.687 714.498,753.684 715.415,749.58 716.333,745.377 717.251,741.08 718.169,736.691 \n",
       "  719.087,732.213 720.004,727.65 720.922,723.005 721.84,718.282 722.758,713.484 723.676,708.614 724.593,703.677 725.511,698.675 726.429,693.614 727.347,688.496 \n",
       "  728.265,683.325 729.182,678.105 730.1,672.84 731.018,667.533 731.936,662.19 732.853,656.813 733.771,651.406 734.689,645.974 735.607,640.52 736.525,635.049 \n",
       "  737.442,629.564 738.36,624.069 739.278,618.568 740.196,613.065 741.114,607.564 742.031,602.069 742.949,596.584 743.867,591.113 744.785,585.658 745.703,580.225 \n",
       "  746.62,574.817 747.538,569.437 748.456,564.089 749.374,558.777 750.292,553.505 751.209,548.275 752.127,543.092 753.045,537.958 753.963,532.877 754.88,527.853 \n",
       "  755.798,522.888 756.716,517.986 757.634,513.149 758.552,508.381 759.469,503.685 760.387,499.063 761.305,494.518 762.223,490.054 763.141,485.671 764.058,481.374 \n",
       "  764.976,477.164 765.894,473.044 766.812,469.016 767.73,465.082 768.647,461.245 769.565,457.506 770.483,453.867 771.401,450.33 772.319,446.897 773.236,443.57 \n",
       "  774.154,440.35 775.072,437.238 775.99,434.235 776.907,431.344 777.825,428.566 778.743,425.9 779.661,423.349 780.579,420.913 781.496,418.594 782.414,416.391 \n",
       "  783.332,414.305 784.25,412.338 785.168,410.489 786.085,408.759 787.003,407.148 787.921,405.656 788.839,404.284 789.757,403.031 790.674,401.898 791.592,400.883 \n",
       "  792.51,399.988 793.428,399.212 794.346,398.554 795.263,398.014 796.181,397.591 797.099,397.286 798.017,397.096 798.934,397.022 799.852,397.062 800.77,397.216 \n",
       "  801.688,397.483 802.606,397.862 803.523,398.351 804.441,398.95 805.359,399.657 806.277,400.471 807.195,401.391 808.112,402.415 809.03,403.543 809.948,404.771 \n",
       "  810.866,406.1 811.784,407.527 812.701,409.051 813.619,410.67 814.537,412.383 815.455,414.187 816.373,416.081 817.29,418.062 818.208,420.13 819.126,422.282 \n",
       "  820.044,424.516 820.962,426.831 821.879,429.223 822.797,431.692 823.715,434.235 824.633,436.85 825.55,439.534 826.468,442.286 827.386,445.104 828.304,447.985 \n",
       "  829.222,450.927 830.139,453.928 831.057,456.985 831.975,460.096 832.893,463.26 833.811,466.473 834.728,469.733 835.646,473.038 836.564,476.386 837.482,479.774 \n",
       "  838.4,483.2 839.317,486.662 840.235,490.156 841.153,493.682 842.071,497.235 842.989,500.815 843.906,504.418 844.824,508.043 845.742,511.686 846.66,515.346 \n",
       "  847.577,519.02 848.495,522.706 849.413,526.401 850.331,530.103 851.249,533.81 852.166,537.519 853.084,541.229 854.002,544.936 854.92,548.638 855.838,552.334 \n",
       "  856.755,556.021 857.673,559.696 858.591,563.358 859.509,567.005 860.427,570.633 861.344,574.241 862.262,577.828 863.18,581.39 864.098,584.925 865.016,588.433 \n",
       "  865.933,591.909 866.851,595.354 867.769,598.764 868.687,602.138 869.604,605.473 870.522,608.769 871.44,612.022 872.358,615.232 873.276,618.397 874.193,621.514 \n",
       "  875.111,624.582 876.029,627.6 876.947,630.566 877.865,633.477 878.782,636.334 879.7,639.134 880.618,641.875 881.536,644.557 882.454,647.178 883.371,649.736 \n",
       "  884.289,652.231 885.207,654.661 886.125,657.025 887.043,659.322 887.96,661.551 888.878,663.711 889.796,665.8 890.714,667.818 891.631,669.764 892.549,671.637 \n",
       "  893.467,673.436 894.385,675.161 895.303,676.81 896.22,678.384 897.138,679.881 898.056,681.302 898.974,682.645 899.892,683.91 900.809,685.097 901.727,686.205 \n",
       "  902.645,687.235 903.563,688.185 904.481,689.056 905.398,689.848 906.316,690.561 907.234,691.193 908.152,691.747 909.07,692.221 909.987,692.616 910.905,692.932 \n",
       "  911.823,693.169 912.741,693.328 913.658,693.408 914.576,693.411 915.494,693.336 916.412,693.184 917.33,692.956 918.247,692.653 919.165,692.274 920.083,691.821 \n",
       "  921.001,691.293 921.919,690.693 922.836,690.021 923.754,689.277 924.672,688.463 925.59,687.579 926.508,686.627 927.425,685.607 928.343,684.521 929.261,683.369 \n",
       "  930.179,682.152 931.097,680.873 932.014,679.532 932.932,678.129 933.85,676.668 934.768,675.148 935.685,673.571 936.603,671.939 937.521,670.252 938.439,668.513 \n",
       "  939.357,666.722 940.274,664.882 941.192,662.993 942.11,661.057 943.028,659.076 943.946,657.052 944.863,654.984 945.781,652.877 946.699,650.73 947.617,648.546 \n",
       "  948.535,646.326 949.452,644.072 950.37,641.785 951.288,639.468 952.206,637.122 953.124,634.748 954.041,632.348 954.959,629.925 955.877,627.479 956.795,625.013 \n",
       "  957.712,622.527 958.63,620.025 959.548,617.507 960.466,614.976 961.384,612.432 962.301,609.878 963.219,607.316 964.137,604.747 965.055,602.173 965.973,599.595 \n",
       "  966.89,597.015 967.808,594.436 968.726,591.858 969.644,589.283 970.562,586.713 971.479,584.15 972.397,581.595 973.315,579.05 974.233,576.515 975.151,573.994 \n",
       "  976.068,571.487 976.986,568.997 977.904,566.523 978.822,564.069 979.739,561.635 980.657,559.223 981.575,556.834 982.493,554.47 983.411,552.132 984.328,549.821 \n",
       "  985.246,547.539 986.164,545.287 987.082,543.067 988,540.879 988.917,538.724 989.835,536.605 990.753,534.522 991.671,532.475 992.589,530.468 993.506,528.499 \n",
       "  994.424,526.571 995.342,524.684 996.26,522.839 997.178,521.038 998.095,519.281 999.013,517.569 999.931,515.902 1000.85,514.283 1001.77,512.71 1002.68,511.186 \n",
       "  1003.6,509.71 1004.52,508.284 1005.44,506.908 1006.36,505.583 1007.27,504.308 1008.19,503.086 1009.11,501.915 1010.03,500.797 1010.94,499.732 1011.86,498.72 \n",
       "  1012.78,497.762 1013.7,496.858 1014.62,496.008 1015.53,495.213 1016.45,494.472 1017.37,493.786 1018.29,493.154 1019.2,492.578 1020.12,492.056 1021.04,491.59 \n",
       "  1021.96,491.178 1022.88,490.821 1023.79,490.519 1024.71,490.272 1025.63,490.079 1026.55,489.94 1027.46,489.855 1028.38,489.824 1029.3,489.846 1030.22,489.922 \n",
       "  1031.14,490.05 1032.05,490.23 1032.97,490.463 1033.89,490.747 1034.81,491.081 1035.72,491.467 1036.64,491.902 1037.56,492.386 1038.48,492.919 1039.4,493.5 \n",
       "  1040.31,494.129 1041.23,494.804 1042.15,495.525 1043.07,496.292 1043.98,497.103 1044.9,497.957 1045.82,498.855 1046.74,499.795 1047.66,500.776 1048.57,501.797 \n",
       "  1049.49,502.858 1050.41,503.957 1051.33,505.094 1052.25,506.268 1053.16,507.478 1054.08,508.722 1055,510 1055.92,511.31 1056.83,512.652 1057.75,514.025 \n",
       "  1058.67,515.428 1059.59,516.859 1060.51,518.317 1061.42,519.802 1062.34,521.311 1063.26,522.845 1064.18,524.402 1065.09,525.981 1066.01,527.581 1066.93,529.2 \n",
       "  1067.85,530.837 1068.77,532.492 1069.68,534.162 1070.6,535.848 1071.52,537.547 1072.44,539.259 1073.35,540.983 1074.27,542.716 1075.19,544.459 1076.11,546.21 \n",
       "  1077.03,547.967 1077.94,549.73 1078.86,551.497 1079.78,553.268 1080.7,555.041 1081.61,556.814 1082.53,558.588 1083.45,560.36 1084.37,562.13 1085.29,563.896 \n",
       "  1086.2,565.657 1087.12,567.413 1088.04,569.162 1088.96,570.902 1089.87,572.634 1090.79,574.356 1091.71,576.067 1092.63,577.765 1093.55,579.45 1094.46,581.121 \n",
       "  1095.38,582.777 1096.3,584.417 1097.22,586.04 1098.13,587.645 1099.05,589.23 1099.97,590.796 1100.89,592.342 1101.81,593.865 1102.72,595.366 1103.64,596.844 \n",
       "  1104.56,598.298 1105.48,599.727 1106.39,601.13 1107.31,602.507 1108.23,603.856 1109.15,605.178 1110.07,606.471 1110.98,607.736 1111.9,608.97 1112.82,610.174 \n",
       "  1113.74,611.347 1114.65,612.488 1115.57,613.597 1116.49,614.673 1117.41,615.717 1118.33,616.726 1119.24,617.702 1120.16,618.643 1121.08,619.549 1122,620.42 \n",
       "  1122.91,621.255 1123.83,622.055 1124.75,622.818 1125.67,623.544 1126.59,624.234 1127.5,624.887 1128.42,625.502 1129.34,626.08 1130.26,626.621 1131.18,627.124 \n",
       "  1132.09,627.589 1133.01,628.017 1133.93,628.406 1134.85,628.758 1135.76,629.071 1136.68,629.347 1137.6,629.585 1138.52,629.786 1139.44,629.948 1140.35,630.073 \n",
       "  1141.27,630.161 1142.19,630.211 1143.11,630.225 1144.02,630.201 1144.94,630.141 1145.86,630.045 1146.78,629.912 1147.7,629.744 1148.61,629.541 1149.53,629.302 \n",
       "  1150.45,629.028 1151.37,628.72 1152.28,628.379 1153.2,628.003 1154.12,627.595 1155.04,627.154 1155.96,626.68 1156.87,626.175 1157.79,625.639 1158.71,625.073 \n",
       "  1159.63,624.476 1160.54,623.85 1161.46,623.194 1162.38,622.511 1163.3,621.8 1164.22,621.061 1165.13,620.296 1166.05,619.506 1166.97,618.69 1167.89,617.849 \n",
       "  1168.8,616.985 1169.72,616.098 1170.64,615.188 1171.56,614.257 1172.48,613.304 1173.39,612.332 1174.31,611.34 1175.23,610.329 1176.15,609.3 1177.06,608.254 \n",
       "  1177.98,607.192 1178.9,606.113 1179.82,605.021 1180.74,603.914 1181.65,602.793 1182.57,601.66 1183.49,600.516 1184.41,599.361 1185.32,598.195 1186.24,597.021 \n",
       "  1187.16,595.837 1188.08,594.647 1189,593.449 1189.91,592.245 1190.83,591.037 1191.75,589.823 1192.67,588.606 1193.58,587.387 1194.5,586.165 1195.42,584.942 \n",
       "  1196.34,583.719 1197.26,582.496 1198.17,581.274 1199.09,580.054 1200.01,578.837 1200.93,577.623 1201.85,576.414 1202.76,575.209 1203.68,574.01 1204.6,572.818 \n",
       "  1205.52,571.632 1206.43,570.455 1207.35,569.286 1208.27,568.127 1209.19,566.977 1210.11,565.839 1211.02,564.711 1211.94,563.595 1212.86,562.493 1213.78,561.403 \n",
       "  1214.69,560.327 1215.61,559.266 1216.53,558.22 1217.45,557.189 1218.37,556.175 1219.28,555.177 1220.2,554.197 1221.12,553.235 1222.04,552.291 1222.95,551.366 \n",
       "  1223.87,550.46 1224.79,549.574 1225.71,548.708 1226.63,547.863 1227.54,547.039 1228.46,546.236 1229.38,545.456 1230.3,544.698 1231.21,543.962 1232.13,543.249 \n",
       "  1233.05,542.56 1233.97,541.894 1234.89,541.252 1235.8,540.634 1236.72,540.04 1237.64,539.471 1238.56,538.927 1239.47,538.408 1240.39,537.914 1241.31,537.446 \n",
       "  1242.23,537.003 1243.15,536.586 1244.06,536.194 1244.98,535.829 1245.9,535.489 1246.82,535.176 1247.73,534.888 1248.65,534.627 1249.57,534.391 1250.49,534.182 \n",
       "  1251.41,533.999 1252.32,533.842 1253.24,533.711 1254.16,533.606 1255.08,533.527 1255.99,533.473 1256.91,533.445 1257.83,533.443 1258.75,533.466 1259.67,533.514 \n",
       "  1260.58,533.587 1261.5,533.685 1262.42,533.807 1263.34,533.954 1264.25,534.125 1265.17,534.32 1266.09,534.538 1267.01,534.78 1267.93,535.045 1268.84,535.332 \n",
       "  1269.76,535.642 1270.68,535.974 1271.6,536.327 1272.52,536.702 1273.43,537.098 1274.35,537.515 1275.27,537.952 1276.19,538.409 1277.1,538.885 1278.02,539.38 \n",
       "  1278.94,539.893 1279.86,540.425 1280.78,540.975 1281.69,541.541 1282.61,542.125 1283.53,542.725 1284.45,543.34 1285.36,543.971 1286.28,544.617 1287.2,545.277 \n",
       "  1288.12,545.951 1289.04,546.638 1289.95,547.338 1290.87,548.05 1291.79,548.774 1292.71,549.509 1293.62,550.254 1294.54,551.01 1295.46,551.775 1296.38,552.55 \n",
       "  1297.3,553.333 1298.21,554.123 1299.13,554.921 1300.05,555.726 1300.97,556.538 1301.88,557.354 1302.8,558.176 1303.72,559.003 1304.64,559.834 1305.56,560.668 \n",
       "  1306.47,561.504 1307.39,562.344 1308.31,563.185 1309.23,564.027 1310.14,564.87 1311.06,565.713 1311.98,566.556 1312.9,567.398 1313.82,568.239 1314.73,569.077 \n",
       "  1315.65,569.913 1316.57,570.746 1317.49,571.575 1318.4,572.4 1319.32,573.221 1320.24,574.037 1321.16,574.847 1322.08,575.651 1322.99,576.448 1323.91,577.239 \n",
       "  1324.83,578.022 1325.75,578.797 1326.66,579.564 1327.58,580.322 1328.5,581.071 1329.42,581.81 1330.34,582.539 1331.25,583.257 1332.17,583.965 1333.09,584.661 \n",
       "  1334.01,585.345 1334.92,586.018 1335.84,586.678 1336.76,587.326 1337.68,587.96 1338.6,588.581 1339.51,589.188 1340.43,589.781 1341.35,590.36 1342.27,590.924 \n",
       "  1343.18,591.474 1344.1,592.008 1345.02,592.527 1345.94,593.03 1346.86,593.517 1347.77,593.988 1348.69,594.443 1349.61,594.881 1350.53,595.303 1351.45,595.708 \n",
       "  1352.36,596.096 1353.28,596.466 1354.2,596.82 1355.12,597.156 1356.03,597.474 1356.95,597.775 1357.87,598.058 1358.79,598.323 1359.71,598.57 1360.62,598.8 \n",
       "  1361.54,599.011 1362.46,599.205 1363.38,599.38 1364.29,599.538 1365.21,599.677 1366.13,599.799 1367.05,599.903 1367.97,599.988 1368.88,600.056 1369.8,600.106 \n",
       "  1370.72,600.138 1371.64,600.153 1372.55,600.15 1373.47,600.13 1374.39,600.092 1375.31,600.037 1376.23,599.966 1377.14,599.877 1378.06,599.771 1378.98,599.649 \n",
       "  1379.9,599.511 1380.81,599.356 1381.73,599.185 1382.65,598.998 1383.57,598.796 1384.49,598.579 1385.4,598.346 1386.32,598.098 1387.24,597.836 1388.16,597.559 \n",
       "  1389.07,597.268 1389.99,596.964 1390.91,596.645 1391.83,596.314 1392.75,595.969 1393.66,595.612 1394.58,595.242 1395.5,594.86 1396.42,594.466 1397.33,594.061 \n",
       "  1398.25,593.645 1399.17,593.217 1400.09,592.78 1401.01,592.332 1401.92,591.875 1402.84,591.408 1403.76,590.932 1404.68,590.447 1405.59,589.954 1406.51,589.453 \n",
       "  1407.43,588.944 1408.35,588.428 1409.27,587.906 1410.18,587.376 1411.1,586.841 1412.02,586.3 1412.94,585.753 1413.85,585.202 1414.77,584.646 1415.69,584.086 \n",
       "  1416.61,583.522 1417.53,582.954 1418.44,582.384 1419.36,581.811 1420.28,581.235 1421.2,580.658 1422.12,580.079 1423.03,579.499 1423.95,578.918 1424.87,578.337 \n",
       "  1425.79,577.756 1426.7,577.175 1427.62,576.595 1428.54,576.016 1429.46,575.438 1430.38,574.862 1431.29,574.289 1432.21,573.718 1433.13,573.15 1434.05,572.585 \n",
       "  1434.96,572.024 1435.88,571.466 1436.8,570.913 1437.72,570.365 1438.64,569.821 1439.55,569.283 1440.47,568.75 1441.39,568.223 1442.31,567.702 1443.22,567.188 \n",
       "  1444.14,566.68 1445.06,566.18 1445.98,565.686 1446.9,565.201 1447.81,564.723 1448.73,564.253 1449.65,563.792 1450.57,563.339 1451.48,562.895 1452.4,562.46 \n",
       "  1453.32,562.034 1454.24,561.618 1455.16,561.212 1456.07,560.816 1456.99,560.429 1457.91,560.053 1458.83,559.688 1459.74,559.333 1460.66,558.989 1461.58,558.656 \n",
       "  1462.5,558.334 1463.42,558.024 1464.33,557.725 1465.25,557.437 1466.17,557.161 1467.09,556.897 1468,556.644 1468.92,556.404 1469.84,556.176 1470.76,555.959 \n",
       "  1471.68,555.755 1472.59,555.563 1473.51,555.384 1474.43,555.216 1475.35,555.061 1476.26,554.919 1477.18,554.789 1478.1,554.671 1479.02,554.566 1479.94,554.473 \n",
       "  1480.85,554.392 1481.77,554.324 1482.69,554.268 1483.61,554.225 1484.52,554.193 1485.44,554.174 1486.36,554.167 1487.28,554.172 1488.2,554.189 1489.11,554.219 \n",
       "  1490.03,554.259 1490.95,554.312 1491.87,554.376 1492.79,554.452 1493.7,554.539 1494.62,554.638 1495.54,554.747 1496.46,554.868 1497.37,554.999 1498.29,555.141 \n",
       "  1499.21,555.294 1500.13,555.457 1501.05,555.63 1501.96,555.814 1502.88,556.007 1503.8,556.21 1504.72,556.423 1505.63,556.645 1506.55,556.876 1507.47,557.115 \n",
       "  1508.39,557.364 1509.31,557.621 1510.22,557.887 1511.14,558.16 1512.06,558.441 1512.98,558.73 1513.89,559.027 1514.81,559.33 1515.73,559.641 1516.65,559.958 \n",
       "  1517.57,560.281 1518.48,560.611 1519.4,560.947 1520.32,561.288 1521.24,561.635 1522.15,561.987 1523.07,562.343 1523.99,562.705 1524.91,563.071 1525.83,563.441 \n",
       "  1526.74,563.815 1527.66,564.193 1528.58,564.574 1529.5,564.958 1530.41,565.345 1531.33,565.734 1532.25,566.126 1533.17,566.52 1534.09,566.915 1535,567.312 \n",
       "  1535.92,567.71 1536.84,568.11 1537.76,568.51 1538.67,568.91 1539.59,569.311 1540.51,569.711 1541.43,570.112 1542.35,570.511 1543.26,570.91 1544.18,571.308 \n",
       "  1545.1,571.705 1546.02,572.1 1546.93,572.493 1547.85,572.884 1548.77,573.272 1549.69,573.659 1550.61,574.042 1551.52,574.423 1552.44,574.8 1553.36,575.174 \n",
       "  1554.28,575.544 1555.19,575.91 1556.11,576.272 1557.03,576.63 1557.95,576.983 1558.87,577.332 1559.78,577.676 1560.7,578.015 1561.62,578.348 1562.54,578.676 \n",
       "  1563.46,578.998 1564.37,579.315 1565.29,579.625 1566.21,579.93 1567.13,580.228 1568.04,580.52 1568.96,580.805 1569.88,581.083 1570.8,581.354 1571.72,581.619 \n",
       "  1572.63,581.876 1573.55,582.126 1574.47,582.369 1575.39,582.604 1576.3,582.832 1577.22,583.051 1578.14,583.263 1579.06,583.468 1579.98,583.664 1580.89,583.852 \n",
       "  1581.81,584.032 1582.73,584.204 1583.65,584.367 1584.56,584.523 1585.48,584.67 1586.4,584.808 1587.32,584.938 1588.24,585.06 1589.15,585.173 1590.07,585.277 \n",
       "  1590.99,585.374 1591.91,585.461 1592.82,585.54 1593.74,585.61 1594.66,585.672 1595.58,585.726 1596.5,585.77 1597.41,585.807 1598.33,585.834 1599.25,585.854 \n",
       "  1600.17,585.865 1601.08,585.867 1602,585.862 1602.92,585.848 1603.84,585.826 1604.76,585.795 1605.67,585.757 1606.59,585.711 1607.51,585.656 1608.43,585.594 \n",
       "  1609.34,585.524 1610.26,585.447 1611.18,585.362 1612.1,585.269 1613.02,585.169 1613.93,585.062 1614.85,584.948 1615.77,584.827 1616.69,584.698 1617.6,584.563 \n",
       "  1618.52,584.422 1619.44,584.274 1620.36,584.119 1621.28,583.958 1622.19,583.791 1623.11,583.618 1624.03,583.44 1624.95,583.255 1625.86,583.065 1626.78,582.87 \n",
       "  1627.7,582.669 1628.62,582.464 1629.54,582.253 1630.45,582.038 1631.37,581.819 1632.29,581.594 1633.21,581.366 1634.12,581.134 1635.04,580.897 1635.96,580.657 \n",
       "  1636.88,580.414 1637.8,580.167 1638.71,579.917 1639.63,579.664 1640.55,579.408 1641.47,579.15 1642.39,578.889 1643.3,578.626 1644.22,578.36 1645.14,578.093 \n",
       "  1646.06,577.824 1646.97,577.554 1647.89,577.282 1648.81,577.009 1649.73,576.735 1650.65,576.461 1651.56,576.185 1652.48,575.909 1653.4,575.633 1654.32,575.357 \n",
       "  1655.23,575.081 1656.15,574.805 1657.07,574.53 1657.99,574.255 1658.91,573.981 1659.82,573.708 1660.74,573.436 1661.66,573.165 1662.58,572.896 1663.49,572.629 \n",
       "  1664.41,572.363 1665.33,572.099 1666.25,571.838 1667.17,571.578 1668.08,571.321 1669,571.067 1669.92,570.815 1670.84,570.566 1671.75,570.32 1672.67,570.077 \n",
       "  1673.59,569.838 1674.51,569.602 1675.43,569.369 1676.34,569.14 1677.26,568.915 1678.18,568.694 1679.1,568.477 1680.01,568.264 1680.93,568.055 1681.85,567.851 \n",
       "  1682.77,567.651 1683.69,567.456 1684.6,567.265 1685.52,567.079 1686.44,566.898 1687.36,566.722 1688.27,566.551 1689.19,566.385 1690.11,566.225 1691.03,566.069 \n",
       "  1691.95,565.919 1692.86,565.774 1693.78,565.635 1694.7,565.501 1695.62,565.373 1696.53,565.25 1697.45,565.133 1698.37,565.022 1699.29,564.917 1700.21,564.817 \n",
       "  1701.12,564.723 1702.04,564.635 1702.96,564.552 1703.88,564.476 1704.79,564.405 1705.71,564.341 1706.63,564.282 1707.55,564.229 1708.47,564.182 1709.38,564.141 \n",
       "  1710.3,564.106 1711.22,564.077 1712.14,564.053 1713.06,564.035 1713.97,564.024 1714.89,564.017 1715.81,564.017 1716.73,564.023 1717.64,564.034 1718.56,564.05 \n",
       "  1719.48,564.073 1720.4,564.101 1721.32,564.134 1722.23,564.173 1723.15,564.217 1724.07,564.267 1724.99,564.321 1725.9,564.381 1726.82,564.447 1727.74,564.517 \n",
       "  1728.66,564.592 1729.58,564.672 1730.49,564.757 1731.41,564.846 1732.33,564.941 1733.25,565.04 1734.16,565.143 1735.08,565.251 1736,565.362 1736.92,565.479 \n",
       "  1737.84,565.599 1738.75,565.723 1739.67,565.851 1740.59,565.983 1741.51,566.119 1742.42,566.258 1743.34,566.401 1744.26,566.546 1745.18,566.696 1746.1,566.848 \n",
       "  1747.01,567.003 1747.93,567.161 1748.85,567.322 1749.77,567.486 1750.68,567.652 1751.6,567.82 1752.52,567.991 1753.44,568.164 1754.36,568.339 1755.27,568.516 \n",
       "  1756.19,568.694 1757.11,568.875 1758.03,569.056 1758.94,569.24 1759.86,569.424 1760.78,569.61 1761.7,569.796 1762.62,569.984 1763.53,570.172 1764.45,570.361 \n",
       "  1765.37,570.551 1766.29,570.741 1767.2,570.931 1768.12,571.121 1769.04,571.312 1769.96,571.502 1770.88,571.692 1771.79,571.882 1772.71,572.071 1773.63,572.26 \n",
       "  1774.55,572.448 1775.46,572.635 1776.38,572.821 1777.3,573.007 1778.22,573.191 1779.14,573.374 1780.05,573.555 1780.97,573.735 1781.89,573.913 1782.81,574.09 \n",
       "  1783.73,574.265 1784.64,574.438 1785.56,574.609 1786.48,574.778 1787.4,574.945 1788.31,575.109 1789.23,575.272 1790.15,575.431 1791.07,575.588 1791.99,575.743 \n",
       "  1792.9,575.894 1793.82,576.043 1794.74,576.189 1795.66,576.333 1796.57,576.473 1797.49,576.61 1798.41,576.743 1799.33,576.874 1800.25,577.001 1801.16,577.125 \n",
       "  1802.08,577.246 1803,577.363 1803.92,577.476 1804.83,577.586 1805.75,577.692 1806.67,577.795 1807.59,577.893 1808.51,577.988 1809.42,578.08 1810.34,578.167 \n",
       "  1811.26,578.251 1812.18,578.33 1813.09,578.406 1814.01,578.478 1814.93,578.545 1815.85,578.609 1816.77,578.669 1817.68,578.724 1818.6,578.776 1819.52,578.824 \n",
       "  1820.44,578.867 1821.35,578.907 1822.27,578.942 1823.19,578.973 1824.11,579.001 1825.03,579.024 1825.94,579.043 1826.86,579.058 1827.78,579.069 1828.7,579.077 \n",
       "  1829.61,579.08 1830.53,579.079 1831.45,579.074 1832.37,579.065 1833.29,579.053 1834.2,579.037 1835.12,579.016 1836.04,578.992 1836.96,578.965 1837.87,578.933 \n",
       "  1838.79,578.898 1839.71,578.859 1840.63,578.817 1841.55,578.771 1842.46,578.722 1843.38,578.669 1844.3,578.613 1845.22,578.554 1846.13,578.491 1847.05,578.426 \n",
       "  1847.97,578.357 1848.89,578.285 1849.81,578.21 1850.72,578.132 1851.64,578.051 1852.56,577.967 1853.48,577.881 1854.39,577.792 1855.31,577.7 1856.23,577.606 \n",
       "  1857.15,577.51 1858.07,577.411 1858.98,577.31 1859.9,577.206 1860.82,577.101 1861.74,576.993 1862.66,576.884 1863.57,576.772 1864.49,576.659 1865.41,576.544 \n",
       "  1866.33,576.428 1867.24,576.31 1868.16,576.19 1869.08,576.069 1870,575.947 1870.92,575.823 1871.83,575.699 1872.75,575.573 1873.67,575.447 1874.59,575.319 \n",
       "  1875.5,575.191 1876.42,575.062 1877.34,574.933 1878.26,574.803 1879.18,574.673 1880.09,574.542 1881.01,574.411 1881.93,574.28 1882.85,574.149 1883.76,574.017 \n",
       "  1884.68,573.886 1885.6,573.755 1886.52,573.625 1887.44,573.494 1888.35,573.364 1889.27,573.235 1890.19,573.106 1891.11,572.978 1892.02,572.85 1892.94,572.724 \n",
       "  1893.86,572.598 1894.78,572.473 1895.7,572.349 1896.61,572.226 1897.53,572.105 1898.45,571.985 1899.37,571.866 1900.28,571.748 1901.2,571.632 1902.12,571.518 \n",
       "  1903.04,571.405 1903.96,571.293 1904.87,571.184 1905.79,571.076 1906.71,570.97 1907.63,570.866 1908.54,570.764 1909.46,570.664 1910.38,570.566 1911.3,570.47 \n",
       "  1912.22,570.376 1913.13,570.284 1914.05,570.195 1914.97,570.108 1915.89,570.023 1916.8,569.94 1917.72,569.86 1918.64,569.783 1919.56,569.708 1920.48,569.635 \n",
       "  1921.39,569.565 1922.31,569.498 1923.23,569.433 1924.15,569.371 1925.06,569.311 1925.98,569.254 1926.9,569.2 1927.82,569.149 1928.74,569.1 1929.65,569.054 \n",
       "  1930.57,569.011 1931.49,568.97 1932.41,568.933 1933.33,568.898 1934.24,568.866 1935.16,568.837 1936.08,568.81 1937,568.786 1937.91,568.766 1938.83,568.747 \n",
       "  1939.75,568.732 1940.67,568.72 1941.59,568.71 1942.5,568.703 1943.42,568.699 1944.34,568.697 1945.26,568.699 1946.17,568.703 1947.09,568.709 1948.01,568.719 \n",
       "  1948.93,568.731 1949.85,568.745 1950.76,568.763 1951.68,568.782 1952.6,568.805 1953.52,568.829 1954.43,568.857 1955.35,568.887 1956.27,568.919 1957.19,568.953 \n",
       "  1958.11,568.99 1959.02,569.03 1959.94,569.071 1960.86,569.115 1961.78,569.161 1962.69,569.209 1963.61,569.259 1964.53,569.311 1965.45,569.366 1966.37,569.422 \n",
       "  1967.28,569.48 1968.2,569.54 1969.12,569.602 1970.04,569.665 1970.95,569.731 1971.87,569.798 1972.79,569.866 1973.71,569.936 1974.63,570.008 1975.54,570.081 \n",
       "  1976.46,570.156 1977.38,570.231 1978.3,570.309 1979.21,570.387 1980.13,570.466 1981.05,570.547 1981.97,570.629 1982.89,570.711 1983.8,570.795 1984.72,570.879 \n",
       "  1985.64,570.965 1986.56,571.051 1987.47,571.138 1988.39,571.225 1989.31,571.313 1990.23,571.401 1991.15,571.49 1992.06,571.58 1992.98,571.669 1993.9,571.759 \n",
       "  1994.82,571.849 1995.73,571.94 1996.65,572.03 1997.57,572.121 1998.49,572.211 1999.41,572.301 2000.32,572.392 2001.24,572.482 2002.16,572.572 2003.08,572.661 \n",
       "  2004,572.75 2004.91,572.839 2005.83,572.927 2006.75,573.015 2007.67,573.102 2008.58,573.189 2009.5,573.275 2010.42,573.36 2011.34,573.444 2012.26,573.528 \n",
       "  2013.17,573.61 2014.09,573.692 2015.01,573.773 2015.93,573.853 2016.84,573.931 2017.76,574.009 2018.68,574.085 2019.6,574.16 2020.52,574.234 2021.43,574.307 \n",
       "  2022.35,574.379 2023.27,574.449 2024.19,574.517 2025.1,574.585 2026.02,574.65 2026.94,574.715 2027.86,574.777 2028.78,574.839 2029.69,574.898 2030.61,574.956 \n",
       "  2031.53,575.013 2032.45,575.067 2033.36,575.12 2034.28,575.172 2035.2,575.221 2036.12,575.269 2037.04,575.315 2037.95,575.359 2038.87,575.402 2039.79,575.442 \n",
       "  2040.71,575.481 2041.62,575.518 2042.54,575.553 2043.46,575.586 2044.38,575.617 2045.3,575.647 2046.21,575.674 2047.13,575.699 2048.05,575.723 2048.97,575.745 \n",
       "  2049.88,575.764 2050.8,575.782 2051.72,575.798 2052.64,575.812 2053.56,575.824 2054.47,575.834 2055.39,575.842 2056.31,575.848 2057.23,575.852 2058.14,575.855 \n",
       "  2059.06,575.855 2059.98,575.854 2060.9,575.851 2061.82,575.846 2062.73,575.839 2063.65,575.83 2064.57,575.819 2065.49,575.807 2066.4,575.793 2067.32,575.777 \n",
       "  2068.24,575.76 2069.16,575.74 2070.08,575.719 2070.99,575.697 2071.91,575.672 2072.83,575.647 2073.75,575.619 2074.67,575.59 2075.58,575.56 2076.5,575.528 \n",
       "  2077.42,575.494 2078.34,575.459 2079.25,575.423 2080.17,575.385 2081.09,575.346 2082.01,575.305 2082.93,575.264 2083.84,575.221 2084.76,575.177 2085.68,575.131 \n",
       "  2086.6,575.085 2087.51,575.037 2088.43,574.989 2089.35,574.939 2090.27,574.888 2091.19,574.837 2092.1,574.784 2093.02,574.731 2093.94,574.677 2094.86,574.622 \n",
       "  2095.77,574.566 2096.69,574.51 2097.61,574.452 2098.53,574.395 2099.45,574.336 2100.36,574.277 2101.28,574.218 2102.2,574.158 2103.12,574.098 2104.03,574.037 \n",
       "  2104.95,573.976 2105.87,573.914 2106.79,573.853 2107.71,573.791 2108.62,573.729 2109.54,573.667 2110.46,573.605 2111.38,573.542 2112.29,573.48 2113.21,573.417 \n",
       "  2114.13,573.355 2115.05,573.293 2115.97,573.231 2116.88,573.169 2117.8,573.108 2118.72,573.046 2119.64,572.985 2120.55,572.924 2121.47,572.864 2122.39,572.804 \n",
       "  2123.31,572.744 2124.23,572.685 2125.14,572.627 2126.06,572.569 2126.98,572.511 2127.9,572.455 2128.81,572.398 2129.73,572.343 2130.65,572.288 2131.57,572.234 \n",
       "  2132.49,572.181 2133.4,572.128 2134.32,572.077 2135.24,572.026 2136.16,571.976 2137.07,571.927 2137.99,571.879 2138.91,571.832 2139.83,571.786 2140.75,571.741 \n",
       "  2141.66,571.697 2142.58,571.654 2143.5,571.612 2144.42,571.571 2145.33,571.531 2146.25,571.493 2147.17,571.455 2148.09,571.419 2149.01,571.384 2149.92,571.35 \n",
       "  2150.84,571.318 2151.76,571.286 2152.68,571.256 2153.6,571.227 2154.51,571.2 2155.43,571.173 2156.35,571.148 2157.27,571.124 2158.18,571.102 2159.1,571.081 \n",
       "  2160.02,571.061 2160.94,571.042 2161.86,571.025 2162.77,571.009 2163.69,570.995 2164.61,570.982 2165.53,570.97 2166.44,570.959 2167.36,570.95 2168.28,570.942 \n",
       "  2169.2,570.936 2170.12,570.93 2171.03,570.926 2171.95,570.924 2172.87,570.922 2173.79,570.922 2174.7,570.924 2175.62,570.926 2176.54,570.93 2177.46,570.935 \n",
       "  2178.38,570.942 2179.29,570.949 2180.21,570.958 2181.13,570.968 2182.05,570.979 2182.96,570.992 2183.88,571.005 2184.8,571.02 2185.72,571.036 2186.64,571.053 \n",
       "  2187.55,571.071 2188.47,571.09 2189.39,571.111 2190.31,571.132 2191.22,571.154 2192.14,571.178 2193.06,571.202 2193.98,571.227 2194.9,571.254 2195.81,571.281 \n",
       "  2196.73,571.309 2197.65,571.338 2198.57,571.368 2199.48,571.398 2200.4,571.43 2201.32,571.462 2202.24,571.495 2203.16,571.529 2204.07,571.563 2204.99,571.598 \n",
       "  2205.91,571.634 2206.83,571.67 2207.74,571.707 2208.66,571.745 2209.58,571.783 2210.5,571.821 2211.42,571.861 2212.33,571.9 2213.25,571.94 2214.17,571.98 \n",
       "  2215.09,572.021 2216,572.062 2216.92,572.104 2217.84,572.145 2218.76,572.187 2219.68,572.229 2220.59,572.272 2221.51,572.314 2222.43,572.357 2223.35,572.4 \n",
       "  2224.27,572.443 2225.18,572.486 2226.1,572.528 2227.02,572.571 2227.94,572.614 2228.85,572.657 2229.77,572.7 2230.69,572.743 2231.61,572.786 2232.53,572.828 \n",
       "  2233.44,572.87 2234.36,572.912 2235.28,572.954 2236.2,572.996 2237.11,573.037 2238.03,573.078 2238.95,573.119 2239.87,573.159 2240.79,573.199 2241.7,573.238 \n",
       "  2242.62,573.277 2243.54,573.316 2244.46,573.354 2245.37,573.392 2246.29,573.429 2247.21,573.465 2248.13,573.501 2249.05,573.537 2249.96,573.572 2250.88,573.606 \n",
       "  2251.8,573.64 2252.72,573.672 2253.63,573.705 2254.55,573.736 2255.47,573.767 2256.39,573.797 2257.31,573.827 2258.22,573.856 2259.14,573.884 2260.06,573.911 \n",
       "  2260.98,573.937 2261.89,573.963 2262.81,573.987 2263.73,574.011 2264.65,574.035 2265.57,574.057 2266.48,574.078 2267.4,574.099 2268.32,574.118 2269.24,574.137 \n",
       "  2270.15,574.155 2271.07,574.172 2271.99,574.188 2272.91,574.204 2273.83,574.218 2274.74,574.232 2275.66,574.244 2276.58,574.256 2277.5,574.266 2278.41,574.276 \n",
       "  2279.33,574.285 2280.25,574.293 2281.17,574.3 2282.09,574.306 2283,574.311 2283.92,574.316 2284.84,574.319 2285.76,574.322 2286.67,574.323 2287.59,574.324 \n",
       "  2288.51,574.324 2289.43,574.323 2290.35,574.321 2291.26,574.318 2292.18,574.314 2293.1,574.309 2294.02,574.304 2294.94,574.298 2295.85,574.29 2296.77,574.283 \n",
       "  2297.69,574.274 \n",
       "  \"/>\n",
       "<polyline clip-path=\"url(#clip0802)\" style=\"stroke:#e26f46; stroke-width:4; stroke-opacity:1; fill:none\" points=\"\n",
       "  462.105,1388.71 463.023,1388.71 463.941,1388.71 464.858,1388.71 465.776,1388.71 466.694,1388.71 467.612,1388.71 468.529,1388.71 469.447,1388.71 470.365,1388.71 \n",
       "  471.283,1388.71 472.201,1388.71 473.118,1388.71 474.036,1388.71 474.954,1388.71 475.872,1388.71 476.79,1388.71 477.707,1388.71 478.625,1388.71 479.543,1388.71 \n",
       "  480.461,1388.71 481.379,1388.71 482.296,1388.71 483.214,1388.71 484.132,1388.71 485.05,1388.71 485.968,1388.71 486.885,1388.71 487.803,1388.71 488.721,1388.71 \n",
       "  489.639,1388.71 490.556,1388.71 491.474,1388.71 492.392,1388.71 493.31,1388.71 494.228,1388.71 495.145,1388.71 496.063,1388.71 496.981,1388.71 497.899,1388.71 \n",
       "  498.817,1388.71 499.734,1388.71 500.652,1388.71 501.57,1388.71 502.488,1388.71 503.406,1388.71 504.323,1388.71 505.241,1388.71 506.159,1388.71 507.077,1388.71 \n",
       "  507.995,1388.71 508.912,1388.71 509.83,1388.71 510.748,1388.71 511.666,1388.71 512.583,1388.71 513.501,1388.71 514.419,1388.71 515.337,1388.71 516.255,1388.71 \n",
       "  517.172,1388.71 518.09,1388.71 519.008,1388.71 519.926,1388.71 520.844,1388.71 521.761,1388.71 522.679,1388.71 523.597,1388.71 524.515,1388.71 525.433,1388.71 \n",
       "  526.35,1388.71 527.268,1388.71 528.186,1388.71 529.104,1388.71 530.022,1388.71 530.939,1388.71 531.857,1388.71 532.775,1388.71 533.693,1388.71 534.61,1388.71 \n",
       "  535.528,1388.71 536.446,1388.71 537.364,1388.71 538.282,1388.71 539.199,1388.71 540.117,1388.71 541.035,1388.71 541.953,1388.71 542.871,1388.71 543.788,1388.71 \n",
       "  544.706,1388.71 545.624,1388.71 546.542,1388.71 547.46,1388.71 548.377,1388.71 549.295,1388.71 550.213,1388.71 551.131,1388.71 552.049,1388.71 552.966,1388.71 \n",
       "  553.884,1388.71 554.802,1388.71 555.72,1388.71 556.637,1388.71 557.555,1388.71 558.473,1388.71 559.391,1388.71 560.309,1388.71 561.226,1388.71 562.144,1388.71 \n",
       "  563.062,1388.71 563.98,1388.71 564.898,1388.71 565.815,1388.71 566.733,1388.71 567.651,1388.71 568.569,1388.71 569.487,1388.71 570.404,1388.71 571.322,1388.71 \n",
       "  572.24,1388.71 573.158,1388.71 574.076,1388.71 574.993,1388.71 575.911,1388.71 576.829,1388.71 577.747,1388.71 578.664,1388.71 579.582,1388.71 580.5,1388.71 \n",
       "  581.418,1388.71 582.336,1388.71 583.253,1388.71 584.171,1388.71 585.089,1388.71 586.007,1388.71 586.925,1388.71 587.842,1388.71 588.76,1388.71 589.678,1388.71 \n",
       "  590.596,1388.71 591.514,1388.71 592.431,1388.71 593.349,1388.71 594.267,1388.71 595.185,1388.71 596.103,1388.71 597.02,1388.71 597.938,1388.71 598.856,1388.71 \n",
       "  599.774,1388.71 600.691,1388.71 601.609,1388.71 602.527,1388.71 603.445,1388.71 604.363,1388.71 605.28,1388.71 606.198,1388.71 607.116,1388.71 608.034,1388.71 \n",
       "  608.952,1388.71 609.869,1388.71 610.787,1388.71 611.705,1388.71 612.623,1388.71 613.541,1388.71 614.458,1388.71 615.376,1388.71 616.294,1388.71 617.212,1388.71 \n",
       "  618.13,1388.71 619.047,1388.71 619.965,1388.71 620.883,1388.71 621.801,1388.71 622.718,1388.71 623.636,1388.71 624.554,1388.71 625.472,1388.71 626.39,1388.71 \n",
       "  627.307,1388.71 628.225,1388.71 629.143,1388.71 630.061,1388.71 630.979,1388.71 631.896,1388.71 632.814,1388.71 633.732,1388.71 634.65,1388.71 635.568,1388.71 \n",
       "  636.485,1388.71 637.403,1388.71 638.321,1388.71 639.239,1388.71 640.157,1388.71 641.074,1388.71 641.992,1388.71 642.91,1388.71 643.828,1388.71 644.745,1388.71 \n",
       "  645.663,1388.71 646.581,1388.71 647.499,1388.71 648.417,1388.71 649.334,1388.71 650.252,1388.71 651.17,1388.71 652.088,1388.71 653.006,1388.71 653.923,1388.71 \n",
       "  654.841,1388.71 655.759,1388.71 656.677,1388.71 657.595,1388.71 658.512,1388.71 659.43,1388.71 660.348,1388.71 661.266,1388.71 662.184,1388.71 663.101,1388.71 \n",
       "  664.019,1388.71 664.937,1388.71 665.855,1388.71 666.772,1388.71 667.69,1388.71 668.608,1388.71 669.526,1388.71 670.444,1388.71 671.361,1388.71 672.279,1388.71 \n",
       "  673.197,1388.71 674.115,1388.71 675.033,1388.71 675.95,1388.71 676.868,1388.71 677.786,1388.71 678.704,1388.71 679.622,1388.71 680.539,1388.71 681.457,1388.71 \n",
       "  682.375,1388.71 683.293,1388.71 684.211,1388.71 685.128,1388.71 686.046,1388.71 686.964,1388.71 687.882,1388.71 688.799,1388.71 689.717,1388.71 690.635,1388.71 \n",
       "  691.553,1388.71 692.471,1388.71 693.388,1388.71 694.306,1388.71 695.224,1388.71 696.142,1388.71 697.06,1388.71 697.977,1388.71 698.895,1388.71 699.813,1388.71 \n",
       "  700.731,1388.71 701.649,1388.71 702.566,1388.71 703.484,1388.71 704.402,1388.71 705.32,1388.71 706.238,1388.71 707.155,1388.71 708.073,1388.71 708.991,1388.71 \n",
       "  709.909,1388.71 710.826,1388.71 711.744,1388.71 712.662,1388.71 713.58,1388.71 714.498,1388.71 715.415,1388.71 716.333,1388.71 717.251,1388.71 718.169,1388.71 \n",
       "  719.087,1388.71 720.004,1388.71 720.922,1388.71 721.84,1388.71 722.758,1388.71 723.676,1388.71 724.593,1388.71 725.511,1388.71 726.429,1388.71 727.347,1388.71 \n",
       "  728.265,1388.71 729.182,1388.71 730.1,1388.71 731.018,1388.71 731.936,1388.71 732.853,1388.71 733.771,1388.71 734.689,1388.71 735.607,1388.71 736.525,1388.71 \n",
       "  737.442,1388.71 738.36,1388.71 739.278,1388.71 740.196,1388.71 741.114,1388.71 742.031,1388.71 742.949,1388.71 743.867,1388.71 744.785,1388.71 745.703,1388.71 \n",
       "  746.62,1388.71 747.538,1388.71 748.456,1388.71 749.374,1388.71 750.292,1388.71 751.209,1388.71 752.127,1388.71 753.045,1388.71 753.963,1388.71 754.88,1388.71 \n",
       "  755.798,1388.71 756.716,1388.71 757.634,1388.71 758.552,1388.71 759.469,1388.71 760.387,1388.71 761.305,1388.71 762.223,1388.71 763.141,1388.71 764.058,1388.71 \n",
       "  764.976,1388.71 765.894,1388.71 766.812,1388.71 767.73,1388.71 768.647,1388.71 769.565,1388.71 770.483,1388.71 771.401,1388.71 772.319,1388.71 773.236,1388.71 \n",
       "  774.154,1388.71 775.072,1388.71 775.99,1388.71 776.907,1388.71 777.825,1388.71 778.743,1388.71 779.661,1388.71 780.579,1388.71 781.496,1388.71 782.414,1388.71 \n",
       "  783.332,1388.71 784.25,1388.71 785.168,1388.71 786.085,1388.71 787.003,1388.71 787.921,1388.71 788.839,1388.71 789.757,1388.71 790.674,1388.71 791.592,1388.71 \n",
       "  792.51,1388.71 793.428,1388.71 794.346,1388.71 795.263,1388.71 796.181,1388.71 797.099,1388.71 798.017,1388.71 798.934,1388.71 799.852,1388.71 800.77,1388.71 \n",
       "  801.688,1388.71 802.606,1388.71 803.523,1388.71 804.441,1388.71 805.359,1388.71 806.277,1388.71 807.195,1388.71 808.112,1388.71 809.03,1388.71 809.948,1388.71 \n",
       "  810.866,1388.71 811.784,1388.71 812.701,1388.71 813.619,1388.71 814.537,1388.71 815.455,1388.71 816.373,1388.71 817.29,1388.71 818.208,1388.71 819.126,1388.71 \n",
       "  820.044,1388.71 820.962,1388.71 821.879,1388.71 822.797,1388.71 823.715,1388.71 824.633,1388.71 825.55,1388.71 826.468,1388.71 827.386,1388.71 828.304,1388.71 \n",
       "  829.222,1388.71 830.139,1388.71 831.057,1388.71 831.975,1388.71 832.893,1388.71 833.811,1388.71 834.728,1388.71 835.646,1388.71 836.564,1388.71 837.482,1388.71 \n",
       "  838.4,1388.71 839.317,1388.71 840.235,1388.71 841.153,1388.71 842.071,1388.71 842.989,1388.71 843.906,1388.71 844.824,1388.71 845.742,1388.71 846.66,1388.71 \n",
       "  847.577,1388.71 848.495,1388.71 849.413,1388.71 850.331,1388.71 851.249,1388.71 852.166,1388.71 853.084,1388.71 854.002,1388.71 854.92,1388.71 855.838,1388.71 \n",
       "  856.755,1388.71 857.673,1388.71 858.591,1388.71 859.509,1388.71 860.427,1388.71 861.344,1388.71 862.262,1388.71 863.18,1388.71 864.098,1388.71 865.016,1388.71 \n",
       "  865.933,1388.71 866.851,1388.71 867.769,1388.71 868.687,1388.71 869.604,1388.71 870.522,1388.71 871.44,1388.71 872.358,1388.71 873.276,1388.71 874.193,1388.71 \n",
       "  875.111,1388.71 876.029,1388.71 876.947,1388.71 877.865,1388.71 878.782,1388.71 879.7,1388.71 880.618,1388.71 881.536,1388.71 882.454,1388.71 883.371,1388.71 \n",
       "  884.289,1388.71 885.207,1388.71 886.125,1388.71 887.043,1388.71 887.96,1388.71 888.878,1388.71 889.796,1388.71 890.714,1388.71 891.631,1388.71 892.549,1388.71 \n",
       "  893.467,1388.71 894.385,1388.71 895.303,1388.71 896.22,1388.71 897.138,1388.71 898.056,1388.71 898.974,1388.71 899.892,1388.71 900.809,1388.71 901.727,1388.71 \n",
       "  902.645,1388.71 903.563,1388.71 904.481,1388.71 905.398,1388.71 906.316,1388.71 907.234,1388.71 908.152,1388.71 909.07,1388.71 909.987,1388.71 910.905,1388.71 \n",
       "  911.823,1388.71 912.741,1388.71 913.658,1388.71 914.576,1388.71 915.494,1388.71 916.412,1388.71 917.33,1388.71 918.247,1388.71 919.165,1388.71 920.083,1388.71 \n",
       "  921.001,1388.71 921.919,1388.71 922.836,1388.71 923.754,1388.71 924.672,1388.71 925.59,1388.71 926.508,1388.71 927.425,1388.71 928.343,1388.71 929.261,1388.71 \n",
       "  930.179,1388.71 931.097,1388.71 932.014,1388.71 932.932,1388.71 933.85,1388.71 934.768,1388.71 935.685,1388.71 936.603,1388.71 937.521,1388.71 938.439,1388.71 \n",
       "  939.357,1388.71 940.274,1388.71 941.192,1388.71 942.11,1388.71 943.028,1388.71 943.946,1388.71 944.863,1388.71 945.781,1388.71 946.699,1388.71 947.617,1388.71 \n",
       "  948.535,1388.71 949.452,1388.71 950.37,1388.71 951.288,1388.71 952.206,1388.71 953.124,1388.71 954.041,1388.71 954.959,1388.71 955.877,1388.71 956.795,1388.71 \n",
       "  957.712,1388.71 958.63,1388.71 959.548,1388.71 960.466,1388.71 961.384,1388.71 962.301,1388.71 963.219,1388.71 964.137,1388.71 965.055,1388.71 965.973,1388.71 \n",
       "  966.89,1388.71 967.808,1388.71 968.726,1388.71 969.644,1388.71 970.562,1388.71 971.479,1388.71 972.397,1388.71 973.315,1388.71 974.233,1388.71 975.151,1388.71 \n",
       "  976.068,1388.71 976.986,1388.71 977.904,1388.71 978.822,1388.71 979.739,1388.71 980.657,1388.71 981.575,1388.71 982.493,1388.71 983.411,1388.71 984.328,1388.71 \n",
       "  985.246,1388.71 986.164,1388.71 987.082,1388.71 988,1388.71 988.917,1388.71 989.835,1388.71 990.753,1388.71 991.671,1388.71 992.589,1388.71 993.506,1388.71 \n",
       "  994.424,1388.71 995.342,1388.71 996.26,1388.71 997.178,1388.71 998.095,1388.71 999.013,1388.71 999.931,1388.71 1000.85,1388.71 1001.77,1388.71 1002.68,1388.71 \n",
       "  1003.6,1388.71 1004.52,1388.71 1005.44,1388.71 1006.36,1388.71 1007.27,1388.71 1008.19,1388.71 1009.11,1388.71 1010.03,1388.71 1010.94,1388.71 1011.86,1388.71 \n",
       "  1012.78,1388.71 1013.7,1388.71 1014.62,1388.71 1015.53,1388.71 1016.45,1388.71 1017.37,1388.71 1018.29,1388.71 1019.2,1388.71 1020.12,1388.71 1021.04,1388.71 \n",
       "  1021.96,1388.71 1022.88,1388.71 1023.79,1388.71 1024.71,1388.71 1025.63,1388.71 1026.55,1388.71 1027.46,1388.71 1028.38,1388.71 1029.3,1388.71 1030.22,1388.71 \n",
       "  1031.14,1388.71 1032.05,1388.71 1032.97,1388.71 1033.89,1388.71 1034.81,1388.71 1035.72,1388.71 1036.64,1388.71 1037.56,1388.71 1038.48,1388.71 1039.4,1388.71 \n",
       "  1040.31,1388.71 1041.23,1388.71 1042.15,1388.71 1043.07,1388.71 1043.98,1388.71 1044.9,1388.71 1045.82,1388.71 1046.74,1388.71 1047.66,1388.71 1048.57,1388.71 \n",
       "  1049.49,1388.71 1050.41,1388.71 1051.33,1388.71 1052.25,1388.71 1053.16,1388.71 1054.08,1388.71 1055,1388.71 1055.92,1388.71 1056.83,1388.71 1057.75,1388.71 \n",
       "  1058.67,1388.71 1059.59,1388.71 1060.51,1388.71 1061.42,1388.71 1062.34,1388.71 1063.26,1388.71 1064.18,1388.71 1065.09,1388.71 1066.01,1388.71 1066.93,1388.71 \n",
       "  1067.85,1388.71 1068.77,1388.71 1069.68,1388.71 1070.6,1388.71 1071.52,1388.71 1072.44,1388.71 1073.35,1388.71 1074.27,1388.71 1075.19,1388.71 1076.11,1388.71 \n",
       "  1077.03,1388.71 1077.94,1388.71 1078.86,1388.71 1079.78,1388.71 1080.7,1388.71 1081.61,1388.71 1082.53,1388.71 1083.45,1388.71 1084.37,1388.71 1085.29,1388.71 \n",
       "  1086.2,1388.71 1087.12,1388.71 1088.04,1388.71 1088.96,1388.71 1089.87,1388.71 1090.79,1388.71 1091.71,1388.71 1092.63,1388.71 1093.55,1388.71 1094.46,1388.71 \n",
       "  1095.38,1388.71 1096.3,1388.71 1097.22,1388.71 1098.13,1388.71 1099.05,1388.71 1099.97,1388.71 1100.89,1388.71 1101.81,1388.71 1102.72,1388.71 1103.64,1388.71 \n",
       "  1104.56,1388.71 1105.48,1388.71 1106.39,1388.71 1107.31,1388.71 1108.23,1388.71 1109.15,1388.71 1110.07,1388.71 1110.98,1388.71 1111.9,1388.71 1112.82,1388.71 \n",
       "  1113.74,1388.71 1114.65,1388.71 1115.57,1388.71 1116.49,1388.71 1117.41,1388.71 1118.33,1388.71 1119.24,1388.71 1120.16,1388.71 1121.08,1388.71 1122,1388.71 \n",
       "  1122.91,1388.71 1123.83,1388.71 1124.75,1388.71 1125.67,1388.71 1126.59,1388.71 1127.5,1388.71 1128.42,1388.71 1129.34,1388.71 1130.26,1388.71 1131.18,1388.71 \n",
       "  1132.09,1388.71 1133.01,1388.71 1133.93,1388.71 1134.85,1388.71 1135.76,1388.71 1136.68,1388.71 1137.6,1388.71 1138.52,1388.71 1139.44,1388.71 1140.35,1388.71 \n",
       "  1141.27,1388.71 1142.19,1388.71 1143.11,1388.71 1144.02,1388.71 1144.94,1388.71 1145.86,1388.71 1146.78,1388.71 1147.7,1388.71 1148.61,1388.71 1149.53,1388.71 \n",
       "  1150.45,1388.71 1151.37,1388.71 1152.28,1388.71 1153.2,1388.71 1154.12,1388.71 1155.04,1388.71 1155.96,1388.71 1156.87,1388.71 1157.79,1388.71 1158.71,1388.71 \n",
       "  1159.63,1388.71 1160.54,1388.71 1161.46,1388.71 1162.38,1388.71 1163.3,1388.71 1164.22,1388.71 1165.13,1388.71 1166.05,1388.71 1166.97,1388.71 1167.89,1388.71 \n",
       "  1168.8,1388.71 1169.72,1388.71 1170.64,1388.71 1171.56,1388.71 1172.48,1388.71 1173.39,1388.71 1174.31,1388.71 1175.23,1388.71 1176.15,1388.71 1177.06,1388.71 \n",
       "  1177.98,1388.71 1178.9,1388.71 1179.82,1388.71 1180.74,1388.71 1181.65,1388.71 1182.57,1388.71 1183.49,1388.71 1184.41,1388.71 1185.32,1388.71 1186.24,1388.71 \n",
       "  1187.16,1388.71 1188.08,1388.71 1189,1388.71 1189.91,1388.71 1190.83,1388.71 1191.75,1388.71 1192.67,1388.71 1193.58,1388.71 1194.5,1388.71 1195.42,1388.71 \n",
       "  1196.34,1388.71 1197.26,1388.71 1198.17,1388.71 1199.09,1388.71 1200.01,1388.71 1200.93,1388.71 1201.85,1388.71 1202.76,1388.71 1203.68,1388.71 1204.6,1388.71 \n",
       "  1205.52,1388.71 1206.43,1388.71 1207.35,1388.71 1208.27,1388.71 1209.19,1388.71 1210.11,1388.71 1211.02,1388.71 1211.94,1388.71 1212.86,1388.71 1213.78,1388.71 \n",
       "  1214.69,1388.71 1215.61,1388.71 1216.53,1388.71 1217.45,1388.71 1218.37,1388.71 1219.28,1388.71 1220.2,1388.71 1221.12,1388.71 1222.04,1388.71 1222.95,1388.71 \n",
       "  1223.87,1388.71 1224.79,1388.71 1225.71,1388.71 1226.63,1388.71 1227.54,1388.71 1228.46,1388.71 1229.38,1388.71 1230.3,1388.71 1231.21,1388.71 1232.13,1388.71 \n",
       "  1233.05,1388.71 1233.97,1388.71 1234.89,1388.71 1235.8,1388.71 1236.72,1388.71 1237.64,1388.71 1238.56,1388.71 1239.47,1388.71 1240.39,1388.71 1241.31,1388.71 \n",
       "  1242.23,1388.71 1243.15,1388.71 1244.06,1388.71 1244.98,1388.71 1245.9,1388.71 1246.82,1388.71 1247.73,1388.71 1248.65,1388.71 1249.57,1388.71 1250.49,1388.71 \n",
       "  1251.41,1388.71 1252.32,1388.71 1253.24,1388.71 1254.16,1388.71 1255.08,1388.71 1255.99,1388.71 1256.91,1388.71 1257.83,1388.71 1258.75,1388.71 1259.67,1388.71 \n",
       "  1260.58,1388.71 1261.5,1388.71 1262.42,1388.71 1263.34,1388.71 1264.25,1388.71 1265.17,1388.71 1266.09,1388.71 1267.01,1388.71 1267.93,1388.71 1268.84,1388.71 \n",
       "  1269.76,1388.71 1270.68,1388.71 1271.6,1388.71 1272.52,1388.71 1273.43,1388.71 1274.35,1388.71 1275.27,1388.71 1276.19,1388.71 1277.1,1388.71 1278.02,1388.71 \n",
       "  1278.94,1388.71 1279.86,1388.71 1280.78,1388.71 1281.69,1388.71 1282.61,1388.71 1283.53,1388.71 1284.45,1388.71 1285.36,1388.71 1286.28,1388.71 1287.2,1388.71 \n",
       "  1288.12,1388.71 1289.04,1388.71 1289.95,1388.71 1290.87,1388.71 1291.79,1388.71 1292.71,1388.71 1293.62,1388.71 1294.54,1388.71 1295.46,1388.71 1296.38,1388.71 \n",
       "  1297.3,1388.71 1298.21,1388.71 1299.13,1388.71 1300.05,1388.71 1300.97,1388.71 1301.88,1388.71 1302.8,1388.71 1303.72,1388.71 1304.64,1388.71 1305.56,1388.71 \n",
       "  1306.47,1388.71 1307.39,1388.71 1308.31,1388.71 1309.23,1388.71 1310.14,1388.71 1311.06,1388.71 1311.98,1388.71 1312.9,1388.71 1313.82,1388.71 1314.73,1388.71 \n",
       "  1315.65,1388.71 1316.57,1388.71 1317.49,1388.71 1318.4,1388.71 1319.32,1388.71 1320.24,1388.71 1321.16,1388.71 1322.08,1388.71 1322.99,1388.71 1323.91,1388.71 \n",
       "  1324.83,1388.71 1325.75,1388.71 1326.66,1388.71 1327.58,1388.71 1328.5,1388.71 1329.42,1388.71 1330.34,1388.71 1331.25,1388.71 1332.17,1388.71 1333.09,1388.71 \n",
       "  1334.01,1388.71 1334.92,1388.71 1335.84,1388.71 1336.76,1388.71 1337.68,1388.71 1338.6,1388.71 1339.51,1388.71 1340.43,1388.71 1341.35,1388.71 1342.27,1388.71 \n",
       "  1343.18,1388.71 1344.1,1388.71 1345.02,1388.71 1345.94,1388.71 1346.86,1388.71 1347.77,1388.71 1348.69,1388.71 1349.61,1388.71 1350.53,1388.71 1351.45,1388.71 \n",
       "  1352.36,1388.71 1353.28,1388.71 1354.2,1388.71 1355.12,1388.71 1356.03,1388.71 1356.95,1388.71 1357.87,1388.71 1358.79,1388.71 1359.71,1388.71 1360.62,1388.71 \n",
       "  1361.54,1388.71 1362.46,1388.71 1363.38,1388.71 1364.29,1388.71 1365.21,1388.71 1366.13,1388.71 1367.05,1388.71 1367.97,1388.71 1368.88,1388.71 1369.8,1388.71 \n",
       "  1370.72,1388.71 1371.64,1388.71 1372.55,1388.71 1373.47,1388.71 1374.39,1388.71 1375.31,1388.71 1376.23,1388.71 1377.14,1388.71 1378.06,1388.71 1378.98,1388.71 \n",
       "  1379.9,1388.71 1380.81,1388.71 1381.73,1388.71 1382.65,1388.71 1383.57,1388.71 1384.49,1388.71 1385.4,1388.71 1386.32,1388.71 1387.24,1388.71 1388.16,1388.71 \n",
       "  1389.07,1388.71 1389.99,1388.71 1390.91,1388.71 1391.83,1388.71 1392.75,1388.71 1393.66,1388.71 1394.58,1388.71 1395.5,1388.71 1396.42,1388.71 1397.33,1388.71 \n",
       "  1398.25,1388.71 1399.17,1388.71 1400.09,1388.71 1401.01,1388.71 1401.92,1388.71 1402.84,1388.71 1403.76,1388.71 1404.68,1388.71 1405.59,1388.71 1406.51,1388.71 \n",
       "  1407.43,1388.71 1408.35,1388.71 1409.27,1388.71 1410.18,1388.71 1411.1,1388.71 1412.02,1388.71 1412.94,1388.71 1413.85,1388.71 1414.77,1388.71 1415.69,1388.71 \n",
       "  1416.61,1388.71 1417.53,1388.71 1418.44,1388.71 1419.36,1388.71 1420.28,1388.71 1421.2,1388.71 1422.12,1388.71 1423.03,1388.71 1423.95,1388.71 1424.87,1388.71 \n",
       "  1425.79,1388.71 1426.7,1388.71 1427.62,1388.71 1428.54,1388.71 1429.46,1388.71 1430.38,1388.71 1431.29,1388.71 1432.21,1388.71 1433.13,1388.71 1434.05,1388.71 \n",
       "  1434.96,1388.71 1435.88,1388.71 1436.8,1388.71 1437.72,1388.71 1438.64,1388.71 1439.55,1388.71 1440.47,1388.71 1441.39,1388.71 1442.31,1388.71 1443.22,1388.71 \n",
       "  1444.14,1388.71 1445.06,1388.71 1445.98,1388.71 1446.9,1388.71 1447.81,1388.71 1448.73,1388.71 1449.65,1388.71 1450.57,1388.71 1451.48,1388.71 1452.4,1388.71 \n",
       "  1453.32,1388.71 1454.24,1388.71 1455.16,1388.71 1456.07,1388.71 1456.99,1388.71 1457.91,1388.71 1458.83,1388.71 1459.74,1388.71 1460.66,1388.71 1461.58,1388.71 \n",
       "  1462.5,1388.71 1463.42,1388.71 1464.33,1388.71 1465.25,1388.71 1466.17,1388.71 1467.09,1388.71 1468,1388.71 1468.92,1388.71 1469.84,1388.71 1470.76,1388.71 \n",
       "  1471.68,1388.71 1472.59,1388.71 1473.51,1388.71 1474.43,1388.71 1475.35,1388.71 1476.26,1388.71 1477.18,1388.71 1478.1,1388.71 1479.02,1388.71 1479.94,1388.71 \n",
       "  1480.85,1388.71 1481.77,1388.71 1482.69,1388.71 1483.61,1388.71 1484.52,1388.71 1485.44,1388.71 1486.36,1388.71 1487.28,1388.71 1488.2,1388.71 1489.11,1388.71 \n",
       "  1490.03,1388.71 1490.95,1388.71 1491.87,1388.71 1492.79,1388.71 1493.7,1388.71 1494.62,1388.71 1495.54,1388.71 1496.46,1388.71 1497.37,1388.71 1498.29,1388.71 \n",
       "  1499.21,1388.71 1500.13,1388.71 1501.05,1388.71 1501.96,1388.71 1502.88,1388.71 1503.8,1388.71 1504.72,1388.71 1505.63,1388.71 1506.55,1388.71 1507.47,1388.71 \n",
       "  1508.39,1388.71 1509.31,1388.71 1510.22,1388.71 1511.14,1388.71 1512.06,1388.71 1512.98,1388.71 1513.89,1388.71 1514.81,1388.71 1515.73,1388.71 1516.65,1388.71 \n",
       "  1517.57,1388.71 1518.48,1388.71 1519.4,1388.71 1520.32,1388.71 1521.24,1388.71 1522.15,1388.71 1523.07,1388.71 1523.99,1388.71 1524.91,1388.71 1525.83,1388.71 \n",
       "  1526.74,1388.71 1527.66,1388.71 1528.58,1388.71 1529.5,1388.71 1530.41,1388.71 1531.33,1388.71 1532.25,1388.71 1533.17,1388.71 1534.09,1388.71 1535,1388.71 \n",
       "  1535.92,1388.71 1536.84,1388.71 1537.76,1388.71 1538.67,1388.71 1539.59,1388.71 1540.51,1388.71 1541.43,1388.71 1542.35,1388.71 1543.26,1388.71 1544.18,1388.71 \n",
       "  1545.1,1388.71 1546.02,1388.71 1546.93,1388.71 1547.85,1388.71 1548.77,1388.71 1549.69,1388.71 1550.61,1388.71 1551.52,1388.71 1552.44,1388.71 1553.36,1388.71 \n",
       "  1554.28,1388.71 1555.19,1388.71 1556.11,1388.71 1557.03,1388.71 1557.95,1388.71 1558.87,1388.71 1559.78,1388.71 1560.7,1388.71 1561.62,1388.71 1562.54,1388.71 \n",
       "  1563.46,1388.71 1564.37,1388.71 1565.29,1388.71 1566.21,1388.71 1567.13,1388.71 1568.04,1388.71 1568.96,1388.71 1569.88,1388.71 1570.8,1388.71 1571.72,1388.71 \n",
       "  1572.63,1388.71 1573.55,1388.71 1574.47,1388.71 1575.39,1388.71 1576.3,1388.71 1577.22,1388.71 1578.14,1388.71 1579.06,1388.71 1579.98,1388.71 1580.89,1388.71 \n",
       "  1581.81,1388.71 1582.73,1388.71 1583.65,1388.71 1584.56,1388.71 1585.48,1388.71 1586.4,1388.71 1587.32,1388.71 1588.24,1388.71 1589.15,1388.71 1590.07,1388.71 \n",
       "  1590.99,1388.71 1591.91,1388.71 1592.82,1388.71 1593.74,1388.71 1594.66,1388.71 1595.58,1388.71 1596.5,1388.71 1597.41,1388.71 1598.33,1388.71 1599.25,1388.71 \n",
       "  1600.17,1388.71 1601.08,1388.71 1602,1388.71 1602.92,1388.71 1603.84,1388.71 1604.76,1388.71 1605.67,1388.71 1606.59,1388.71 1607.51,1388.71 1608.43,1388.71 \n",
       "  1609.34,1388.71 1610.26,1388.71 1611.18,1388.71 1612.1,1388.71 1613.02,1388.71 1613.93,1388.71 1614.85,1388.71 1615.77,1388.71 1616.69,1388.71 1617.6,1388.71 \n",
       "  1618.52,1388.71 1619.44,1388.71 1620.36,1388.71 1621.28,1388.71 1622.19,1388.71 1623.11,1388.71 1624.03,1388.71 1624.95,1388.71 1625.86,1388.71 1626.78,1388.71 \n",
       "  1627.7,1388.71 1628.62,1388.71 1629.54,1388.71 1630.45,1388.71 1631.37,1388.71 1632.29,1388.71 1633.21,1388.71 1634.12,1388.71 1635.04,1388.71 1635.96,1388.71 \n",
       "  1636.88,1388.71 1637.8,1388.71 1638.71,1388.71 1639.63,1388.71 1640.55,1388.71 1641.47,1388.71 1642.39,1388.71 1643.3,1388.71 1644.22,1388.71 1645.14,1388.71 \n",
       "  1646.06,1388.71 1646.97,1388.71 1647.89,1388.71 1648.81,1388.71 1649.73,1388.71 1650.65,1388.71 1651.56,1388.71 1652.48,1388.71 1653.4,1388.71 1654.32,1388.71 \n",
       "  1655.23,1388.71 1656.15,1388.71 1657.07,1388.71 1657.99,1388.71 1658.91,1388.71 1659.82,1388.71 1660.74,1388.71 1661.66,1388.71 1662.58,1388.71 1663.49,1388.71 \n",
       "  1664.41,1388.71 1665.33,1388.71 1666.25,1388.71 1667.17,1388.71 1668.08,1388.71 1669,1388.71 1669.92,1388.71 1670.84,1388.71 1671.75,1388.71 1672.67,1388.71 \n",
       "  1673.59,1388.71 1674.51,1388.71 1675.43,1388.71 1676.34,1388.71 1677.26,1388.71 1678.18,1388.71 1679.1,1388.71 1680.01,1388.71 1680.93,1388.71 1681.85,1388.71 \n",
       "  1682.77,1388.71 1683.69,1388.71 1684.6,1388.71 1685.52,1388.71 1686.44,1388.71 1687.36,1388.71 1688.27,1388.71 1689.19,1388.71 1690.11,1388.71 1691.03,1388.71 \n",
       "  1691.95,1388.71 1692.86,1388.71 1693.78,1388.71 1694.7,1388.71 1695.62,1388.71 1696.53,1388.71 1697.45,1388.71 1698.37,1388.71 1699.29,1388.71 1700.21,1388.71 \n",
       "  1701.12,1388.71 1702.04,1388.71 1702.96,1388.71 1703.88,1388.71 1704.79,1388.71 1705.71,1388.71 1706.63,1388.71 1707.55,1388.71 1708.47,1388.71 1709.38,1388.71 \n",
       "  1710.3,1388.71 1711.22,1388.71 1712.14,1388.71 1713.06,1388.71 1713.97,1388.71 1714.89,1388.71 1715.81,1388.71 1716.73,1388.71 1717.64,1388.71 1718.56,1388.71 \n",
       "  1719.48,1388.71 1720.4,1388.71 1721.32,1388.71 1722.23,1388.71 1723.15,1388.71 1724.07,1388.71 1724.99,1388.71 1725.9,1388.71 1726.82,1388.71 1727.74,1388.71 \n",
       "  1728.66,1388.71 1729.58,1388.71 1730.49,1388.71 1731.41,1388.71 1732.33,1388.71 1733.25,1388.71 1734.16,1388.71 1735.08,1388.71 1736,1388.71 1736.92,1388.71 \n",
       "  1737.84,1388.71 1738.75,1388.71 1739.67,1388.71 1740.59,1388.71 1741.51,1388.71 1742.42,1388.71 1743.34,1388.71 1744.26,1388.71 1745.18,1388.71 1746.1,1388.71 \n",
       "  1747.01,1388.71 1747.93,1388.71 1748.85,1388.71 1749.77,1388.71 1750.68,1388.71 1751.6,1388.71 1752.52,1388.71 1753.44,1388.71 1754.36,1388.71 1755.27,1388.71 \n",
       "  1756.19,1388.71 1757.11,1388.71 1758.03,1388.71 1758.94,1388.71 1759.86,1388.71 1760.78,1388.71 1761.7,1388.71 1762.62,1388.71 1763.53,1388.71 1764.45,1388.71 \n",
       "  1765.37,1388.71 1766.29,1388.71 1767.2,1388.71 1768.12,1388.71 1769.04,1388.71 1769.96,1388.71 1770.88,1388.71 1771.79,1388.71 1772.71,1388.71 1773.63,1388.71 \n",
       "  1774.55,1388.71 1775.46,1388.71 1776.38,1388.71 1777.3,1388.71 1778.22,1388.71 1779.14,1388.71 1780.05,1388.71 1780.97,1388.71 1781.89,1388.71 1782.81,1388.71 \n",
       "  1783.73,1388.71 1784.64,1388.71 1785.56,1388.71 1786.48,1388.71 1787.4,1388.71 1788.31,1388.71 1789.23,1388.71 1790.15,1388.71 1791.07,1388.71 1791.99,1388.71 \n",
       "  1792.9,1388.71 1793.82,1388.71 1794.74,1388.71 1795.66,1388.71 1796.57,1388.71 1797.49,1388.71 1798.41,1388.71 1799.33,1388.71 1800.25,1388.71 1801.16,1388.71 \n",
       "  1802.08,1388.71 1803,1388.71 1803.92,1388.71 1804.83,1388.71 1805.75,1388.71 1806.67,1388.71 1807.59,1388.71 1808.51,1388.71 1809.42,1388.71 1810.34,1388.71 \n",
       "  1811.26,1388.71 1812.18,1388.71 1813.09,1388.71 1814.01,1388.71 1814.93,1388.71 1815.85,1388.71 1816.77,1388.71 1817.68,1388.71 1818.6,1388.71 1819.52,1388.71 \n",
       "  1820.44,1388.71 1821.35,1388.71 1822.27,1388.71 1823.19,1388.71 1824.11,1388.71 1825.03,1388.71 1825.94,1388.71 1826.86,1388.71 1827.78,1388.71 1828.7,1388.71 \n",
       "  1829.61,1388.71 1830.53,1388.71 1831.45,1388.71 1832.37,1388.71 1833.29,1388.71 1834.2,1388.71 1835.12,1388.71 1836.04,1388.71 1836.96,1388.71 1837.87,1388.71 \n",
       "  1838.79,1388.71 1839.71,1388.71 1840.63,1388.71 1841.55,1388.71 1842.46,1388.71 1843.38,1388.71 1844.3,1388.71 1845.22,1388.71 1846.13,1388.71 1847.05,1388.71 \n",
       "  1847.97,1388.71 1848.89,1388.71 1849.81,1388.71 1850.72,1388.71 1851.64,1388.71 1852.56,1388.71 1853.48,1388.71 1854.39,1388.71 1855.31,1388.71 1856.23,1388.71 \n",
       "  1857.15,1388.71 1858.07,1388.71 1858.98,1388.71 1859.9,1388.71 1860.82,1388.71 1861.74,1388.71 1862.66,1388.71 1863.57,1388.71 1864.49,1388.71 1865.41,1388.71 \n",