4Misc_Start4Platform@ ROGIXXz;,P VT$mL,winspoolM9-107 SHARP MX-4141N PS132.246.76.222M9-107 SHARP MX-4141N PS߀ odXXLetterSCDM$,SS0HM9-107 SHARP MX-4141N PS߀ odXXLetter4I #d$dd****d ܑddo ddL o222222< $ ROGIXXz;,P VT$mL,winspoolM9-107 SHARP MX-4141N PS132.246.76.222M9-107 SHARP MX-4141N PS߀ odXXLetterSCDM$,SS0HM9-107 SHARP MX-4141N PS߀ odXXLetter4I #d$dd****d ܑddo ddL o222222< $ ROGIXXp 1< VT$mL,winspoolM9-107 SHARP MX-4141N PS132.246.76.222M9-107 SHARP MX-4141N PS߀ odXXLetterSCDM$,SS0HM9-107 SHARP MX-4141N PS߀ odXXLetter4I #d$dd****d ܑddo ddL o222222< $  ROGIXXz;,P VT$mL,winspoolM9-107 SHARP MX-4141N PS132.246.76.222M9-107 SHARP MX-4141N PS߀ odXXLetterSCDM$,SS0HM9-107 SHARP MX-4141N PS߀ odXXLetter4I #d$dd****d ܑddo ddL o222222< $^Graph*@@wwwwww?wwwwww? xWDashSettings#  !>_Normal@ Arial<HHHH$$>_Normal@ Arial<HHHH$$444444 Normal@ Consolas<HHHH$$4 4 4 4 4 4 shome?7d"^RC:Users:corbinj:Desktop:C:Users:corbinj:Desktop,hWOؿOP(6v,h/PVO`WO(*wl `WO<(*w,hWOL*whWOntlhnt1$Z\! h/ؿO\Hu4.hntcc 5RecentWindowsallDetailsfullhistoryGraph0:yw vs xwGraph2:yw vs xw;Peak0_without vs xw;...Igor Reference.ihfTable1:woC.ldUsing Igor.ihf 4Misc_EndXOPState_StartData Browseresults_with[*][0] GizmoPeakFunctions2Table1:woC.ldScatNLaywoC.ldUsing Igor.ihfaData Browserroot4XOPState_EndV_Flagmake_and_fit_test_data() withC[0]= {176.417,136.908,8.03353,473.063,126.224,8.90722,-7.67906,131.253,-3.21967} woC[0]= {156.447,132.204,10.0537,462.703,129.075,9.38499,-38.229,118.472,-6.58927} make_and_fit_test_data() withC[0]= {124.071,140,7.07106,354.492,130,7.07107,177.245,120,7.07107} woC[0]= {125.505,137.476,5.95915,-1.37866e+10,-9.70614e+06,1.87513e+07,8.30368e+11,4.3072e+09,-4.74528e+15} SavePICT/O/P=home/E=-5/B=144 as "fit_with_or_without_constraint.png" !(N@ MՀxw(N????BBBBBBBBBBBBBBBBBBBBBBBBBBBBCCCCCCCCC C C C C CCCCCCCCCCCCCCCCCCC C!C"C#C$C%C&C'C(C)C*C+C,C-C.C/C0C1C2C3C4C5C6C7C8C9C:C;CC?C@CACBCCCDCECFCGCHCICJCKCLCMCNCOCPCQCRCSCTCUCVCWCXCYCZC[C\C]C^C_C`CaCbCcC@MՀyw(N????>><Ɍ>`>I?K?uۋ?_??S!@WoN@x@@@@~AVAG5AMAmfAF1AAmARlA;?A AA A|A{-AAA-A{5A2LAzA%A[AAMAȀAhA@PA8A6"A AN@ @@M@j@>U<@1@6?B?I?F?x?->e>:xF>>=6e==B<^<~<@;hL5;:K c:9K9 968V 8@7%;7ʹ{6y5S5\145$4H3u2|A2l1WJ0=0K׏/".*.c^-d,+(+fJ=*t) (J'&Ս &)#%:$ R#Gg"{y! G*Ifm"9s\_.I/12V..<^D+] 0 eTC $Kd`M cw(N ????M.]Ҙd0/.}' 1L6@M:tY}cڪlVXt6z>`8?ŵC`vڋ(fR5ڌ)w!&*-G01-344155556:6P6a6p6}666666666666666777"7*737<7D7M7V7^77F Mconstraints(N????MK0 > 0[@M W_sigma(N ????G@ua ݂?^OR? g?ץ.Xd?]q3?-.uWE?\9??~JPנ?$dM withC(N????M-EBCC{pBCmB A=AcAxlocsigma$WUdM woC(N????@M.0C]T>>VͿ>>$?T??ר??Vy?f%@Ч7@ X@Xz@ڏ@Tܢ@! @Y@@z@@"AA A Af AN AU^A@3@i@@l@3S@Y/@w@6#g@ÙE@o&@ @?̹??l?8?nm?Z>VY>r>\2>,>X=!=~d5=<;)|M)T(t-(F']&U% (+%@|$[2#y#zR:"@Y!# [%[J~G+%,<*h#m#U`.A/>19>OL|>>>?VC?|?ġ?V?D@@UC@l@Z7@@@{@RAAl+AY?A[SAgAUyAAA>zAAAA^AwAAڙASIAmAzARTyAnfAdRA>Av*A cA-AYq@W@/@p@ik@A@@2?QW?o?z?B? ?>>y>j7>0/>#g==?=]=ʷ!y7k[cܽ<9-京>6 3t۷(Z/ǵ)ǃBHѡıEKM׽ 3cT ժzy è|,? Ģ7m+@՞:{".UUs/:}KFeޏ54'ъgNPn氃Vف.@HO9Tw@MՀplottedResults_with(N????"Q=ra=1j>>b?%?У=@C:@b! A\AjANB49BrBBBw C'C&BCsYCUjCfMtC5tClC:\CWFCa+CqC fBBBwBBjBصAskAA@ZM@?w??n>g=Ri= <3!6<G;:r>:t 9+8B8ߓA7z65e4,3r2/2 1. 0/@-3,h+2|* ))M('?%'+$7P#@" K8c+cYFztnE1sN c' r 98jo,L MiK566 %7RL8959nQ:=:?;c4#q>g+??@`h@:@pA?nAԳA@B9B~BBBnC6/&CCC_CvCMC8C)ÊCpCJCXlC*RC-Q6CCAjBB:B]BNBAAwDAx@V@ܙ7@?n??I># >=T=z7Ԇ654ԣ4dF3p2&q1C0I/.֪-=,+ v*l)ٗ(42'[&2%$ $"f !d $OФ2 G4;9}wHa? !  wFVkZ9!J}Aqp Ř13¬4j%5*6%'7891r9R:;P:<<p=<>~>e??6Va@_@1AoVA8A A%BeB~B"RB@B cBpC!CBrBzKB/B,`kB,"*BNAA]Y9As|@=m@X?5t?PT>*J>큩=X=+gJ<5;*]:a999+8hr:7\~?6?:5‹+4L34-1V05/H.#u-+ H*o(b!q'=%}i$z"m=6! IJ&bsf̖*ZfLe j ^d>^_60f4d Mն ew(N ???? .........*fdTXET????fullhistoryfullhistoryfullhistorymake xw= p+100 make yw display yw vs xw yw=0 yw = 20*exp(-(xw - 130)^2 / 10^2) yw += 10*exp(-(xw - 120)^2 / 10^2) yw += 7*exp(-(xw - 140)^2 / 10^2) make_and_fit_test_data() cw[0]= {-115.577,133.655,6.32857,851.153,129.399,9.27524,79.0518,125.122,-5.79601} cw[0]= {248.532,121.623,7.42073,269.898,135.961,8.00571,137.698,129.595,6.084} make_and_fit_test_data() cw[0]= {188.634,121.023,7.56074,326.631,133.762,8.99091,141.908,129.634,7.10228} cw[0]= {93.039,129.572,5.79725,312.49,135.017,8.24895,-250.56,121.786,-7.4843} make_and_fit_test_data() cw[0]= {254.417,121.75,7.45479,260.068,136.199,7.9405,141.498,129.696,6.09703} cw[0]= {8.71804e+07,1.16016e+09,3.58728e+08,58982.6,2736.91,320.379,-58337.9,7403.28,-722.828} make_and_fit_test_data() cw[0]= {221.756,133.109,9.47495,265.378,126.746,9.23311,172.037,127.762,10.1445} cw[0]= {31.3552,130.442,4.91898,835.24,129.334,9.32575,208.672,129.885,-7.5671} make_and_fit_test_data() cw[0]= {33.538,141.785,0.0461877,29.5014,129.74,8.27743,11.5679,132.982,8.30045} cw[0]= {-33.6093,124.68,5.79088,504.419,130.785,9.7205,-186.984,124.161,-8.18099} make_and_fit_test_data() cw[0]= {259.444,121.854,7.48327,244.938,136.573,7.84098,151.553,129.817,6.14093} cw[0]= {-3201.79,164.784,7.8527,3656.78,202.703,29.4069,2.06524e+36,-4.30404e+33,-2.18059e+08} make_and_fit_test_data() cw[0]= {205.499,121.542,7.73573,303.692,133.958,9.01567,148.078,130.109,7.28369} cw[0]= {234.962,121.285,7.33045,241.168,136.629,7.84401,179.77,129.561,6.29081} make_and_fit_test_data() cw[0]= {196.104,121.912,7.94748,296.002,133.134,9.37239,165.32,130.802,7.86748} cw[0]= {47.747,129.077,5.35237,432.621,132.701,8.89102,-175.992,120.701,-7.27783} make_and_fit_test_data() cw[0]= {-0.548548,139.566,-0.054616,434.294,129.636,10.1267,225.213,128.129,9.65482} cw[0]= {190.26,121.652,7.83806,326.875,133.041,9.31372,140.259,130.516,7.58283} make_and_fit_test_data() cw[0]= {-0.548548,139.566,-0.054616,434.294,129.636,10.1267,225.213,128.129,9.65482} cw[0]= {190.26,121.652,7.83806,326.875,133.041,9.31372,140.259,130.516,7.58283} make_and_fit_test_data() cw[0]= {-0.548548,139.566,-0.054616,434.294,129.636,10.1267,225.213,128.129,9.65482} cw[0]= {190.26,121.652,7.83806,326.875,133.041,9.31372,140.259,130.516,7.58283} make_and_fit_test_data() cw[0]= {-0.548548,139.566,-0.054616,434.294,129.636,10.1267,225.213,128.129,9.65482} cw[0]= {190.26,121.652,7.83806,326.875,133.041,9.31372,140.259,130.516,7.58283} make_and_fit_test_data() cw[0]= {-0.548548,139.566,-0.054616,434.294,129.636,10.1267,225.213,128.129,9.65482} cw[0]= {190.26,121.652,7.83806,326.875,133.041,9.31372,140.259,130.516,7.58283} make_and_fit_test_data() cw[0]= {-0.548548,139.566,-0.054616,434.294,129.636,10.1267,225.213,128.129,9.65482} cw[0]= {190.26,121.652,7.83806,326.875,133.041,9.31372,140.259,130.516,7.58283} make_and_fit_test_data() cw[0]= {-0.548548,139.566,-0.054616,434.294,129.636,10.1267,225.213,128.129,9.65482} cw[0]= {190.26,121.652,7.83806,326.875,133.041,9.31372,140.259,130.516,7.58283} make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() edit woc.ld make_and_fit_test_data() append plottedResults_without[][0] /TN=Peak0 append plottedResults_without[][0] /TN=Peak1 append plottedResults_without[][1] /TN=Peak1 remove peak1 append plottedResults_without[][2] /TN=Peak2 vs xw remove peak0 remove peak1 RemoveFromGraph Peak2 append plottedResults_without[][0] /TN=Peak0 vs xw display yw vs xw appendtograph plottedResults_without[][0] /TN=Peak0 vs xw appendtograph plottedResults_without[][1] /TN=Peak1 vs xw appendtograph plottedResults_without[][2] /TN=Peak2 vs xw appendtograph plottedResults_with[][0] /TN=Peak0 vs xw appendtograph plottedResults_with[][1] /TN=Peak1 vs xw appendtograph plottedResults_with[][2] /TN=Peak2 vs xw display yw vs xw appendtograph plottedResults_without[][0] /TN=Peak0_without vs xw appendtograph plottedResults_without[][1] /TN=Peak1_without vs xw appendtograph plottedResults_without[][2] /TN=Peak2_without vs xw appendtograph plottedResults_with[][0] /TN=Peak0_with vs xw appendtograph plottedResults_with[][1] /TN=Peak1_with vs xw appendtograph plottedResults_with[][2] /TN=Peak2_with vs xw ModifyGraph lsize(yw)=3,rgb(yw)=(0,0,0),rgb(Peak0_with)=(0,15872,65280);DelayUpdate ModifyGraph rgb(Peak1_with)=(0,15872,65280),rgb(Peak2_with)=(0,15872,65280) make /o xw= p+100 make /o yw= 0 yw += 20*exp(-(xw - 130)^2 / 10^2) yw += 10*exp(-(xw - 120)^2 / 10^2) yw += 7*exp(-(xw - 140)^2 / 10^2) make_and_fit_test_data() make_and_fit_test_data() legend Legend/C/N=text0/J "\\s(yw) yw\r\\s(Peak0_with) Peak0_with\r\\s(Peak1_with) Peak1_with\r\\s(Peak2_with) Peak2_with\r\\s(Peak0_without) Peak0_without";DelayUpdate AppendText/N=text0 "\\s(Peak1_without) Peak1_without\r\\s(Peak2_without) Peak2_without\r...constraints" make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() Legend/C/N=text0/J "\\s(yw) yw\r\\s(Peak0_with) Peak0_with constraints\r\\s(Peak1_with) Peak1_with\r\\s(Peak2_with) Peak2_with\r\\s(Peak0_without) Peak0_without";DelayUpdate AppendText/N=text0 "\\s(Peak1_without) Peak1_without\r\\s(Peak2_without) Peak2_without" make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() 0 0 make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() make_and_fit_test_data() withC[0]= {176.417,136.908,8.03353,473.063,126.224,8.90722,-7.67906,131.253,-3.21967} woC[0]= {156.447,132.204,10.0537,462.703,129.075,9.38499,-38.229,118.472,-6.58927} d TXET????GallDetailsallDetails//#pragma rtGlobals=3 // Use modern global access method and strict wave access. // //// A 'prototype' function for when functions take functions as arguments. //function jFit_proto(coeffwave,x) // // For use with FUNCREF // wave coeffwave // variable x //end // //function GaussSum(w, x) //wave w //variable x // // variable out, j, nGauss= 3 // for (j=0; j 0"} // simple constraint wave for tests // // // ??? differnt results here! // if (withConstraint) // FuncFit /NTHR=0 /N/Q/W=2 GaussSum, cw, ydata /X=xdata /D=fitResult /C=constraints // else // FuncFit /NTHR=0 /N/Q/W=2 GaussSum, cw, ydata /X=xdata /D=fitResult // endif // // if (1) // add epsilon wave // duplicate /o cw, ew // ew= 1e-10 // // copy-paste of above lines // if (withConstraint) // FuncFit /NTHR=0 /N/Q/W=2 GaussSum, cw, ydata /X=xdata /D=fitResult /C=constraints /E=ew // else // FuncFit /NTHR=0 /N/Q/W=2 GaussSum, cw, ydata /X=xdata /D=fitResult /E=ew // endif // endif // // // print V_FitError // return cw //end // // //function make_and_fit_test_data() // setrandomseed 0.2 // // make /o xw= p+100 // make /o yw= 0 // yw += 20*exp(-(xw - 130)^2 / 10^2) // yw += 10*exp(-(xw - 120)^2 / 10^2) // yw += 7*exp(-(xw - 140)^2 / 10^2) // // // these give different results...??? only difference is at "???" above // wave cw= fitGaussSum(ydata=yw, xdata=xw, withConstraint=0) // duplicate /o cw, withC // wave cw= fitGaussSum(ydata=yw, xdata=xw, withConstraint=1) // duplicate /o cw, woC // // // variable plotTheResults= 1 // if (plotTheResults) // // undo flatten // Redimension /N=(3, 3) withC; matrixop /o withC= withC^t // SetDimLabel 1, 0, A, withC // SetDimLabel 1, 1, xloc, withC // SetDimLabel 1, 2, sigma, withC // // Redimension /N=(3, 3) woC; matrixop /o woC= woC^t // SetDimLabel 1, 0, A, woC // SetDimLabel 1, 1, xloc, woC // SetDimLabel 1, 2, sigma, woC // // make /o/n=(numpnts(yw), 3) plottedResults_without // plottedResults_without[][0]= woC[0][%A] / woC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - woC[0][%xloc])^2 / (2 * woC[0][%sigma]^2)) // plottedResults_without[][1]= woC[1][%A] / woC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - woC[1][%xloc])^2 / (2 * woC[1][%sigma]^2)) // plottedResults_without[][2]= woC[2][%A] / woC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - woC[2][%xloc])^2 / (2 * woC[2][%sigma]^2)) // // make /o/n=(numpnts(yw), 3) plottedResults_with // plottedResults_without[][0]= withC[0][%A] / withC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - withC[0][%xloc])^2 / (2 * withC[0][%sigma]^2)) // plottedResults_without[][1]= withC[1][%A] / withC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - withC[1][%xloc])^2 / (2 * withC[1][%sigma]^2)) // plottedResults_without[][2]= withC[2][%A] / withC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - withC[2][%xloc])^2 / (2 * withC[2][%sigma]^2)) // // display /K=1 yw vs xw // appendtograph plottedResults_without[][0] /TN=Peak0_without vs xw // appendtograph plottedResults_without[][1] /TN=Peak1_without vs xw // appendtograph plottedResults_without[][2] /TN=Peak2_without vs xw // appendtograph plottedResults_with[][0] /TN=Peak0_with vs xw // appendtograph plottedResults_with[][1] /TN=Peak1_with vs xw // appendtograph plottedResults_with[][2] /TN=Peak2_with vs xw // endif //end  // Platform=WindowsNT, IGORVersion=6.372, architecture=Intel, systemTextEncoding="Windows-1252", historyTextEncoding="Windows-1252", procwinTextEncoding="Windows-1252" Silent 101 // use | as bitwise or -- not comment. DefaultFont "Arial" String/G root:gWMSetNextTextFilesTextEncoding = "Windows-1252" // Text encoding for "allDetails". Used by Igor Pro 7. OpenProc/W=(242.25,215,741.75,509.75)/J=24614 "allDetails"; MoveWindow/P=$"allDetails" 0, 0, 0, 0 // Minimize the window. String/G root:gWMSetNextTextFilesTextEncoding = "Windows-1252" // Text encoding for fullhistory. Used by Igor Pro 7. OpenNotebook/N=fullhistory/W=(5.25,42.5,504.75,337.25)/J=18360 "fullhistory"; MoveWindow 0, 0, 0, 0 // Minimize the window. MoveWindow/P 731.25,189.5,1230.75,484.25 MoveWindow/C 6.75,574.25,1374.75,716.75 Table1() Graph0() Graph2() KillStrings/Z root:gWMSetNextTextFilesTextEncoding Window Graph2() : Graph PauseUpdate; Silent 1 // building window... Display /W=(180,215,635.25,516.5) yw vs xw AppendToGraph plottedResults_without[*][0]/TN=Peak0_without vs xw AppendToGraph plottedResults_without[*][1]/TN=Peak1_without vs xw AppendToGraph plottedResults_without[*][2]/TN=Peak2_without vs xw AppendToGraph plottedResults_with[*][0]/TN=Peak0_with vs xw AppendToGraph plottedResults_with[*][1]/TN=Peak1_with vs xw AppendToGraph plottedResults_with[*][2]/TN=Peak2_with vs xw ModifyGraph lSize(yw)=3 ModifyGraph rgb(yw)=(0,0,0),rgb(Peak0_with)=(0,15872,65280),rgb(Peak1_with)=(0,15872,65280) ModifyGraph rgb(Peak2_with)=(0,15872,65280) Legend/C/N=text0/J "\\s(yw) yw\r\\s(Peak0_with) Peak0_with constraints\r\\s(Peak1_with) Peak1_with\r\\s(Peak2_with) Peak2_with\r\\s(Peak0_without) Peak0_without" AppendText "\\s(Peak1_without) Peak1_without\r\\s(Peak2_without) Peak2_without" EndMacro Window Graph0() : Graph PauseUpdate; Silent 1 // building window... Display /W=(304.5,161.75,699,370.25) yw vs xw EndMacro Window Table1() : Table PauseUpdate; Silent 1 // building window... Edit/W=(5.25,42.5,510,236.75) woC.ld ModifyTable format(Point)=1 EndMacro {#pragma rtGlobals=3 // Use modern global access method and strict wave access. // A 'prototype' function for when functions take functions as arguments. function jFit_proto(coeffwave,x) // For use with FUNCREF wave coeffwave variable x end function GaussSum(w, x) wave w variable x variable out, j, nGauss= 3 for (j=0; j 0"} // simple constraint wave for tests // ??? differnt results here! if (withConstraint) FuncFit /NTHR=0 /N/Q/W=2 GaussSum, cw, ydata /X=xdata /D=fitResult /C=constraints else FuncFit /NTHR=0 /N/Q/W=2 GaussSum, cw, ydata /X=xdata /D=fitResult endif return cw end function make_and_fit_test_data() setrandomseed 0.3 make /o xw= p+100 make /o yw= 0 variable sigma_for_all= 10 yw += 20*exp(-(xw - 130)^2 / sigma_for_all^2) // NB: here without /sqrt(2pi) yw += 10*exp(-(xw - 120)^2 / sigma_for_all^2) yw += 7*exp(-(xw - 140)^2 / sigma_for_all^2) // these give different results...??? only difference is at "???" above wave cw= fitGaussSum(ydata=yw, xdata=xw, withConstraint=0) duplicate /o cw, withC wave cw= fitGaussSum(ydata=yw, xdata=xw, withConstraint=1) duplicate /o cw, woC print withC print woC variable plotTheResults= 1 if (plotTheResults) // undo flatten Redimension /N=(3, 3) withC; matrixop /o withC= withC^t SetDimLabel 1, 0, A, withC SetDimLabel 1, 1, xloc, withC SetDimLabel 1, 2, sigma, withC Redimension /N=(3, 3) woC; matrixop /o woC= woC^t SetDimLabel 1, 0, A, woC SetDimLabel 1, 1, xloc, woC SetDimLabel 1, 2, sigma, woC make /o/n=(numpnts(yw), 3) plottedResults_without plottedResults_without[][0]= woC[0][%A] / woC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - woC[0][%xloc])^2 / (2 * woC[0][%sigma]^2)) plottedResults_without[][1]= woC[1][%A] / woC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - woC[1][%xloc])^2 / (2 * woC[1][%sigma]^2)) plottedResults_without[][2]= woC[2][%A] / woC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - woC[2][%xloc])^2 / (2 * woC[2][%sigma]^2)) make /o/n=(numpnts(yw), 3) plottedResults_with plottedResults_without[][0]= withC[0][%A] / withC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - withC[0][%xloc])^2 / (2 * withC[0][%sigma]^2)) plottedResults_without[][1]= withC[1][%A] / withC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - withC[1][%xloc])^2 / (2 * withC[1][%sigma]^2)) plottedResults_without[][2]= withC[2][%A] / withC[0][%sigma] / sqrt(2*pi) * exp(-(xw[p] - withC[2][%xloc])^2 / (2 * withC[2][%sigma]^2)) display /K=1 yw vs xw appendtograph plottedResults_without[][0] /TN=Peak0_without vs xw appendtograph plottedResults_without[][1] /TN=Peak1_without vs xw appendtograph plottedResults_without[][2] /TN=Peak2_without vs xw appendtograph plottedResults_with[][0] /TN=Peak0_with vs xw appendtograph plottedResults_with[][1] /TN=Peak1_with vs xw appendtograph plottedResults_with[][2] /TN=Peak2_with vs xw endif end