SHARE
TWEET

Untitled

a guest Oct 30th, 2015 120 Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. data = {{30.`, 7.203774779333104`*^7}, {30.1`,
  2.  
  3.    7.252478307216305`*^7}, {30.2`, 7.301185887293766`*^7}, {30.3`,
  4.  
  5.    7.349167691889875`*^7}, {30.4`, 7.399421616562334`*^7}, {30.5`,
  6.  
  7.    7.448534216744526`*^7}, {30.6`, 7.497201601371181`*^7}, {30.7`,
  8.  
  9.    7.547000023885964`*^7}, {30.8`, 7.597276907993355`*^7}, {30.9`,
  10.  
  11.    7.646281349372171`*^7}, {31.`, 7.696372067480722`*^7}, {31.1`,
  12.  
  13.    7.746546829656683`*^7}, {31.2`, 7.7972789507134`*^7}, {31.3`,
  14.  
  15.    7.848877393695667`*^7}, {31.4`, 7.900242156431046`*^7}, {31.5`,
  16.  
  17.    7.951206918431002`*^7}, {31.6`, 8.001578238132471`*^7}, {31.7`,
  18.  
  19.    8.054842269392538`*^7}, {31.8`, 8.105356380942908`*^7}, {31.9`,
  20.  
  21.    8.158952975048122`*^7}, {32.`, 8.210701487498114`*^7}, {32.1`,
  22.  
  23.    8.262426654586336`*^7}, {32.2`, 8.316928219055673`*^7}, {32.3`,
  24.  
  25.    8.368938097536127`*^7}, {32.4`, 8.423644950100565`*^7}, {32.5`,
  26.  
  27.    8.476303956033066`*^7}, {32.6`, 8.530743653126228`*^7}, {32.7`,
  28.  
  29.    8.58526662435352`*^7}, {32.8`, 8.639106328655201`*^7}, {32.9`,
  30.  
  31.    8.694213360985236`*^7}, {33.`, 8.749825197687946`*^7}, {33.1`,
  32.  
  33.    8.804845797051753`*^7}, {33.2`, 8.86127965266481`*^7}, {33.3`,
  34.  
  35.    8.916119638447481`*^7}, {33.4`, 8.973660954401416`*^7}, {33.5`,
  36.  
  37.    9.03107757190635`*^7}, {33.6`, 9.086740013475901`*^7}, {33.7`,
  38.  
  39.    9.144590855277804`*^7}, {33.8`, 9.20239409597612`*^7}, {33.9`,
  40.  
  41.    9.262266607977717`*^7}, {34.`, 9.320018565458618`*^7}, {34.1`,
  42.  
  43.    9.37937737178022`*^7}, {34.2`, 9.439758446742903`*^7}, {34.3`,
  44.  
  45.    9.50100704030957`*^7}, {34.4`, 9.561379139501218`*^7}, {34.5`,
  46.  
  47.    9.624574524551912`*^7}, {34.6`, 9.686430795145366`*^7}, {34.7`,
  48.  
  49.    9.75086589796861`*^7}, {34.8`, 9.81462740078182`*^7}, {34.9`,
  50.  
  51.    9.881039421974972`*^7}, {35.`, 9.947494151687416`*^7}, {35.1`,
  52.  
  53.    1.0017916095007841`*^8}, {35.2`, 1.0088396460330561`*^8}, {35.3`,
  54.  
  55.    1.0161647225364536`*^8}, {35.4`, 1.0236887098190545`*^8}, {35.5`,
  56.  
  57.    1.031659269882836`*^8}, {35.6`, 1.0399947109283426`*^8}, {35.7`,
  58.  
  59.    1.0484471814162257`*^8}, {35.8`, 1.0575920267025879`*^8}, {35.9`,
  60.  
  61.    1.0671975088920957`*^8}, {36.`, 1.0773376020462047`*^8}, {36.1`,
  62.  
  63.    1.0880155404222648`*^8}, {36.2`, 1.099497050517212`*^8}, {36.3`,
  64.  
  65.    1.1112743434496064`*^8}, {36.4`, 1.1240604531195396`*^8}, {36.5`,
  66.  
  67.    1.1371877659067413`*^8}, {36.6`, 1.1510642888135706`*^8}, {36.7`,
  68.  
  69.    1.1653752334186544`*^8}, {36.8`, 1.1801996007894506`*^8}, {36.9`,
  70.  
  71.    1.1953209974397758`*^8}, {37.`, 1.2106263658229922`*^8}, {37.1`,
  72.  
  73.    1.2261959910837291`*^8}, {37.2`, 1.2417161846062526`*^8}, {37.3`,
  74.  
  75.    1.2570098354724145`*^8}, {37.4`, 1.2718312442073257`*^8}, {37.5`,
  76.  
  77.    1.2863627654768747`*^8}, {37.6`, 1.3000757635085316`*^8}, {37.7`,
  78.  
  79.    1.3128890519390689`*^8}, {37.8`, 1.3248453608852671`*^8}, {37.9`,
  80.  
  81.    1.3356713306223537`*^8}, {38.`, 1.3455798479605815`*^8}, {38.1`,
  82.  
  83.    1.3544232938307893`*^8}, {38.2`, 1.36207984344155`*^8}, {38.3`,
  84.  
  85.    1.3689957727896827`*^8}, {38.4`, 1.3748396720574304`*^8}, {38.5`,
  86.  
  87.    1.379862589688045`*^8}, {38.6`, 1.3841481591848132`*^8}, {38.7`,
  88.  
  89.    1.3882704044962057`*^8}, {38.8`, 1.3918103353564966`*^8}, {38.9`,
  90.  
  91.    1.3951091445327854`*^8}, {39.`, 1.398770618762684`*^8}, {39.1`,
  92.  
  93.    1.402328987082155`*^8}, {39.2`, 1.406331752775201`*^8}, {39.3`,
  94.  
  95.    1.4107499098296505`*^8}, {39.4`, 1.4157517560517234`*^8}, {39.5`,
  96.  
  97.    1.42135249548077`*^8}, {39.6`, 1.427619104561011`*^8}, {39.7`,
  98.  
  99.    1.434620741473293`*^8}, {39.8`, 1.4425550011580953`*^8}, {39.9`,
  100.  
  101.    1.451080997372032`*^8}, {40.`, 1.460274391829685`*^8}, {40.1`,
  102.  
  103.    1.4704059019867703`*^8}, {40.2`, 1.4808947956960264`*^8}, {40.3`,
  104.  
  105.    1.4922586105683982`*^8}, {40.4`, 1.5041833824736422`*^8}, {40.5`,
  106.  
  107.    1.516552507324408`*^8}, {40.6`, 1.5296589978234124`*^8}, {40.7`,
  108.  
  109.    1.5429365109355915`*^8}, {40.8`, 1.5567547480488497`*^8}, {40.9`,
  110.  
  111.    1.570980903198671`*^8}, {41.`, 1.5855118369710433`*^8}, {41.1`,
  112.  
  113.    1.6006695009150204`*^8}, {41.2`, 1.615732454206531`*^8}, {41.3`,
  114.  
  115.    1.6314217052281234`*^8}, {41.4`, 1.6471425348362955`*^8}, {41.5`,
  116.  
  117.    1.663309180151197`*^8}, {41.6`, 1.6796406025811493`*^8}, {41.7`,
  118.  
  119.    1.6962553449974087`*^8}, {41.8`, 1.7130733761090565`*^8}, {41.9`,
  120.  
  121.    1.730336542829305`*^8}, {42.`, 1.7475357107690114`*^8}, {42.1`,
  122.  
  123.    1.765175115963545`*^8}, {42.2`, 1.7830670596485522`*^8}, {42.3`,
  124.  
  125.    1.801123980348885`*^8}, {42.4`, 1.8193151350275412`*^8}, {42.5`,
  126.  
  127.    1.837714201279694`*^8}, {42.6`, 1.8563922543349117`*^8}, {42.7`,
  128.  
  129.    1.875468295478556`*^8}, {42.8`, 1.8943975419183657`*^8}, {42.9`,
  130.  
  131.    1.9138082744216394`*^8}, {43.`, 1.9332794959105176`*^8}, {43.1`,
  132.  
  133.    1.9529416561835557`*^8}, {43.2`, 1.9729174340013802`*^8}, {43.3`,
  134.  
  135.    1.9930558684445086`*^8}, {43.4`, 2.0132774300021586`*^8}, {43.5`,
  136.  
  137.    2.0335974216727114`*^8}, {43.6`, 2.054015009494567`*^8}, {43.7`,
  138.  
  139.    2.0747153005369624`*^8}, {43.8`, 2.0957196847553968`*^8}, {43.9`,
  140.  
  141.    2.1167140344728222`*^8}, {44.`, 2.137718800927287`*^8}, {44.1`,
  142.  
  143.    2.159028699677915`*^8}, {44.2`, 2.180355476995228`*^8}, {44.3`,
  144.  
  145.    2.2018497451847026`*^8}, {44.4`, 2.2231831238825506`*^8}, {44.5`,
  146.  
  147.    2.2449320656847304`*^8}, {44.6`, 2.2665178896336156`*^8}, {44.7`,
  148.  
  149.    2.288514384068139`*^8}, {44.8`, 2.3102158933570683`*^8}, {44.9`,
  150.  
  151.    2.3319905140780103`*^8}, {45.`, 2.3538982955529702`*^8}, {45.1`,
  152.  
  153.    2.3756453188722473`*^8}, {45.2`, 2.3976424959402618`*^8}, {45.3`,
  154.  
  155.    2.4196625445161235`*^8}, {45.4`, 2.4414164810844252`*^8}, {45.5`,
  156.  
  157.    2.4631146621622288`*^8}, {45.6`, 2.4850006332077453`*^8}, {45.7`,
  158.  
  159.    2.506777724299708`*^8}, {45.8`, 2.528298277117795`*^8}, {45.9`,
  160.  
  161.    2.549837638761892`*^8}, {46.`, 2.571549174245451`*^8}, {46.1`,
  162.  
  163.    2.5927650097765744`*^8}, {46.2`, 2.614159022212292`*^8}, {46.3`,
  164.  
  165.    2.6351491216254857`*^8}, {46.400000000000006`,
  166.  
  167.    2.656130720167915`*^8}, {46.5`, 2.676993607703125`*^8}, {46.6`,
  168.  
  169.    2.697453341583533`*^8}, {46.7`, 2.717991659860158`*^8}, {46.8`,
  170.  
  171.    2.738124497316633`*^8}, {46.900000000000006`,
  172.  
  173.    2.7580743364059085`*^8}, {47.`, 2.777887868370922`*^8}, {47.1`,
  174.  
  175.    2.797347506091462`*^8}, {47.2`, 2.816794268008462`*^8}, {47.3`,
  176.  
  177.    2.835709813301522`*^8}, {47.400000000000006`,
  178.  
  179.    2.85443491403207`*^8}, {47.5`, 2.8729090950057536`*^8}, {47.6`,
  180.  
  181.    2.8908443705797243`*^8}, {47.7`, 2.9086483592138124`*^8}, {47.8`,
  182.  
  183.    2.92616691674387`*^8}, {47.900000000000006`,
  184.  
  185.    2.9433509685529596`*^8}, {48.`, 2.960088737198468`*^8}, {48.1`,
  186.  
  187.    2.976565234506772`*^8}, {48.2`, 2.9925258891573817`*^8}, {48.3`,
  188.  
  189.    3.008277593335896`*^8}, {48.400000000000006`,
  190.  
  191.    3.023595999719454`*^8}, {48.5`, 3.0385356075930405`*^8}, {48.6`,
  192.  
  193.    3.053191558830993`*^8}, {48.7`, 3.0672176011885536`*^8}, {48.8`,
  194.  
  195.    3.0809176116870475`*^8}, {48.900000000000006`,
  196.  
  197.    3.0941555142050624`*^8}, {49.`, 3.106976561130896`*^8}, {49.1`,
  198.  
  199.    3.119210770647695`*^8}, {49.2`, 3.1312904585429984`*^8}, {49.3`,
  200.  
  201.    3.14269326294585`*^8}, {49.400000000000006`,
  202.  
  203.    3.1534993015989745`*^8}, {49.5`, 3.164288056676036`*^8}, {49.6`,
  204.  
  205.    3.174125897385491`*^8}, {49.7`, 3.1837495088897604`*^8}, {49.8`,
  206.  
  207.    3.193133049216981`*^8}, {49.900000000000006`,
  208.  
  209.    3.2019425708446777`*^8}, {50.`, 3.2101704306454825`*^8}, {50.1`,
  210.  
  211.    3.2178362134887046`*^8}, {50.2`, 3.224967487788363`*^8}, {50.3`,
  212.  
  213.    3.231829002710494`*^8}, {50.400000000000006`,
  214.  
  215.    3.2381834720160854`*^8}, {50.5`, 3.244295964255826`*^8}, {50.6`,
  216.  
  217.    3.249753718344813`*^8}, {50.7`, 3.254545514881043`*^8}, {50.8`,
  218.  
  219.    3.2591837400298154`*^8}, {50.900000000000006`,
  220.  
  221.    3.2632365933676875`*^8}, {51.`, 3.2667997897341466`*^8}, {51.1`,
  222.  
  223.    3.2703092651426727`*^8}, {51.2`, 3.272835357111885`*^8}, {51.3`,
  224.  
  225.    3.275149703003217`*^8}, {51.400000000000006`,
  226.  
  227.    3.277259862151018`*^8}, {51.5`, 3.27876579075336`*^8}, {51.6`,
  228.  
  229.    3.2799119733837414`*^8}, {51.7`, 3.280517264635551`*^8}, {51.8`,
  230.  
  231.    3.281051864182986`*^8}, {51.900000000000006`,
  232.  
  233.    3.281093417519217`*^8}, {52.`, 3.280684897143639`*^8}, {52.1`,
  234.  
  235.    3.279915117797926`*^8}, {52.2`, 3.2788317811779594`*^8}, {52.3`,
  236.  
  237.    3.277580039236555`*^8}, {52.400000000000006`,
  238.  
  239.    3.275948217595974`*^8}, {52.5`, 3.2741968660497904`*^8}, {52.6`,
  240.  
  241.    3.271747210017073`*^8}, {52.7`, 3.26958004513807`*^8}, {52.8`,
  242.  
  243.    3.266619178914003`*^8}, {52.900000000000006`,
  244.  
  245.    3.2638326351653737`*^8}, {53.`, 3.2604991073091185`*^8}, {53.1`,
  246.  
  247.    3.257298898482738`*^8}, {53.2`, 3.253993078430637`*^8}, {53.3`,
  248.  
  249.    3.250215389229545`*^8}, {53.400000000000006`,
  250.  
  251.    3.246508945613699`*^8}, {53.5`, 3.2425771317495406`*^8}, {53.6`,
  252.  
  253.    3.238621442780316`*^8}, {53.7`, 3.234660569237573`*^8}, {53.8`,
  254.  
  255.    3.2304347565793437`*^8}, {53.900000000000006`,
  256.  
  257.    3.226579982807768`*^8}, {54.`, 3.2225651841191244`*^8}, {54.1`,
  258.  
  259.    3.218416511829087`*^8}, {54.2`, 3.2144594671490043`*^8}, {54.3`,
  260.  
  261.    3.2103956410249674`*^8}, {54.400000000000006`,
  262.  
  263.    3.2067161343084925`*^8}, {54.5`, 3.2031151794254494`*^8}, {54.6`,
  264.  
  265.    3.1995676323550916`*^8}, {54.7`, 3.19626771998548`*^8}, {54.8`,
  266.  
  267.    3.193194423603016`*^8}, {54.900000000000006`,
  268.  
  269.    3.1905377917943543`*^8}, {55.`, 3.187621522016274`*^8}, {55.1`,
  270.  
  271.    3.185288615087083`*^8}, {55.2`, 3.183411334022303`*^8}, {55.3`,
  272.  
  273.    3.1814248803313965`*^8}, {55.400000000000006`,
  274.  
  275.    3.179967526419423`*^8}, {55.5`, 3.1788441370664203`*^8}, {55.6`,
  276.  
  277.    3.178214843787729`*^8}, {55.7`, 3.177761088569737`*^8}, {55.8`,
  278.  
  279.    3.1777585708835393`*^8}, {55.900000000000006`,
  280.  
  281.    3.178249325119825`*^8}, {56.`, 3.1789864177878904`*^8}, {56.1`,
  282.  
  283.    3.1800775050714767`*^8}, {56.2`, 3.181463200220724`*^8}, {56.3`,
  284.  
  285.    3.183489146761801`*^8}, {56.400000000000006`,
  286.  
  287.    3.1856665057856274`*^8}, {56.5`, 3.1880351280853844`*^8}, {56.6`,
  288.  
  289.    3.1907354706502706`*^8}, {56.7`, 3.194194364617434`*^8}, {56.8`,
  290.  
  291.    3.197451130817305`*^8}, {56.900000000000006`,
  292.  
  293.    3.2012117524567086`*^8}, {57.`, 3.204891793031598`*^8}, {57.1`,
  294.  
  295.    3.20900370470142`*^8}, {57.2`, 3.2133001073449385`*^8}, {57.3`,
  296.  
  297.    3.2178043207863295`*^8}, {57.400000000000006`,
  298.  
  299.    3.222354304574492`*^8}, {57.5`, 3.226777640787773`*^8}, {57.6`,
  300.  
  301.    3.231446544168732`*^8}, {57.7`, 3.235705329086409`*^8}, {57.8`,
  302.  
  303.    3.2402358090043026`*^8}, {57.900000000000006`,
  304.  
  305.    3.244647752492141`*^8}, {58.`, 3.2488508892401576`*^8}, {58.1`,
  306.  
  307.    3.252863569579624`*^8}, {58.2`, 3.2567307851998484`*^8}, {58.3`,
  308.  
  309.    3.260121485618781`*^8}, {58.400000000000006`,
  310.  
  311.    3.2637577520916504`*^8}, {58.5`, 3.266698262690229`*^8}, {58.6`,
  312.  
  313.    3.2693902222056985`*^8}, {58.7`, 3.271843402403688`*^8}, {58.8`,
  314.  
  315.    3.274008029045284`*^8}, {58.900000000000006`,
  316.  
  317.    3.275676434283537`*^8}, {59.`, 3.276976502336794`*^8}, {59.1`,
  318.  
  319.    3.278111028639613`*^8}, {59.2`, 3.278833145123387`*^8}, {59.3`,
  320.  
  321.    3.279203837939683`*^8}, {59.400000000000006`,
  322.  
  323.    3.2793486463030225`*^8}, {59.5`, 3.2793169009993476`*^8}, {59.6`,
  324.  
  325.    3.278800781384124`*^8}, {59.7`, 3.278085422278357`*^8}, {59.8`,
  326.  
  327.    3.2773228249207693`*^8}, {59.900000000000006`,
  328.  
  329.    3.276425686405634`*^8}, {60.`, 3.27525355089541`*^8}, {60.1`,
  330.  
  331.    3.274205735459841`*^8}, {60.2`, 3.273102998587829`*^8}, {60.3`,
  332.  
  333.    3.271910896669299`*^8}, {60.400000000000006`,
  334.  
  335.    3.270905043842246`*^8}, {60.5`, 3.2700377989076877`*^8}, {60.6`,
  336.  
  337.    3.269513295509556`*^8}, {60.7`, 3.2689788756182164`*^8}, {60.8`,
  338.  
  339.    3.269143593662879`*^8}, {60.900000000000006`,
  340.  
  341.    3.269166613731415`*^8}, {61.`, 3.26971477573183`*^8}, {61.1`,
  342.  
  343.    3.2710141498697007`*^8}, {61.2`, 3.272539281600619`*^8}, {61.3`,
  344.  
  345.    3.274463006551376`*^8}, {61.400000000000006`,
  346.  
  347.    3.2772409906367373`*^8}, {61.5`, 3.280420169744154`*^8}, {61.6`,
  348.  
  349.    3.283980366076435`*^8}, {61.7`, 3.288302842707272`*^8}, {61.8`,
  350.  
  351.    3.2934515797862315`*^8}, {61.900000000000006`,
  352.  
  353.    3.299096974195862`*^8}, {62.`, 3.305271441580234`*^8}, {62.1`,
  354.  
  355.    3.3122888418547016`*^8}, {62.2`,
  356.  
  357.    3.31988917254585`*^8}, {62.300000000000004`,
  358.  
  359.    3.3283874468914557`*^8}, {62.4`, 3.337092715026067`*^8}, {62.5`,
  360.  
  361.    3.3468427326762915`*^8}, {62.6`, 3.356743378003942`*^8}, {62.7`,
  362.  
  363.    3.367263448689508`*^8}, {62.800000000000004`,
  364.  
  365.    3.378299690685078`*^8}, {62.9`, 3.389910674491396`*^8}, {63.`,
  366.  
  367.    3.402016763027503`*^8}, {63.1`, 3.414245297564836`*^8}, {63.2`,
  368.  
  369.    3.426780800382867`*^8}, {63.300000000000004`,
  370.  
  371.    3.4396742322083825`*^8}, {63.4`, 3.452682035823037`*^8}, {63.5`,
  372.  
  373.    3.46605361211071`*^8}, {63.6`, 3.479388082514252`*^8}, {63.7`,
  374.  
  375.    3.4925874871254766`*^8}, {63.800000000000004`,
  376.  
  377.    3.506170420950347`*^8}, {63.9`, 3.519154600804556`*^8}, {64.`,
  378.  
  379.    3.5322981180072695`*^8}, {64.1`, 3.545289680868335`*^8}, {64.2`,
  380.  
  381.    3.5581905592917013`*^8}, {64.30000000000001`,
  382.  
  383.    3.5704211897827905`*^8}, {64.4`, 3.582551007439326`*^8}, {64.5`,
  384.  
  385.    3.5944058332260954`*^8}, {64.6`, 3.605756569776652`*^8}, {64.7`,
  386.  
  387.    3.6164205471404`*^8}, {64.80000000000001`,
  388.  
  389.    3.6267334190372324`*^8}, {64.9`, 3.636912133768186`*^8}, {65.`,
  390.  
  391.    3.6463633306688815`*^8}, {65.1`, 3.655268696122972`*^8}, {65.2`,
  392.  
  393.    3.6635543329660434`*^8}, {65.30000000000001`,
  394.  
  395.    3.6713083830029327`*^8}, {65.4`, 3.6786690749869096`*^8}, {65.5`,
  396.  
  397.    3.6855780815393907`*^8}, {65.6`, 3.691908229707714`*^8}, {65.7`,
  398.  
  399.    3.697965049602041`*^8}, {65.80000000000001`,
  400.  
  401.    3.703591469341711`*^8}, {65.9`, 3.708577605620679`*^8}, {66.`,
  402.  
  403.    3.713301150699523`*^8}, {66.1`, 3.7174966664643246`*^8}, {66.2`,
  404.  
  405.    3.72155464312756`*^8}, {66.30000000000001`,
  406.  
  407.    3.725194219931153`*^8}, {66.4`, 3.728704857414911`*^8}, {66.5`,
  408.  
  409.    3.7319707439847356`*^8}, {66.6`, 3.735130511447414`*^8}, {66.7`,
  410.  
  411.    3.738044033955617`*^8}, {66.80000000000001`,
  412.  
  413.    3.740787563024863`*^8}, {66.9`, 3.743872606550009`*^8}, {67.`,
  414.  
  415.    3.7466644826032037`*^8}, {67.1`, 3.7494796011228526`*^8}, {67.2`,
  416.  
  417.    3.752567597003231`*^8}, {67.30000000000001`,
  418.  
  419.    3.755528597048888`*^8}, {67.4`, 3.758867756212609`*^8}, {67.5`,
  420.  
  421.    3.762253642893481`*^8}, {67.6`, 3.765563102468596`*^8}, {67.7`,
  422.  
  423.    3.76929418561745`*^8}, {67.80000000000001`,
  424.  
  425.    3.7734675463116574`*^8}, {67.9`, 3.777793327652246`*^8}, {68.`,
  426.  
  427.    3.78229775978888`*^8}, {68.1`, 3.786886041650724`*^8}, {68.2`,
  428.  
  429.    3.7921789819393176`*^8}, {68.30000000000001`,
  430.  
  431.    3.797400500297583`*^8}, {68.4`, 3.803015063139415`*^8}, {68.5`,
  432.  
  433.    3.809133938161798`*^8}, {68.6`, 3.815470864869022`*^8}, {68.7`,
  434.  
  435.    3.8222110275067735`*^8}, {68.80000000000001`,
  436.  
  437.    3.829037876579778`*^8}, {68.9`, 3.835887650962112`*^8}, {69.`,
  438.  
  439.    3.8435872719788545`*^8}, {69.1`, 3.851099881198329`*^8}, {69.2`,
  440.  
  441.    3.859342022241043`*^8}, {69.30000000000001`,
  442.  
  443.    3.867606783579225`*^8}, {69.4`, 3.8759238864472085`*^8}, {69.5`,
  444.  
  445.    3.8845543983685386`*^8}, {69.6`, 3.89349742833892`*^8}, {69.7`,
  446.  
  447.    3.9029700599955183`*^8}, {69.80000000000001`,
  448.  
  449.    3.912164360108417`*^8}, {69.9`, 3.9216870483815044`*^8}, {70.`,
  450.  
  451.    3.931483602407363`*^8}, {70.1`, 3.9410969655725956`*^8}, {70.2`,
  452.  
  453.    3.9513894195117193`*^8}, {70.30000000000001`,
  454.  
  455.    3.9616302770693266`*^8}, {70.4`, 3.9716179276966923`*^8}, {70.5`,
  456.  
  457.    3.982128004055699`*^8}, {70.6`, 3.9927235371773595`*^8}, {70.7`,
  458.  
  459.    4.003509165551405`*^8}, {70.80000000000001`,
  460.  
  461.    4.014030954849812`*^8}, {70.9`, 4.024962310210849`*^8}, {71.`,
  462.  
  463.    4.0358209453883123`*^8}, {71.1`, 4.0467931066218626`*^8}, {71.2`,
  464.  
  465.    4.057559335818988`*^8}, {71.30000000000001`,
  466.  
  467.    4.068661360490459`*^8}, {71.4`, 4.080052928244571`*^8}, {71.5`,
  468.  
  469.    4.091169005766451`*^8}, {71.6`, 4.102270502761557`*^8}, {71.7`,
  470.  
  471.    4.1138578091937834`*^8}, {71.80000000000001`,
  472.  
  473.    4.1251155210459965`*^8}, {71.9`, 4.1364941770349723`*^8}, {72.`,
  474.  
  475.    4.1479014369837517`*^8}, {72.1`, 4.1593269639449173`*^8}, {72.2`,
  476.  
  477.    4.170742870689594`*^8}, {72.30000000000001`,
  478.  
  479.    4.1822452256960845`*^8}, {72.4`, 4.19377280532779`*^8}, {72.5`,
  480.  
  481.    4.2054366824108684`*^8}, {72.6`, 4.2166842641395384`*^8}, {72.7`,
  482.  
  483.    4.2283143004773027`*^8}, {72.80000000000001`,
  484.  
  485.    4.240186452268878`*^8}, {72.9`, 4.251858930902198`*^8}, {73.`,
  486.  
  487.    4.2632043093512297`*^8}, {73.1`, 4.2749669162602633`*^8}, {73.2`,
  488.  
  489.    4.286488051759319`*^8}, {73.30000000000001`,
  490.  
  491.    4.2984178745028996`*^8}, {73.4`, 4.310295223728404`*^8}, {73.5`,
  492.  
  493.    4.321735402830664`*^8}, {73.6`, 4.333782728263144`*^8}, {73.7`,
  494.  
  495.    4.345358988467979`*^8}, {73.80000000000001`,
  496.  
  497.    4.357223893549589`*^8}, {73.9`, 4.368848682802905`*^8}, {74.`,
  498.  
  499.    4.381027789385323`*^8}, {74.1`, 4.392741358383033`*^8}, {74.2`,
  500.  
  501.    4.404523585717674`*^8}, {74.30000000000001`,
  502.  
  503.    4.4165790947883683`*^8}, {74.4`, 4.428364190545716`*^8}, {74.5`,
  504.  
  505.    4.4402239358321327`*^8}, {74.6`, 4.4520446091178155`*^8}, {74.7`,
  506.  
  507.    4.4640519786811954`*^8}, {74.80000000000001`,
  508.  
  509.    4.476036989583306`*^8}, {74.9`, 4.488212002259784`*^8}, {75.`,
  510.  
  511.    4.50016702284048`*^8}, {75.1`, 4.511972553212277`*^8}, {75.2`,
  512.  
  513.    4.524220942450661`*^8}, {75.30000000000001`,
  514.  
  515.    4.535864223500876`*^8}, {75.4`, 4.548080748479822`*^8}, {75.5`,
  516.  
  517.    4.5601746125887114`*^8}, {75.6`, 4.572276632620131`*^8}, {75.7`,
  518.  
  519.    4.5842204294312906`*^8}, {75.80000000000001`,
  520.  
  521.    4.596595601715562`*^8}, {75.9`, 4.6086719915453845`*^8}, {76.`,
  522.  
  523.    4.621003597526813`*^8}, {76.1`, 4.633181225397551`*^8}, {76.2`,
  524.  
  525.    4.645075736797058`*^8}, {76.30000000000001`,
  526.  
  527.    4.6572917749778366`*^8}, {76.4`, 4.669677776682115`*^8}, {76.5`,
  528.  
  529.    4.6817393303444105`*^8}, {76.6`, 4.694097331758887`*^8}, {76.7`,
  530.  
  531.    4.7063634216148126`*^8}, {76.80000000000001`,
  532.  
  533.    4.718550648902903`*^8}, {76.9`, 4.730765984916197`*^8}, {77.`,
  534.  
  535.    4.743005877022079`*^8}, {77.1`, 4.7555780128627634`*^8}, {77.2`,
  536.  
  537.    4.767942632996616`*^8}, {77.30000000000001`,
  538.  
  539.    4.7802375760961807`*^8}, {77.4`, 4.7925397668717337`*^8}, {77.5`,
  540.  
  541.    4.8048861965555114`*^8}, {77.6`, 4.8175653709807914`*^8}, {77.7`,
  542.  
  543.    4.829827580688622`*^8}, {77.80000000000001`,
  544.  
  545.    4.8422945178250486`*^8}, {77.9`, 4.8547249122644436`*^8}, {78.`,
  546.  
  547.    4.8672796681303275`*^8}, {78.1`, 4.8798126343066156`*^8}, {78.2`,
  548.  
  549.    4.892411314706424`*^8}, {78.30000000000001`,
  550.  
  551.    4.9046098222035885`*^8}, {78.4`, 4.9171998456332064`*^8}, {78.5`,
  552.  
  553.    4.92980620158521`*^8}, {78.6`, 4.9423632445393294`*^8}, {78.7`,
  554.  
  555.    4.954999596248343`*^8}, {78.80000000000001`,
  556.  
  557.    4.967454296502248`*^8}, {78.9`, 4.97999207575818`*^8}, {79.`,
  558.  
  559.    4.992957422144184`*^8}, {79.1`, 5.005565145880107`*^8}, {79.2`,
  560.  
  561.    5.0182334144736385`*^8}, {79.30000000000001`,
  562.  
  563.    5.0308015555361265`*^8}, {79.4`, 5.043309905755629`*^8}, {79.5`,
  564.  
  565.    5.0560503568740994`*^8}, {79.6`, 5.068821803821288`*^8}, {79.7`,
  566.  
  567.    5.0815950435918`*^8}, {79.80000000000001`,
  568.  
  569.    5.094525165929398`*^8}, {79.9`, 5.1070616158582973`*^8}, {80.`,
  570.  
  571.    5.12012990408671`*^8}}
  572.  
  573.  
  574.  
  575.  
  576. peakfunc[A_, \[Mu]_, \[Sigma]_, x_] =
  577.  
  578.  A^2 E^(-((x - \[Mu])^2/(2 \[Sigma]^2)))
  579.  
  580. Clear[model]
  581.  
  582.  
  583.  
  584.  
  585. model[data_, n_] :=
  586.  
  587.  Module[{dataconfig, modelfunc, objfunc, fitvar, fitres},
  588.  
  589.   dataconfig = {A[#], \[Mu][#], \[Sigma][#]} & /@ Range[n];
  590.  
  591.   modelfunc = peakfunc[##, fitvar] & @@@ dataconfig // Total;
  592.  
  593.   objfunc =
  594.  
  595.    Total[((Sqrt[data[[All, 2]]])/
  596.  
  597.        data[[All,
  598.  
  599.         1]]) (data[[All, 2]] - (modelfunc /. fitvar -> # &) /@
  600.  
  601.          data[[All, 1]])^2];
  602.  
  603.   FindMinimum[objfunc, Flatten@dataconfig]]
  604.  
  605.  
  606.  
  607.  
  608. Clear[modelvalue]
  609.  
  610. modelvalue[data_, n_] /; NumericQ[n] :=
  611.  
  612.  If[n >= 1, model[data, n][[1]], 0]
  613.  
  614.  
  615.  
  616.  
  617. fitres = ReleaseHold[
  618.  
  619.    Hold[{Round[n], model[data, Round[n]]}] /.
  620.  
  621.     FindMinimum[modelvalue[data, Round[n]], {n, 3},
  622.  
  623.       Method -> "PrincipalAxis"][[2]]] // Quiet
  624.  
  625.  
  626.  
  627.  
  628. resfunc =
  629.  
  630.   peakfunc[A[#], \[Mu][#], \[Sigma][#], x] & /@ Range[fitres[[1]]] /.
  631.  
  632.    fitres[[2, 2]];
  633.  
  634.  
  635.  
  636.  
  637. Show@{Plot[Evaluate[resfunc], {x, 30, 80},
  638.  
  639.    PlotStyle -> ({Directive[Dashed, Thick,
  640.  
  641.          ColorData["Rainbow"][#]]} & /@
  642.  
  643.       Rescale[Range[Length[resfunc]]]), PlotRange -> All,
  644.  
  645.    Frame -> True, Axes -> False],
  646.  
  647.   Plot[Evaluate[Total@resfunc], {x, 30, 80},
  648.  
  649.    PlotStyle -> Directive[Thick, Red], PlotRange -> All,
  650.  
  651.    Frame -> True, Axes -> False],
  652.  
  653.   Graphics[{PointSize[.003], Gray, Line@data}]}
  654.  
  655.  
  656.  
  657.  
  658. {#, modelvalue[data, #]} & /@ Range[1, 13] //
  659.  
  660.   ListLogPlot[#, Joined -> True] & // Quiet
  661.  
  662.  
  663.  
  664.  
  665. dataflat = Flatten@data;
  666.  
  667. (*Residuals plot*)
  668.  
  669. ListPlot[Transpose[{dataflat[[
  670.  
  671.     1 ;; -1 ;;
  672.  
  673.      2]], (dataflat[[2 ;; -1 ;; 2]] -
  674.  
  675.       Table[Evaluate[Total@resfunc], {x, 30, 80, 0.1}])/
  676.  
  677.     Sqrt[data[[All, 2]]]}], Filling -> Axis]
RAW Paste Data
We use cookies for various purposes including analytics. By continuing to use Pastebin, you agree to our use of cookies as described in the Cookies Policy. OK, I Understand
 
Top