(*^ ::[ frontEndVersion = "Microsoft Windows Mathematica Notebook Front End Version 2.2"; microsoftWindowsStandardFontEncoding; fontset = title, "Arial", 24, L0, center, nohscroll, bold; fontset = subtitle, "Arial", 18, L0, center, nohscroll, bold; fontset = subsubtitle, "Arial", 14, L0, center, nohscroll, bold; fontset = section, "Arial", 14, L0, bold, grayBox; fontset = subsection, "Arial", 12, L0, bold, blackBox; fontset = subsubsection, "Arial", 10, L0, bold, whiteBox; fontset = text, "Arial", 12, L0; fontset = smalltext, "Arial", 10, L0; fontset = input, "Courier New", 10, L0, nowordwrap; fontset = output, "Courier New", 10, L0, nowordwrap, bold; fontset = message, "Courier New", 10, L0, nowordwrap, R65280; fontset = print, "Courier New", 10, L0, nowordwrap; fontset = info, "Courier New", 10, L0, nowordwrap; fontset = postscript, "Courier New", 8, L0, nowordwrap; fontset = name, "Arial", 10, L0, nohscroll, italic, B65280; fontset = header, "Times New Roman", 10, L0, right, nohscroll; fontset = footer, "Times New Roman", 10, L0, right, nohscroll; fontset = help, "Arial", 10, L0, nohscroll; fontset = clipboard, "Arial", 12, L0, nohscroll; fontset = completions, "Arial", 12, L0, nowordwrap, nohscroll; fontset = graphics, "Courier New", 10, L0, nowordwrap, nohscroll; fontset = special1, "Arial", 12, L0, nowordwrap, nohscroll; fontset = special2, "Arial", 12, L0, center, nowordwrap, nohscroll; fontset = special3, "Arial", 12, L0, right, nowordwrap, nohscroll; fontset = special4, "Arial", 12, L0, nowordwrap, nohscroll; fontset = special5, "Arial", 12, L0, nowordwrap, nohscroll; fontset = leftheader, "Arial", 12, L0, nowordwrap, nohscroll; fontset = leftfooter, "Arial", 12, L0, nowordwrap, nohscroll; fontset = reserved1, "Courier New", 10, L0, nowordwrap, nohscroll;] :[font = input; nowordwrap; ] <All,PlotStyle->{Red,PointSize[0.02]}] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 299; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0845207 0.0535705 0.608848 0.148629 [ [(2.5)] .04941 .60885 0 2 Msboxa [(7.5)] .31726 .60885 0 2 Msboxa [(10)] .45118 .60885 0 2 Msboxa [(12.5)] .58511 .60885 0 2 Msboxa [(15)] .71904 .60885 0 2 Msboxa [(17.5)] .85296 .60885 0 2 Msboxa [(20)] .98689 .60885 0 2 Msboxa [(-4)] .17083 .01433 1 0 Msboxa [(-3)] .17083 .16296 1 0 Msboxa [(-2)] .17083 .31159 1 0 Msboxa [(-1)] .17083 .46022 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .04941 .60885 m .04941 .6151 L s P [(2.5)] .04941 .60885 0 2 Mshowa p .002 w .31726 .60885 m .31726 .6151 L s P [(7.5)] .31726 .60885 0 2 Mshowa p .002 w .45118 .60885 m .45118 .6151 L s P [(10)] .45118 .60885 0 2 Mshowa p .002 w .58511 .60885 m .58511 .6151 L s P [(12.5)] .58511 .60885 0 2 Mshowa p .002 w .71904 .60885 m .71904 .6151 L s P [(15)] .71904 .60885 0 2 Mshowa p .002 w .85296 .60885 m .85296 .6151 L s P [(17.5)] .85296 .60885 0 2 Mshowa p .002 w .98689 .60885 m .98689 .6151 L s P [(20)] .98689 .60885 0 2 Mshowa p .001 w .07619 .60885 m .07619 .6126 L s P p .001 w .10298 .60885 m .10298 .6126 L s P p .001 w .12976 .60885 m .12976 .6126 L s P p .001 w .15655 .60885 m .15655 .6126 L s P p .001 w .21012 .60885 m .21012 .6126 L s P p .001 w .2369 .60885 m .2369 .6126 L s P p .001 w .26369 .60885 m .26369 .6126 L s P p .001 w .29047 .60885 m .29047 .6126 L s P p .001 w .34404 .60885 m .34404 .6126 L s P p .001 w .37083 .60885 m .37083 .6126 L s P p .001 w .39761 .60885 m .39761 .6126 L s P p .001 w .4244 .60885 m .4244 .6126 L s P p .001 w .47797 .60885 m .47797 .6126 L s P p .001 w .50475 .60885 m .50475 .6126 L s P p .001 w .53154 .60885 m .53154 .6126 L s P p .001 w .55832 .60885 m .55832 .6126 L s P p .001 w .6119 .60885 m .6119 .6126 L s P p .001 w .63868 .60885 m .63868 .6126 L s P p .001 w .66547 .60885 m .66547 .6126 L s P p .001 w .69225 .60885 m .69225 .6126 L s P p .001 w .74582 .60885 m .74582 .6126 L s P p .001 w .77261 .60885 m .77261 .6126 L s P p .001 w .79939 .60885 m .79939 .6126 L s P p .001 w .82618 .60885 m .82618 .6126 L s P p .001 w .87975 .60885 m .87975 .6126 L s P p .001 w .90653 .60885 m .90653 .6126 L s P p .001 w .93332 .60885 m .93332 .6126 L s P p .001 w .9601 .60885 m .9601 .6126 L s P p .001 w .02262 .60885 m .02262 .6126 L s P p .002 w 0 .60885 m 1 .60885 L s P p .002 w .18333 .01433 m .18958 .01433 L s P [(-4)] .17083 .01433 1 0 Mshowa p .002 w .18333 .16296 m .18958 .16296 L s P [(-3)] .17083 .16296 1 0 Mshowa p .002 w .18333 .31159 m .18958 .31159 L s P [(-2)] .17083 .31159 1 0 Mshowa p .002 w .18333 .46022 m .18958 .46022 L s P [(-1)] .17083 .46022 1 0 Mshowa p .001 w .18333 .04406 m .18708 .04406 L s P p .001 w .18333 .07378 m .18708 .07378 L s P p .001 w .18333 .10351 m .18708 .10351 L s P p .001 w .18333 .13323 m .18708 .13323 L s P p .001 w .18333 .19269 m .18708 .19269 L s P p .001 w .18333 .22241 m .18708 .22241 L s P p .001 w .18333 .25214 m .18708 .25214 L s P p .001 w .18333 .28186 m .18708 .28186 L s P p .001 w .18333 .34132 m .18708 .34132 L s P p .001 w .18333 .37104 m .18708 .37104 L s P p .001 w .18333 .40077 m .18708 .40077 L s P p .001 w .18333 .43049 m .18708 .43049 L s P p .001 w .18333 .48994 m .18708 .48994 L s P p .001 w .18333 .51967 m .18708 .51967 L s P p .001 w .18333 .5494 m .18708 .5494 L s P p .001 w .18333 .57912 m .18708 .57912 L s P p .002 w .18333 0 m .18333 .61803 L s P P p 1 0 0 r p .02 w .36063 .5754 Mdot .53988 .59222 Mdot .83916 .60147 Mdot .78425 .60047 Mdot .64001 .59663 Mdot .15212 .48692 Mdot .88916 .60224 Mdot .67029 .59763 Mdot .62106 .59594 Mdot .24164 .54552 Mdot .93827 .60288 Mdot .43544 .58457 Mdot .60025 .59511 Mdot .88246 .60214 Mdot .74798 .59969 Mdot .12082 .44612 Mdot .13058 .46078 Mdot .95559 .60309 Mdot .22831 .53988 Mdot .30792 .56551 Mdot .59414 .59486 Mdot .0889 .37953 Mdot .64277 .59673 Mdot .64035 .59664 Mdot .06514 .2997 Mdot .74478 .59962 Mdot .30405 .56462 Mdot .64658 .59686 Mdot .19608 .52271 Mdot .67787 .59786 Mdot .95466 .60308 Mdot .39257 .57985 Mdot .38313 .57863 Mdot .32629 .5694 Mdot .72725 .5992 Mdot .67283 .5977 Mdot .68474 .59806 Mdot .91494 .60259 Mdot .40331 .58115 Mdot .29098 .5614 Mdot .22985 .54056 Mdot .59999 .5951 Mdot .26208 .55293 Mdot .23004 .54065 Mdot .86584 .6019 Mdot .89064 .60226 Mdot .16121 .49593 Mdot .24043 .54504 Mdot .46336 .58705 Mdot .31238 .5665 Mdot .08267 .36185 Mdot .04402 .18832 Mdot .83971 .60148 Mdot .43092 .58413 Mdot .60025 .59511 Mdot .22349 .53765 Mdot .7213 .59905 Mdot .31974 .56807 Mdot .58271 .59435 Mdot .10193 .41087 Mdot .14197 .47554 Mdot .77644 .60031 Mdot .08193 .35962 Mdot .75503 .59985 Mdot .22911 .54023 Mdot .34674 .57315 Mdot .61575 .59574 Mdot .14587 .48008 Mdot .88134 .60213 Mdot .81301 .60102 Mdot .32125 .56838 Mdot .7337 .59936 Mdot .20189 .52624 Mdot .94415 .60295 Mdot .14259 .47627 Mdot .21174 .53175 Mdot .90037 .6024 Mdot .85189 .60168 Mdot .30724 .56535 Mdot .72864 .59923 Mdot .41279 .58223 Mdot .90252 .60242 Mdot .11416 .43484 Mdot .19144 .51972 Mdot .89916 .60238 Mdot .23417 .54244 Mdot .77919 .60037 Mdot .04505 .19504 Mdot .86571 .6019 Mdot .18005 .51172 Mdot .58109 .59428 Mdot .43829 .58484 Mdot .19092 .51938 Mdot .52266 .59123 Mdot .09938 .40527 Mdot .35384 .57433 Mdot .18711 .5168 Mdot .14229 .47592 Mdot .39481 .58013 Mdot .11986 .44457 Mdot .51599 .59082 Mdot .9267 .60274 Mdot .92179 .60268 Mdot .51103 .59051 Mdot .75349 .59982 Mdot .97526 .60331 Mdot .67672 .59782 Mdot .64404 .59677 Mdot .51223 .59058 Mdot .52436 .59133 Mdot .14004 .47319 Mdot .42712 .58374 Mdot .61093 .59555 Mdot .4873 .5889 Mdot .67603 .5978 Mdot .96316 .60317 Mdot .90944 .60252 Mdot .97479 .6033 Mdot .25533 .55063 Mdot .35337 .57425 Mdot .72999 .59927 Mdot .8159 .60107 Mdot .61257 .59561 Mdot .4114 .58208 Mdot .02426 .01965 Mdot .91099 .60254 Mdot .33511 .57108 Mdot .93472 .60284 Mdot .41257 .58221 Mdot .94682 .60299 Mdot .17986 .51157 Mdot .61607 .59575 Mdot .43676 .58469 Mdot .78706 .60053 Mdot .81107 .60098 Mdot .39828 .58055 Mdot .80943 .60095 Mdot .47418 .58792 Mdot .49694 .58958 Mdot .52635 .59145 Mdot .10699 .42136 Mdot .56642 .59359 Mdot .91403 .60258 Mdot .93724 .60287 Mdot .49989 .58978 Mdot .17279 .50605 Mdot .67571 .59779 Mdot .7614 .59999 Mdot .85316 .6017 Mdot .30984 .56594 Mdot .37884 .57805 Mdot .77038 .60019 Mdot .29447 .56229 Mdot .22658 .53909 Mdot .36733 .57642 Mdot .27623 .55734 Mdot .85433 .60172 Mdot .12606 .45424 Mdot .20124 .52586 Mdot .35154 .57395 Mdot .33497 .57106 Mdot .06433 .29629 Mdot .94301 .60294 Mdot .28084 .55866 Mdot .67117 .59765 Mdot .07688 .34357 Mdot .32646 .56943 Mdot .78754 .60054 Mdot .45902 .58669 Mdot .6841 .59804 Mdot .28039 .55853 Mdot .35518 .57454 Mdot .93232 .60281 Mdot .31347 .56674 Mdot .97619 .60332 Mdot .05525 .25386 Mdot .32764 .56966 Mdot .96526 .6032 Mdot .02381 .01472 Mdot .20034 .52532 Mdot .55168 .59285 Mdot .50313 .58999 Mdot .29391 .56215 Mdot .17057 .50422 Mdot .34848 .57344 Mdot .84399 .60155 Mdot .35918 .57518 Mdot .81136 .60099 Mdot .05415 .24814 Mdot .83298 .60137 Mdot .05296 .2418 Mdot .65526 .59715 Mdot .48817 .58896 Mdot .17475 .50763 Mdot .21688 .53442 Mdot .0268 .04648 Mdot .85529 .60173 Mdot .1697 .50348 Mdot .51873 .59099 Mdot .34523 .57289 Mdot P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] No Input Form was saved for this expression. ;[o] -Graphics- :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] First test, radius fixed, fit for Hamaker and angle, the goal is A->1, a->0.9714: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={R->100}; sol=Sort[Join[NonlinearFit[dat1,ptmodel/.cond,d,{a,A}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp1,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-100)/100,(A R-100)/100}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.608613 0.145296 [ [(2.5)] .05026 .60861 0 2 Msboxa [(7.5)] .31481 .60861 0 2 Msboxa [(10)] .44709 .60861 0 2 Msboxa [(12.5)] .57937 .60861 0 2 Msboxa [(15)] .71164 .60861 0 2 Msboxa [(17.5)] .84392 .60861 0 2 Msboxa [(20)] .97619 .60861 0 2 Msboxa [(-4)] .17004 .02743 1 0 Msboxa [(-3)] .17004 .17272 1 0 Msboxa [(-2)] .17004 .31802 1 0 Msboxa [(-1)] .17004 .46332 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .60861 m .05026 .61486 L s P [(2.5)] .05026 .60861 0 2 Mshowa p .002 w .31481 .60861 m .31481 .61486 L s P [(7.5)] .31481 .60861 0 2 Mshowa p .002 w .44709 .60861 m .44709 .61486 L s P [(10)] .44709 .60861 0 2 Mshowa p .002 w .57937 .60861 m .57937 .61486 L s P [(12.5)] .57937 .60861 0 2 Mshowa p .002 w .71164 .60861 m .71164 .61486 L s P [(15)] .71164 .60861 0 2 Mshowa p .002 w .84392 .60861 m .84392 .61486 L s P [(17.5)] .84392 .60861 0 2 Mshowa p .002 w .97619 .60861 m .97619 .61486 L s P [(20)] .97619 .60861 0 2 Mshowa p .001 w .07672 .60861 m .07672 .61236 L s P p .001 w .10317 .60861 m .10317 .61236 L s P p .001 w .12963 .60861 m .12963 .61236 L s P p .001 w .15608 .60861 m .15608 .61236 L s P p .001 w .20899 .60861 m .20899 .61236 L s P p .001 w .23545 .60861 m .23545 .61236 L s P p .001 w .2619 .60861 m .2619 .61236 L s P p .001 w .28836 .60861 m .28836 .61236 L s P p .001 w .34127 .60861 m .34127 .61236 L s P p .001 w .36772 .60861 m .36772 .61236 L s P p .001 w .39418 .60861 m .39418 .61236 L s P p .001 w .42063 .60861 m .42063 .61236 L s P p .001 w .47354 .60861 m .47354 .61236 L s P p .001 w .5 .60861 m .5 .61236 L s P p .001 w .52646 .60861 m .52646 .61236 L s P p .001 w .55291 .60861 m .55291 .61236 L s P p .001 w .60582 .60861 m .60582 .61236 L s P p .001 w .63228 .60861 m .63228 .61236 L s P p .001 w .65873 .60861 m .65873 .61236 L s P p .001 w .68519 .60861 m .68519 .61236 L s P p .001 w .7381 .60861 m .7381 .61236 L s P p .001 w .76455 .60861 m .76455 .61236 L s P p .001 w .79101 .60861 m .79101 .61236 L s P p .001 w .81746 .60861 m .81746 .61236 L s P p .001 w .87037 .60861 m .87037 .61236 L s P p .001 w .89683 .60861 m .89683 .61236 L s P p .001 w .92328 .60861 m .92328 .61236 L s P p .001 w .94974 .60861 m .94974 .61236 L s P p .001 w .02381 .60861 m .02381 .61236 L s P p .002 w 0 .60861 m 1 .60861 L s P p .002 w .18254 .02743 m .18879 .02743 L s P [(-4)] .17004 .02743 1 0 Mshowa p .002 w .18254 .17272 m .18879 .17272 L s P [(-3)] .17004 .17272 1 0 Mshowa p .002 w .18254 .31802 m .18879 .31802 L s P [(-2)] .17004 .31802 1 0 Mshowa p .002 w .18254 .46332 m .18879 .46332 L s P [(-1)] .17004 .46332 1 0 Mshowa p .001 w .18254 .05649 m .18629 .05649 L s P p .001 w .18254 .08555 m .18629 .08555 L s P p .001 w .18254 .11461 m .18629 .11461 L s P p .001 w .18254 .14366 m .18629 .14366 L s P p .001 w .18254 .20178 m .18629 .20178 L s P p .001 w .18254 .23084 m .18629 .23084 L s P p .001 w .18254 .2599 m .18629 .2599 L s P p .001 w .18254 .28896 m .18629 .28896 L s P p .001 w .18254 .34708 m .18629 .34708 L s P p .001 w .18254 .37614 m .18629 .37614 L s P p .001 w .18254 .4052 m .18629 .4052 L s P p .001 w .18254 .43426 m .18629 .43426 L s P p .001 w .18254 .49238 m .18629 .49238 L s P p .001 w .18254 .52143 m .18629 .52143 L s P p .001 w .18254 .55049 m .18629 .55049 L s P p .001 w .18254 .57955 m .18629 .57955 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .11301 L .04365 .18889 L .05357 .24868 L .06349 .29661 L .07341 .33562 L .08333 .3678 L .10317 .41726 L .12302 .45301 L .14286 .47966 L .1627 .50006 L .18254 .51602 L .20238 .52873 L .22222 .53902 L .24206 .54746 L .2619 .55447 L .30159 .56533 L .34127 .57326 L .38095 .57922 L .42063 .5838 L .46032 .5874 L .5 .59028 L .53968 .59262 L .57937 .59455 L .61905 .59615 L .65873 .5975 L .69841 .59864 L .7381 .59961 L .77778 .60046 L .81746 .60119 L .85714 .60183 L .89683 .60239 L .93651 .60288 L .97619 .60332 L s P P p 1 0 0 r p .02 w .25223 .5512 Mdot .731 .59945 Mdot .52788 .59198 Mdot .03592 .13139 Mdot .61243 .5959 Mdot .11067 .43207 Mdot .90303 .60247 Mdot .50855 .59083 Mdot .15337 .49112 Mdot .62353 .59631 Mdot .93741 .60289 Mdot .88588 .60224 Mdot .70899 .59891 Mdot .17869 .51321 Mdot .20492 .53017 Mdot .16659 .5035 Mdot .56869 .59407 Mdot .07973 .35679 Mdot .64184 .59695 Mdot .54737 .59303 Mdot .51228 .59106 Mdot .90004 .60243 Mdot .55106 .59321 Mdot .03796 .14772 Mdot .81826 .6012 Mdot .83546 .60149 Mdot .63715 .59679 Mdot .69269 .59848 Mdot .09228 .3922 Mdot .63677 .59678 Mdot .42662 .5844 Mdot .11321 .4367 Mdot .89053 .6023 Mdot .18936 .52071 Mdot .36328 .57676 Mdot .17521 .51055 Mdot .52121 .59159 Mdot .31593 .56848 Mdot .86025 .60187 Mdot .59943 .59539 Mdot .79178 .60073 Mdot .93885 .60291 Mdot .38896 .58024 Mdot .84946 .60171 Mdot .47462 .58851 Mdot .68391 .59824 Mdot .13515 .47018 Mdot .53251 .59223 Mdot .59013 .59501 Mdot .24039 .5468 Mdot .70608 .59884 Mdot .61714 .59608 Mdot .89049 .6023 Mdot .65448 .59736 Mdot .52534 .59183 Mdot .78121 .60052 Mdot .07452 .33954 Mdot .12112 .45003 Mdot .91556 .60263 Mdot .5648 .59388 Mdot .34924 .57459 Mdot .34822 .57442 Mdot .61935 .59616 Mdot .70726 .59887 Mdot .70505 .59881 Mdot .07659 .3466 Mdot .20352 .52938 Mdot .27157 .55746 Mdot .1528 .49054 Mdot .17426 .50982 Mdot .21343 .53472 Mdot .58566 .59482 Mdot .84039 .60157 Mdot .26125 .55425 Mdot .24485 .54852 Mdot .53018 .59211 Mdot .51599 .59128 Mdot .4774 .58872 Mdot .33628 .57239 Mdot .22963 .54236 Mdot .20729 .53147 Mdot .07918 .35505 Mdot .97004 .60325 Mdot .46971 .58814 Mdot .15347 .49122 Mdot .23227 .54349 Mdot .88481 .60222 Mdot .67828 .59808 Mdot .39263 .58069 Mdot .22469 .54016 Mdot .08016 .35815 Mdot .97125 .60327 Mdot .84907 .6017 Mdot .91588 .60263 Mdot .85521 .6018 Mdot .46998 .58816 Mdot .72555 .59932 Mdot .26152 .55434 Mdot .81693 .60118 Mdot .55943 .59363 Mdot .44811 .58638 Mdot .35854 .57605 Mdot .39176 .58058 Mdot .36348 .57679 Mdot .51889 .59146 Mdot .42379 .58412 Mdot .63986 .59688 Mdot .93001 .6028 Mdot .8905 .6023 Mdot .6923 .59847 Mdot .09738 .4045 Mdot .15668 .49441 Mdot .62546 .59638 Mdot .77433 .60039 Mdot .2575 .55302 Mdot .43465 .58517 Mdot .71481 .59906 Mdot .03467 .12105 Mdot .09876 .40765 Mdot .08086 .36031 Mdot .18445 .51737 Mdot .36475 .57698 Mdot .11088 .43246 Mdot .23188 .54333 Mdot .46161 .58751 Mdot .43374 .58508 Mdot .22792 .54161 Mdot .97528 .60331 Mdot .4893 .58957 Mdot .96442 .60319 Mdot .41435 .58315 Mdot .90738 .60252 Mdot .25372 .55173 Mdot .56644 .59396 Mdot .16665 .50355 Mdot .35836 .57603 Mdot .68123 .59816 Mdot .90081 .60244 Mdot .47712 .5887 Mdot .90172 .60245 Mdot .96401 .60319 Mdot .91349 .6026 Mdot .57346 .59428 Mdot .02992 .07837 Mdot .34356 .57365 Mdot .43968 .58563 Mdot .20072 .52777 Mdot .10512 .42128 Mdot .49568 .59 Mdot .1805 .51455 Mdot .04237 .18014 Mdot .25498 .55216 Mdot .05456 .25391 Mdot .31768 .56884 Mdot .45728 .58716 Mdot .31157 .56756 Mdot .13754 .47322 Mdot .84808 .60169 Mdot .91301 .60259 Mdot .76677 .60024 Mdot .44114 .58577 Mdot .61007 .59581 Mdot .42257 .584 Mdot .37891 .57895 Mdot .39182 .58059 Mdot .08504 .37275 Mdot .91073 .60257 Mdot .74966 .59987 Mdot .797 .60082 Mdot .87777 .60213 Mdot .86019 .60187 Mdot .13481 .46974 Mdot .44286 .58592 Mdot .50093 .59035 Mdot .0441 .19189 Mdot .14583 .48307 Mdot .62846 .59649 Mdot .08461 .37152 Mdot .69392 .59852 Mdot .31114 .56746 Mdot .87311 .60206 Mdot .40956 .58263 Mdot .03673 .13799 Mdot .29855 .56462 Mdot .57006 .59413 Mdot .77381 .60038 Mdot .54978 .59315 Mdot .65178 .59728 Mdot .8975 .60239 Mdot .59099 .59505 Mdot .22739 .54137 Mdot .30366 .56581 Mdot .33048 .57133 Mdot .87029 .60202 Mdot .31756 .56881 Mdot .59555 .59524 Mdot .72368 .59928 Mdot .79793 .60084 Mdot .19772 .52599 Mdot .16996 .50635 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] -8 -8 Errors:{-0.0000769532, 2.89402 10 , 0, 2.89402 10 } :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.057038651197105, A -> 1.000000028940222, R -> 100} ;[o] {a -> 1.05704, A -> 1., R -> 100} :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] angle fixed, fit for Hamaker and radius, the goal is A->1, R->20: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={a->1.05712}; sol=Sort[Join[NonlinearFit[dat1,ptmodel/.cond,d,{{R,50,150},{A,0.5,1.5}}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp1,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-100)/100,(A R-100)/100}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 299; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.608648 0.145406 [ [(2.5)] .05026 .60865 0 2 Msboxa [(7.5)] .31481 .60865 0 2 Msboxa [(10)] .44709 .60865 0 2 Msboxa [(12.5)] .57937 .60865 0 2 Msboxa [(15)] .71164 .60865 0 2 Msboxa [(17.5)] .84392 .60865 0 2 Msboxa [(20)] .97619 .60865 0 2 Msboxa [(-4)] .17004 .02702 1 0 Msboxa [(-3)] .17004 .17243 1 0 Msboxa [(-2)] .17004 .31784 1 0 Msboxa [(-1)] .17004 .46324 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .60865 m .05026 .6149 L s P [(2.5)] .05026 .60865 0 2 Mshowa p .002 w .31481 .60865 m .31481 .6149 L s P [(7.5)] .31481 .60865 0 2 Mshowa p .002 w .44709 .60865 m .44709 .6149 L s P [(10)] .44709 .60865 0 2 Mshowa p .002 w .57937 .60865 m .57937 .6149 L s P [(12.5)] .57937 .60865 0 2 Mshowa p .002 w .71164 .60865 m .71164 .6149 L s P [(15)] .71164 .60865 0 2 Mshowa p .002 w .84392 .60865 m .84392 .6149 L s P [(17.5)] .84392 .60865 0 2 Mshowa p .002 w .97619 .60865 m .97619 .6149 L s P [(20)] .97619 .60865 0 2 Mshowa p .001 w .07672 .60865 m .07672 .6124 L s P p .001 w .10317 .60865 m .10317 .6124 L s P p .001 w .12963 .60865 m .12963 .6124 L s P p .001 w .15608 .60865 m .15608 .6124 L s P p .001 w .20899 .60865 m .20899 .6124 L s P p .001 w .23545 .60865 m .23545 .6124 L s P p .001 w .2619 .60865 m .2619 .6124 L s P p .001 w .28836 .60865 m .28836 .6124 L s P p .001 w .34127 .60865 m .34127 .6124 L s P p .001 w .36772 .60865 m .36772 .6124 L s P p .001 w .39418 .60865 m .39418 .6124 L s P p .001 w .42063 .60865 m .42063 .6124 L s P p .001 w .47354 .60865 m .47354 .6124 L s P p .001 w .5 .60865 m .5 .6124 L s P p .001 w .52646 .60865 m .52646 .6124 L s P p .001 w .55291 .60865 m .55291 .6124 L s P p .001 w .60582 .60865 m .60582 .6124 L s P p .001 w .63228 .60865 m .63228 .6124 L s P p .001 w .65873 .60865 m .65873 .6124 L s P p .001 w .68519 .60865 m .68519 .6124 L s P p .001 w .7381 .60865 m .7381 .6124 L s P p .001 w .76455 .60865 m .76455 .6124 L s P p .001 w .79101 .60865 m .79101 .6124 L s P p .001 w .81746 .60865 m .81746 .6124 L s P p .001 w .87037 .60865 m .87037 .6124 L s P p .001 w .89683 .60865 m .89683 .6124 L s P p .001 w .92328 .60865 m .92328 .6124 L s P p .001 w .94974 .60865 m .94974 .6124 L s P p .001 w .02381 .60865 m .02381 .6124 L s P p .002 w 0 .60865 m 1 .60865 L s P p .002 w .18254 .02702 m .18879 .02702 L s P [(-4)] .17004 .02702 1 0 Mshowa p .002 w .18254 .17243 m .18879 .17243 L s P [(-3)] .17004 .17243 1 0 Mshowa p .002 w .18254 .31784 m .18879 .31784 L s P [(-2)] .17004 .31784 1 0 Mshowa p .002 w .18254 .46324 m .18879 .46324 L s P [(-1)] .17004 .46324 1 0 Mshowa p .001 w .18254 .0561 m .18629 .0561 L s P p .001 w .18254 .08519 m .18629 .08519 L s P p .001 w .18254 .11427 m .18629 .11427 L s P p .001 w .18254 .14335 m .18629 .14335 L s P p .001 w .18254 .20151 m .18629 .20151 L s P p .001 w .18254 .23059 m .18629 .23059 L s P p .001 w .18254 .25967 m .18629 .25967 L s P p .001 w .18254 .28875 m .18629 .28875 L s P p .001 w .18254 .34692 m .18629 .34692 L s P p .001 w .18254 .376 m .18629 .376 L s P p .001 w .18254 .40508 m .18629 .40508 L s P p .001 w .18254 .43416 m .18629 .43416 L s P p .001 w .18254 .49232 m .18629 .49232 L s P p .001 w .18254 .5214 m .18629 .5214 L s P p .001 w .18254 .55049 m .18629 .55049 L s P p .001 w .18254 .57957 m .18629 .57957 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .11294 L .04365 .18878 L .05357 .24853 L .06349 .29644 L .07341 .33544 L .08333 .3676 L .10317 .41706 L .12302 .45281 L .14286 .47947 L .1627 .49988 L .18254 .51585 L .20238 .52857 L .22222 .53886 L .24206 .54731 L .2619 .55433 L .30159 .56521 L .34127 .57316 L .38095 .57913 L .42063 .58372 L .46032 .58734 L .5 .59023 L .53968 .59257 L .57937 .5945 L .61905 .59611 L .65873 .59746 L .69841 .59861 L .7381 .59959 L .77778 .60044 L .81746 .60117 L .85714 .60182 L .89683 .60238 L .93651 .60288 L .97619 .60332 L s P P p 1 0 0 r p .02 w .35766 .57593 Mdot .53469 .59238 Mdot .83029 .60143 Mdot .77605 .60045 Mdot .63359 .5967 Mdot .15171 .48936 Mdot .87967 .60218 Mdot .6635 .59767 Mdot .61487 .59602 Mdot .24013 .54669 Mdot .92817 .60281 Mdot .43154 .58489 Mdot .59432 .59521 Mdot .87304 .60209 Mdot .74022 .59969 Mdot .1208 .44945 Mdot .13044 .46379 Mdot .94528 .60301 Mdot .22697 .54117 Mdot .30559 .56625 Mdot .58829 .59496 Mdot .08927 .3843 Mdot .63631 .59679 Mdot .63392 .59671 Mdot .0658 .3062 Mdot .73707 .59962 Mdot .30177 .56538 Mdot .64008 .59692 Mdot .19513 .52438 Mdot .67098 .59789 Mdot .94436 .603 Mdot .3892 .58028 Mdot .37988 .57909 Mdot .32374 .57005 Mdot .71975 .59921 Mdot .666 .59775 Mdot .67776 .59809 Mdot .90513 .60253 Mdot .39981 .58155 Mdot .28886 .56223 Mdot .22848 .54184 Mdot .59407 .5952 Mdot .26032 .55394 Mdot .22867 .54193 Mdot .85663 .60185 Mdot .88113 .6022 Mdot .16069 .49818 Mdot .23894 .54622 Mdot .45911 .58733 Mdot .31 .56722 Mdot .08312 .36701 Mdot .04495 .19724 Mdot .83082 .60144 Mdot .42707 .58446 Mdot .59432 .59521 Mdot .22221 .53899 Mdot .71387 .59906 Mdot .31727 .56876 Mdot .577 .59447 Mdot .10214 .41497 Mdot .14169 .47823 Mdot .76833 .6003 Mdot .08239 .36483 Mdot .74719 .59985 Mdot .22775 .54152 Mdot .34393 .57372 Mdot .60963 .59582 Mdot .14554 .48267 Mdot .87194 .60207 Mdot .80446 .60099 Mdot .31876 .56906 Mdot .72613 .59936 Mdot .20087 .52783 Mdot .93398 .60288 Mdot .1423 .47894 Mdot .2106 .53323 Mdot .89074 .60234 Mdot .84286 .60164 Mdot .30492 .5661 Mdot .72113 .59924 Mdot .40917 .58261 Mdot .89286 .60236 Mdot .11422 .43841 Mdot .19054 .52146 Mdot .88954 .60232 Mdot .23275 .54368 Mdot .77105 .60035 Mdot .04596 .20381 Mdot .85651 .60185 Mdot .1793 .51362 Mdot .5754 .5944 Mdot .43436 .58516 Mdot .19004 .52112 Mdot .51769 .59141 Mdot .09963 .40948 Mdot .35095 .57488 Mdot .18628 .5186 Mdot .14201 .4786 Mdot .39141 .58055 Mdot .11986 .44793 Mdot .5111 .59101 Mdot .91675 .60267 Mdot .91189 .60261 Mdot .5062 .5907 Mdot .74567 .59981 Mdot .96471 .60323 Mdot .66985 .59786 Mdot .63756 .59683 Mdot .50738 .59078 Mdot .51937 .59151 Mdot .13978 .47593 Mdot .42332 .58409 Mdot .60487 .59564 Mdot .48276 .58913 Mdot .66917 .59784 Mdot .95275 .6031 Mdot .89969 .60245 Mdot .96424 .60322 Mdot .25365 .55169 Mdot .35048 .5748 Mdot .72246 .59928 Mdot .80731 .60104 Mdot .60649 .5957 Mdot .4078 .58246 Mdot .02543 .03223 Mdot .90123 .60247 Mdot .33245 .5717 Mdot .92466 .60277 Mdot .40895 .58258 Mdot .93662 .60291 Mdot .17911 .51348 Mdot .60994 .59583 Mdot .43284 .58502 Mdot .77883 .60051 Mdot .80254 .60095 Mdot .39484 .58097 Mdot .80092 .60092 Mdot .4698 .58817 Mdot .49228 .58979 Mdot .52133 .59162 Mdot .10714 .42522 Mdot .5609 .59372 Mdot .90423 .60251 Mdot .92715 .6028 Mdot .4952 .58999 Mdot .17213 .50808 Mdot .66885 .59783 Mdot .75348 .59998 Mdot .84411 .60166 Mdot .30749 .56667 Mdot .37564 .57852 Mdot .76235 .60017 Mdot .29231 .5631 Mdot .22525 .5404 Mdot .36427 .57692 Mdot .27429 .55825 Mdot .84526 .60167 Mdot .12597 .45739 Mdot .20023 .52746 Mdot .34868 .57451 Mdot .33231 .57168 Mdot .065 .30287 Mdot .93286 .60287 Mdot .27884 .55955 Mdot .66436 .5977 Mdot .0774 .34913 Mdot .3239 .57008 Mdot .7793 .60052 Mdot .45483 .58697 Mdot .67714 .59808 Mdot .27841 .55942 Mdot .35227 .57509 Mdot .92229 .60274 Mdot .31108 .56746 Mdot .96562 .60324 Mdot .05604 .26136 Mdot .32507 .57031 Mdot .95483 .60312 Mdot .02498 .0274 Mdot .19934 .52693 Mdot .54635 .593 Mdot .4984 .5902 Mdot .29175 .56296 Mdot .16994 .50629 Mdot .34565 .57401 Mdot .83505 .60151 Mdot .35622 .57571 Mdot .80282 .60096 Mdot .05495 .25576 Mdot .82418 .60133 Mdot .05378 .24955 Mdot .64865 .5972 Mdot .48362 .58919 Mdot .17406 .50962 Mdot .21568 .53583 Mdot .02794 .05848 Mdot .84622 .60169 Mdot .16907 .50557 Mdot .51381 .59118 Mdot .34244 .57347 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{0., -0.0998713, 0.10818, -0.00249542} :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.05712, A -> 0.900128720251945, R -> 110.8179925918564} ;[o] {a -> 1.05712, A -> 0.900129, R -> 110.818} :[font = subsubsection; inactive; ] the results are not very good, but the product R A is good: :[font = input; startGroup; nowordwrap; ] R A/.% :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] 99.7504578525973 ;[o] 99.7505 :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] fixing hamaker and fitting for the radius and angle is almost equivalent to fixing the radius, so it is not worth doing again. :[font = subsubsection; inactive; ] now everything free: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={}; sol=Sort[Join[NonlinearFit[dat1,ptmodel/.cond,d,{{a,0.9,1.1},{R,50,150},{A,0.5,1.5}}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp1,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-100)/100,(A R-100)/100}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 299; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.608602 0.145323 [ [(2.5)] .05026 .6086 0 2 Msboxa [(7.5)] .31481 .6086 0 2 Msboxa [(10)] .44709 .6086 0 2 Msboxa [(12.5)] .57937 .6086 0 2 Msboxa [(15)] .71164 .6086 0 2 Msboxa [(17.5)] .84392 .6086 0 2 Msboxa [(20)] .97619 .6086 0 2 Msboxa [(-4)] .17004 .02731 1 0 Msboxa [(-3)] .17004 .17263 1 0 Msboxa [(-2)] .17004 .31796 1 0 Msboxa [(-1)] .17004 .46328 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .6086 m .05026 .61485 L s P [(2.5)] .05026 .6086 0 2 Mshowa p .002 w .31481 .6086 m .31481 .61485 L s P [(7.5)] .31481 .6086 0 2 Mshowa p .002 w .44709 .6086 m .44709 .61485 L s P [(10)] .44709 .6086 0 2 Mshowa p .002 w .57937 .6086 m .57937 .61485 L s P [(12.5)] .57937 .6086 0 2 Mshowa p .002 w .71164 .6086 m .71164 .61485 L s P [(15)] .71164 .6086 0 2 Mshowa p .002 w .84392 .6086 m .84392 .61485 L s P [(17.5)] .84392 .6086 0 2 Mshowa p .002 w .97619 .6086 m .97619 .61485 L s P [(20)] .97619 .6086 0 2 Mshowa p .001 w .07672 .6086 m .07672 .61235 L s P p .001 w .10317 .6086 m .10317 .61235 L s P p .001 w .12963 .6086 m .12963 .61235 L s P p .001 w .15608 .6086 m .15608 .61235 L s P p .001 w .20899 .6086 m .20899 .61235 L s P p .001 w .23545 .6086 m .23545 .61235 L s P p .001 w .2619 .6086 m .2619 .61235 L s P p .001 w .28836 .6086 m .28836 .61235 L s P p .001 w .34127 .6086 m .34127 .61235 L s P p .001 w .36772 .6086 m .36772 .61235 L s P p .001 w .39418 .6086 m .39418 .61235 L s P p .001 w .42063 .6086 m .42063 .61235 L s P p .001 w .47354 .6086 m .47354 .61235 L s P p .001 w .5 .6086 m .5 .61235 L s P p .001 w .52646 .6086 m .52646 .61235 L s P p .001 w .55291 .6086 m .55291 .61235 L s P p .001 w .60582 .6086 m .60582 .61235 L s P p .001 w .63228 .6086 m .63228 .61235 L s P p .001 w .65873 .6086 m .65873 .61235 L s P p .001 w .68519 .6086 m .68519 .61235 L s P p .001 w .7381 .6086 m .7381 .61235 L s P p .001 w .76455 .6086 m .76455 .61235 L s P p .001 w .79101 .6086 m .79101 .61235 L s P p .001 w .81746 .6086 m .81746 .61235 L s P p .001 w .87037 .6086 m .87037 .61235 L s P p .001 w .89683 .6086 m .89683 .61235 L s P p .001 w .92328 .6086 m .92328 .61235 L s P p .001 w .94974 .6086 m .94974 .61235 L s P p .001 w .02381 .6086 m .02381 .61235 L s P p .002 w 0 .6086 m 1 .6086 L s P p .002 w .18254 .02731 m .18879 .02731 L s P [(-4)] .17004 .02731 1 0 Mshowa p .002 w .18254 .17263 m .18879 .17263 L s P [(-3)] .17004 .17263 1 0 Mshowa p .002 w .18254 .31796 m .18879 .31796 L s P [(-2)] .17004 .31796 1 0 Mshowa p .002 w .18254 .46328 m .18879 .46328 L s P [(-1)] .17004 .46328 1 0 Mshowa p .001 w .18254 .05638 m .18629 .05638 L s P p .001 w .18254 .08544 m .18629 .08544 L s P p .001 w .18254 .1145 m .18629 .1145 L s P p .001 w .18254 .14357 m .18629 .14357 L s P p .001 w .18254 .2017 m .18629 .2017 L s P p .001 w .18254 .23076 m .18629 .23076 L s P p .001 w .18254 .25983 m .18629 .25983 L s P p .001 w .18254 .28889 m .18629 .28889 L s P p .001 w .18254 .34702 m .18629 .34702 L s P p .001 w .18254 .37609 m .18629 .37609 L s P p .001 w .18254 .40515 m .18629 .40515 L s P p .001 w .18254 .43421 m .18629 .43421 L s P p .001 w .18254 .49234 m .18629 .49234 L s P p .001 w .18254 .52141 m .18629 .52141 L s P p .001 w .18254 .55047 m .18629 .55047 L s P p .001 w .18254 .57954 m .18629 .57954 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .11291 L .04365 .18872 L .05357 .24845 L .06349 .29636 L .07341 .33535 L .08333 .36751 L .10317 .41697 L .12302 .45272 L .14286 .47938 L .1627 .4998 L .18254 .51577 L .20238 .52849 L .22222 .5388 L .24206 .54725 L .2619 .55427 L .30159 .56516 L .34127 .57311 L .38095 .57909 L .42063 .58369 L .46032 .58731 L .5 .5902 L .53968 .59256 L .57937 .59449 L .61905 .5961 L .65873 .59745 L .69841 .5986 L .7381 .59959 L .77778 .60044 L .81746 .60117 L .85714 .60181 L .89683 .60238 L .93651 .60288 L .97619 .60332 L s P P p 1 0 0 r p .02 w .35766 .5759 Mdot .53469 .59234 Mdot .83029 .60139 Mdot .77605 .60041 Mdot .63359 .59666 Mdot .15171 .48938 Mdot .87967 .60214 Mdot .6635 .59763 Mdot .61487 .59598 Mdot .24013 .54668 Mdot .92817 .60277 Mdot .43154 .58486 Mdot .59432 .59517 Mdot .87304 .60205 Mdot .74022 .59965 Mdot .1208 .44949 Mdot .13044 .46383 Mdot .94528 .60297 Mdot .22697 .54116 Mdot .30559 .56623 Mdot .58829 .59492 Mdot .08927 .38438 Mdot .63631 .59675 Mdot .63392 .59667 Mdot .0658 .30633 Mdot .73707 .59958 Mdot .30177 .56536 Mdot .64008 .59688 Mdot .19513 .52438 Mdot .67098 .59785 Mdot .94436 .60296 Mdot .3892 .58025 Mdot .37988 .57906 Mdot .32374 .57003 Mdot .71975 .59917 Mdot .666 .59771 Mdot .67776 .59805 Mdot .90513 .60248 Mdot .39981 .58152 Mdot .28886 .56221 Mdot .22848 .54183 Mdot .59407 .59516 Mdot .26032 .55393 Mdot .22867 .54192 Mdot .85663 .60181 Mdot .88113 .60216 Mdot .16069 .49819 Mdot .23894 .54621 Mdot .45911 .58729 Mdot .31 .5672 Mdot .08312 .3671 Mdot .04495 .19743 Mdot .83082 .6014 Mdot .42707 .58443 Mdot .59432 .59517 Mdot .22221 .53899 Mdot .71387 .59902 Mdot .31727 .56873 Mdot .577 .59443 Mdot .10214 .41503 Mdot .14169 .47826 Mdot .76833 .60025 Mdot .08239 .36492 Mdot .74719 .59981 Mdot .22775 .54151 Mdot .34393 .5737 Mdot .60963 .59578 Mdot .14554 .4827 Mdot .87194 .60203 Mdot .80446 .60095 Mdot .31876 .56904 Mdot .72613 .59932 Mdot .20087 .52783 Mdot .93398 .60284 Mdot .1423 .47897 Mdot .2106 .53322 Mdot .89074 .60229 Mdot .84286 .60159 Mdot .30492 .56608 Mdot .72113 .5992 Mdot .40917 .58258 Mdot .89286 .60232 Mdot .11422 .43846 Mdot .19054 .52146 Mdot .88954 .60228 Mdot .23275 .54367 Mdot .77105 .60031 Mdot .04596 .204 Mdot .85651 .6018 Mdot .1793 .51363 Mdot .5754 .59436 Mdot .43436 .58513 Mdot .19004 .52113 Mdot .51769 .59137 Mdot .09963 .40955 Mdot .35095 .57485 Mdot .18628 .51861 Mdot .14201 .47863 Mdot .39141 .58052 Mdot .11986 .44797 Mdot .5111 .59097 Mdot .91675 .60263 Mdot .91189 .60257 Mdot .5062 .59067 Mdot .74567 .59977 Mdot .96471 .60319 Mdot .66985 .59782 Mdot .63756 .59679 Mdot .50738 .59074 Mdot .51937 .59147 Mdot .13978 .47596 Mdot .42332 .58406 Mdot .60487 .5956 Mdot .48276 .58909 Mdot .66917 .5978 Mdot .95275 .60305 Mdot .89969 .60241 Mdot .96424 .60318 Mdot .25365 .55168 Mdot .35048 .57477 Mdot .72246 .59923 Mdot .80731 .601 Mdot .60649 .59566 Mdot .4078 .58242 Mdot .02543 .03252 Mdot .90123 .60243 Mdot .33245 .57168 Mdot .92466 .60273 Mdot .40895 .58255 Mdot .93662 .60287 Mdot .17911 .51349 Mdot .60994 .59579 Mdot .43284 .58498 Mdot .77883 .60046 Mdot .80254 .60091 Mdot .39484 .58094 Mdot .80092 .60088 Mdot .4698 .58814 Mdot .49228 .58976 Mdot .52133 .59159 Mdot .10714 .42528 Mdot .5609 .59368 Mdot .90423 .60247 Mdot .92715 .60276 Mdot .4952 .58995 Mdot .17213 .50809 Mdot .66885 .59779 Mdot .75348 .59994 Mdot .84411 .60161 Mdot .30749 .56665 Mdot .37564 .57849 Mdot .76235 .60013 Mdot .29231 .56308 Mdot .22525 .54039 Mdot .36427 .57689 Mdot .27429 .55823 Mdot .84526 .60163 Mdot .12597 .45743 Mdot .20023 .52746 Mdot .34868 .57448 Mdot .33231 .57165 Mdot .065 .30299 Mdot .93286 .60283 Mdot .27884 .55953 Mdot .66436 .59766 Mdot .0774 .34923 Mdot .3239 .57006 Mdot .7793 .60047 Mdot .45483 .58694 Mdot .67714 .59803 Mdot .27841 .55941 Mdot .35227 .57506 Mdot .92229 .6027 Mdot .31108 .56743 Mdot .96562 .6032 Mdot .05604 .26151 Mdot .32507 .57029 Mdot .95483 .60308 Mdot .02498 .02769 Mdot .19934 .52693 Mdot .54635 .59296 Mdot .4984 .59017 Mdot .29175 .56294 Mdot .16994 .5063 Mdot .34565 .57398 Mdot .83505 .60147 Mdot .35622 .57568 Mdot .80282 .60092 Mdot .05495 .25591 Mdot .82418 .60129 Mdot .05378 .24971 Mdot .64865 .59716 Mdot .48362 .58916 Mdot .17406 .50963 Mdot .21568 .53583 Mdot .02794 .05875 Mdot .84622 .60165 Mdot .16907 .50558 Mdot .51381 .59114 Mdot .34244 .57344 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{0.0967424, -0.154929, 0.179961, -0.00284875} :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.159388319433704, A -> 0.845071046765538, R -> 117.9961443506819} ;[o] {a -> 1.15939, A -> 0.845071, R -> 117.996} :[font = subsubsection; inactive; ] and the product A R is very good: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] R A/.% :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] 99.7151252207282 ;[o] 99.7151 :[font = subsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 12; ] dat2, data with 10% noise : :[font = input; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] dat2=Table[Module[{temp}, temp=Random[Real,{2,20}]; {temp,noise*0.4+ForceParametricTip[1.05712,temp,1,100]}], {x,1,200}]; :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] lp2=ListPlot[dat2,PlotRange->All,PlotStyle->{Red,PointSize[0.02]}] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 299; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0908631 0.0535817 0.554558 0.143563 [ [(2.5)] .04309 .55456 0 2 Msboxa [(7.5)] .311 .55456 0 2 Msboxa [(10)] .44495 .55456 0 2 Msboxa [(12.5)] .57891 .55456 0 2 Msboxa [(15)] .71286 .55456 0 2 Msboxa [(17.5)] .84682 .55456 0 2 Msboxa [(20)] .98077 .55456 0 2 Msboxa [(-3)] .16455 .12387 1 0 Msboxa [(-2)] .16455 .26743 1 0 Msboxa [(-1)] .16455 .411 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .04309 .55456 m .04309 .56081 L s P [(2.5)] .04309 .55456 0 2 Mshowa p .002 w .311 .55456 m .311 .56081 L s P [(7.5)] .311 .55456 0 2 Mshowa p .002 w .44495 .55456 m .44495 .56081 L s P [(10)] .44495 .55456 0 2 Mshowa p .002 w .57891 .55456 m .57891 .56081 L s P [(12.5)] .57891 .55456 0 2 Mshowa p .002 w .71286 .55456 m .71286 .56081 L s P [(15)] .71286 .55456 0 2 Mshowa p .002 w .84682 .55456 m .84682 .56081 L s P [(17.5)] .84682 .55456 0 2 Mshowa p .002 w .98077 .55456 m .98077 .56081 L s P [(20)] .98077 .55456 0 2 Mshowa p .001 w .06988 .55456 m .06988 .55831 L s P p .001 w .09667 .55456 m .09667 .55831 L s P p .001 w .12346 .55456 m .12346 .55831 L s P p .001 w .15025 .55456 m .15025 .55831 L s P p .001 w .20384 .55456 m .20384 .55831 L s P p .001 w .23063 .55456 m .23063 .55831 L s P p .001 w .25742 .55456 m .25742 .55831 L s P p .001 w .28421 .55456 m .28421 .55831 L s P p .001 w .33779 .55456 m .33779 .55831 L s P p .001 w .36458 .55456 m .36458 .55831 L s P p .001 w .39137 .55456 m .39137 .55831 L s P p .001 w .41816 .55456 m .41816 .55831 L s P p .001 w .47175 .55456 m .47175 .55831 L s P p .001 w .49854 .55456 m .49854 .55831 L s P p .001 w .52533 .55456 m .52533 .55831 L s P p .001 w .55212 .55456 m .55212 .55831 L s P p .001 w .6057 .55456 m .6057 .55831 L s P p .001 w .63249 .55456 m .63249 .55831 L s P p .001 w .65928 .55456 m .65928 .55831 L s P p .001 w .68607 .55456 m .68607 .55831 L s P p .001 w .73965 .55456 m .73965 .55831 L s P p .001 w .76644 .55456 m .76644 .55831 L s P p .001 w .79324 .55456 m .79324 .55831 L s P p .001 w .82003 .55456 m .82003 .55831 L s P p .001 w .87361 .55456 m .87361 .55831 L s P p .001 w .9004 .55456 m .9004 .55831 L s P p .001 w .92719 .55456 m .92719 .55831 L s P p .001 w .95398 .55456 m .95398 .55831 L s P p .001 w .0163 .55456 m .0163 .55831 L s P p .002 w 0 .55456 m 1 .55456 L s P p .002 w .17705 .12387 m .1833 .12387 L s P [(-3)] .16455 .12387 1 0 Mshowa p .002 w .17705 .26743 m .1833 .26743 L s P [(-2)] .16455 .26743 1 0 Mshowa p .002 w .17705 .411 m .1833 .411 L s P [(-1)] .16455 .411 1 0 Mshowa p .001 w .17705 .15258 m .1808 .15258 L s P p .001 w .17705 .18129 m .1808 .18129 L s P p .001 w .17705 .21001 m .1808 .21001 L s P p .001 w .17705 .23872 m .1808 .23872 L s P p .001 w .17705 .29615 m .1808 .29615 L s P p .001 w .17705 .32486 m .1808 .32486 L s P p .001 w .17705 .35357 m .1808 .35357 L s P p .001 w .17705 .38228 m .1808 .38228 L s P p .001 w .17705 .43971 m .1808 .43971 L s P p .001 w .17705 .46842 m .1808 .46842 L s P p .001 w .17705 .49713 m .1808 .49713 L s P p .001 w .17705 .52585 m .1808 .52585 L s P p .001 w .17705 .09516 m .1808 .09516 L s P p .001 w .17705 .06644 m .1808 .06644 L s P p .001 w .17705 .03773 m .1808 .03773 L s P p .001 w .17705 .00902 m .1808 .00902 L s P p .001 w .17705 .58327 m .1808 .58327 L s P p .001 w .17705 .61198 m .1808 .61198 L s P p .002 w .17705 0 m .17705 .61803 L s P P p 1 0 0 r p .02 w .23053 .52891 Mdot .66284 .59983 Mdot .61377 .53701 Mdot .32753 .48267 Mdot .41478 .56448 Mdot .04721 .20203 Mdot .03974 .13678 Mdot .62627 .57378 Mdot .87864 .56036 Mdot .58271 .57309 Mdot .25091 .45459 Mdot .36613 .57459 Mdot .03668 .11739 Mdot .68406 .59621 Mdot .40369 .56665 Mdot .37283 .55451 Mdot .96968 .49512 Mdot .34191 .49346 Mdot .94624 .57654 Mdot .69641 .54392 Mdot .0839 .29431 Mdot .13 .44343 Mdot .81376 .49196 Mdot .74465 .57701 Mdot .79338 .57337 Mdot .07689 .35425 Mdot .40599 .57473 Mdot .68483 .49396 Mdot .41709 .57254 Mdot .35922 .55688 Mdot .45161 .49067 Mdot .64358 .52028 Mdot .38401 .58339 Mdot .52988 .54846 Mdot .55102 .59424 Mdot .766 .58259 Mdot .73841 .51645 Mdot .70541 .6019 Mdot .34872 .50369 Mdot .03688 .19299 Mdot .9124 .54568 Mdot .65843 .50217 Mdot .47709 .51709 Mdot .03115 .1419 Mdot .10938 .37145 Mdot .48205 .50661 Mdot .53914 .52374 Mdot .69682 .53696 Mdot .7815 .50325 Mdot .97218 .54194 Mdot .44908 .56469 Mdot .9516 .54917 Mdot .51745 .51489 Mdot .30947 .49846 Mdot .05666 .29963 Mdot .29462 .51379 Mdot .92805 .56157 Mdot .79335 .52074 Mdot .40521 .49914 Mdot .11283 .43379 Mdot .60449 .498 Mdot .12142 .39654 Mdot .17171 .44112 Mdot .15059 .4877 Mdot .63503 .48996 Mdot .82189 .55195 Mdot .59468 .49327 Mdot .54357 .48361 Mdot .6474 .53977 Mdot .73099 .57351 Mdot .25849 .47009 Mdot .63447 .58915 Mdot .63478 .49776 Mdot .52935 .53824 Mdot .47937 .56568 Mdot .39506 .53875 Mdot .82511 .57322 Mdot .55395 .48857 Mdot .24674 .51407 Mdot .02668 .01472 Mdot .5801 .50326 Mdot .27646 .483 Mdot .33791 .57092 Mdot .62277 .52984 Mdot .6839 .58316 Mdot .10972 .43454 Mdot .22083 .43835 Mdot .69543 .5209 Mdot .37649 .50777 Mdot .15778 .41786 Mdot .14606 .46036 Mdot .1474 .40374 Mdot .54673 .54254 Mdot .85171 .59412 Mdot .22511 .49878 Mdot .24524 .49851 Mdot .52198 .56592 Mdot .15182 .45666 Mdot .31745 .53565 Mdot .43717 .51212 Mdot .92173 .52849 Mdot .29569 .46439 Mdot .79197 .56063 Mdot .16459 .49882 Mdot .26155 .45544 Mdot .29365 .56699 Mdot .05273 .23075 Mdot .06471 .27557 Mdot .51152 .56762 Mdot .89366 .5889 Mdot .21037 .4358 Mdot .47279 .53756 Mdot .26941 .5355 Mdot .19341 .4659 Mdot .45821 .49312 Mdot .04512 .19664 Mdot .21297 .43279 Mdot .73223 .55158 Mdot .17654 .46811 Mdot .68382 .49576 Mdot .64579 .5756 Mdot .77093 .51465 Mdot .45172 .55484 Mdot .31444 .53825 Mdot .19861 .52478 Mdot .13733 .40074 Mdot .72116 .57663 Mdot .10851 .41332 Mdot .52449 .56297 Mdot .35705 .48487 Mdot .36425 .48584 Mdot .654 .49673 Mdot .69924 .53183 Mdot .86385 .59032 Mdot .26382 .52581 Mdot .56571 .50218 Mdot .08151 .36209 Mdot .44468 .57896 Mdot .34112 .46324 Mdot .35247 .48529 Mdot .7974 .52234 Mdot .97619 .49296 Mdot .44944 .48881 Mdot .33849 .56753 Mdot .73098 .57168 Mdot .45541 .48132 Mdot .48346 .48145 Mdot .87047 .59127 Mdot .41825 .5016 Mdot .44209 .58648 Mdot .09343 .32974 Mdot .10592 .41657 Mdot .27681 .55966 Mdot .1105 .42329 Mdot .80814 .58664 Mdot .75279 .58795 Mdot .09346 .31473 Mdot .31367 .55183 Mdot .59077 .49372 Mdot .42397 .58139 Mdot .18882 .50921 Mdot .96265 .60312 Mdot .11169 .40283 Mdot .87303 .51113 Mdot .81565 .56156 Mdot .77883 .53463 Mdot .02381 .07495 Mdot .04235 .18004 Mdot .91112 .56734 Mdot .70944 .56934 Mdot .33665 .52824 Mdot .30177 .54375 Mdot .16413 .48116 Mdot .31989 .55089 Mdot .06875 .25275 Mdot .42763 .54491 Mdot .23387 .49622 Mdot .62957 .51098 Mdot .22636 .51841 Mdot .60352 .56388 Mdot .29601 .46384 Mdot .87485 .60332 Mdot .94014 .54991 Mdot .58938 .5079 Mdot .7923 .57739 Mdot .28579 .49875 Mdot .73986 .5585 Mdot .83894 .58076 Mdot .52229 .48986 Mdot .22567 .49692 Mdot .31222 .49571 Mdot .60292 .58428 Mdot .03251 .08401 Mdot .70884 .58975 Mdot .07315 .33378 Mdot .13576 .44647 Mdot .26162 .55649 Mdot .83074 .5208 Mdot .50253 .52058 Mdot .97258 .54684 Mdot P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] No Input Form was saved for this expression. ;[o] -Graphics- :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] First test, radius fixed, fit for Hamaker and angle, the goal is A->1, a->1.05712: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={R->100}; sol=Sort[Join[NonlinearFit[dat2,ptmodel/.cond,d,{{a,0.9,1.1},{A,0.5,1.5}}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp2,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-100)/100,(A R-100)/100}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 299; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.558003 0.133422 [ [(2.5)] .05026 .558 0 2 Msboxa [(7.5)] .31481 .558 0 2 Msboxa [(10)] .44709 .558 0 2 Msboxa [(12.5)] .57937 .558 0 2 Msboxa [(15)] .71164 .558 0 2 Msboxa [(17.5)] .84392 .558 0 2 Msboxa [(20)] .97619 .558 0 2 Msboxa [(-4)] .17004 .02432 1 0 Msboxa [(-3)] .17004 .15774 1 0 Msboxa [(-2)] .17004 .29116 1 0 Msboxa [(-1)] .17004 .42458 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .558 m .05026 .56425 L s P [(2.5)] .05026 .558 0 2 Mshowa p .002 w .31481 .558 m .31481 .56425 L s P [(7.5)] .31481 .558 0 2 Mshowa p .002 w .44709 .558 m .44709 .56425 L s P [(10)] .44709 .558 0 2 Mshowa p .002 w .57937 .558 m .57937 .56425 L s P [(12.5)] .57937 .558 0 2 Mshowa p .002 w .71164 .558 m .71164 .56425 L s P [(15)] .71164 .558 0 2 Mshowa p .002 w .84392 .558 m .84392 .56425 L s P [(17.5)] .84392 .558 0 2 Mshowa p .002 w .97619 .558 m .97619 .56425 L s P [(20)] .97619 .558 0 2 Mshowa p .001 w .07672 .558 m .07672 .56175 L s P p .001 w .10317 .558 m .10317 .56175 L s P p .001 w .12963 .558 m .12963 .56175 L s P p .001 w .15608 .558 m .15608 .56175 L s P p .001 w .20899 .558 m .20899 .56175 L s P p .001 w .23545 .558 m .23545 .56175 L s P p .001 w .2619 .558 m .2619 .56175 L s P p .001 w .28836 .558 m .28836 .56175 L s P p .001 w .34127 .558 m .34127 .56175 L s P p .001 w .36772 .558 m .36772 .56175 L s P p .001 w .39418 .558 m .39418 .56175 L s P p .001 w .42063 .558 m .42063 .56175 L s P p .001 w .47354 .558 m .47354 .56175 L s P p .001 w .5 .558 m .5 .56175 L s P p .001 w .52646 .558 m .52646 .56175 L s P p .001 w .55291 .558 m .55291 .56175 L s P p .001 w .60582 .558 m .60582 .56175 L s P p .001 w .63228 .558 m .63228 .56175 L s P p .001 w .65873 .558 m .65873 .56175 L s P p .001 w .68519 .558 m .68519 .56175 L s P p .001 w .7381 .558 m .7381 .56175 L s P p .001 w .76455 .558 m .76455 .56175 L s P p .001 w .79101 .558 m .79101 .56175 L s P p .001 w .81746 .558 m .81746 .56175 L s P p .001 w .87037 .558 m .87037 .56175 L s P p .001 w .89683 .558 m .89683 .56175 L s P p .001 w .92328 .558 m .92328 .56175 L s P p .001 w .94974 .558 m .94974 .56175 L s P p .001 w .02381 .558 m .02381 .56175 L s P p .002 w 0 .558 m 1 .558 L s P p .002 w .18254 .02432 m .18879 .02432 L s P [(-4)] .17004 .02432 1 0 Mshowa p .002 w .18254 .15774 m .18879 .15774 L s P [(-3)] .17004 .15774 1 0 Mshowa p .002 w .18254 .29116 m .18879 .29116 L s P [(-2)] .17004 .29116 1 0 Mshowa p .002 w .18254 .42458 m .18879 .42458 L s P [(-1)] .17004 .42458 1 0 Mshowa p .001 w .18254 .051 m .18629 .051 L s P p .001 w .18254 .07768 m .18629 .07768 L s P p .001 w .18254 .10437 m .18629 .10437 L s P p .001 w .18254 .13105 m .18629 .13105 L s P p .001 w .18254 .18442 m .18629 .18442 L s P p .001 w .18254 .21111 m .18629 .21111 L s P p .001 w .18254 .23779 m .18629 .23779 L s P p .001 w .18254 .26447 m .18629 .26447 L s P p .001 w .18254 .31784 m .18629 .31784 L s P p .001 w .18254 .34453 m .18629 .34453 L s P p .001 w .18254 .37121 m .18629 .37121 L s P p .001 w .18254 .3979 m .18629 .3979 L s P p .001 w .18254 .45127 m .18629 .45127 L s P p .001 w .18254 .47795 m .18629 .47795 L s P p .001 w .18254 .50463 m .18629 .50463 L s P p .001 w .18254 .53132 m .18629 .53132 L s P p .001 w .18254 .58469 m .18629 .58469 L s P p .001 w .18254 .61137 m .18629 .61137 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .10466 L .04365 .17409 L .05357 .2288 L .06349 .27265 L .07341 .30835 L .08333 .33779 L .10317 .38305 L .12302 .41575 L .14286 .44014 L .1627 .45881 L .18254 .4734 L .20238 .48503 L .22222 .49444 L .24206 .50216 L .2619 .50857 L .30159 .51851 L .34127 .52576 L .38095 .53121 L .42063 .5354 L .46032 .53869 L .5 .54132 L .53968 .54346 L .57937 .54521 L .61905 .54668 L .65873 .5479 L .69841 .54894 L .7381 .54983 L .77778 .5506 L .81746 .55127 L .85714 .55185 L .89683 .55236 L .93651 .55281 L .97619 .55321 L s P P p 1 0 0 r p .02 w .23535 .53417 Mdot .66225 .60008 Mdot .61379 .5417 Mdot .33114 .49119 Mdot .41729 .56722 Mdot .05434 .23037 Mdot .04695 .16973 Mdot .62614 .57587 Mdot .87534 .56339 Mdot .58312 .57523 Mdot .25548 .4651 Mdot .36925 .57662 Mdot .04394 .15172 Mdot .6832 .59671 Mdot .40634 .56924 Mdot .37587 .55796 Mdot .96524 .50276 Mdot .34534 .50122 Mdot .9421 .57843 Mdot .6954 .54811 Mdot .09057 .31614 Mdot .13609 .45472 Mdot .81127 .49982 Mdot .74303 .57887 Mdot .79115 .57549 Mdot .08364 .37184 Mdot .40862 .57675 Mdot .68396 .50169 Mdot .41957 .57472 Mdot .36243 .56016 Mdot .45367 .49862 Mdot .64323 .52614 Mdot .38691 .5848 Mdot .53095 .55233 Mdot .55182 .59488 Mdot .76411 .58405 Mdot .73687 .52259 Mdot .70428 .602 Mdot .35206 .51073 Mdot .04413 .22198 Mdot .90868 .54975 Mdot .65789 .50932 Mdot .47882 .52318 Mdot .03847 .17449 Mdot .11572 .38783 Mdot .48372 .51345 Mdot .54009 .52936 Mdot .6958 .54164 Mdot .77941 .51032 Mdot .96771 .54628 Mdot .45116 .56742 Mdot .94739 .553 Mdot .51867 .52114 Mdot .31331 .50587 Mdot .06366 .32108 Mdot .29864 .52011 Mdot .92413 .56452 Mdot .79111 .52658 Mdot .40784 .5065 Mdot .11913 .44577 Mdot .60462 .50544 Mdot .12761 .41115 Mdot .17727 .45258 Mdot .15641 .49586 Mdot .63479 .49797 Mdot .8193 .55558 Mdot .59494 .50105 Mdot .54447 .49207 Mdot .647 .54426 Mdot .72954 .57561 Mdot .26297 .4795 Mdot .63423 .59015 Mdot .63453 .50522 Mdot .53043 .54284 Mdot .48108 .56834 Mdot .39782 .54332 Mdot .82248 .57535 Mdot .55472 .49668 Mdot .25136 .52037 Mdot .03406 .05629 Mdot .58055 .51033 Mdot .28071 .4915 Mdot .34139 .57321 Mdot .62267 .53503 Mdot .68304 .58459 Mdot .11606 .44646 Mdot .22578 .45 Mdot .69442 .52672 Mdot .37949 .51452 Mdot .16352 .43096 Mdot .15194 .47046 Mdot .15326 .41784 Mdot .54759 .54683 Mdot .84874 .59477 Mdot .23001 .50616 Mdot .24988 .50592 Mdot .52315 .56856 Mdot .15763 .46702 Mdot .32118 .54043 Mdot .4394 .51856 Mdot .91789 .53378 Mdot .2997 .4742 Mdot .78976 .56365 Mdot .17024 .5062 Mdot .26598 .46589 Mdot .29769 .56956 Mdot .05978 .25707 Mdot .07162 .29872 Mdot .51282 .57014 Mdot .89017 .58992 Mdot .21545 .44764 Mdot .47458 .54221 Mdot .27375 .54029 Mdot .1987 .47561 Mdot .46018 .5009 Mdot .05226 .22536 Mdot .21801 .44484 Mdot .73077 .55523 Mdot .18204 .47766 Mdot .68296 .50336 Mdot .64541 .57756 Mdot .76898 .52091 Mdot .45377 .55826 Mdot .31821 .54285 Mdot .20383 .53033 Mdot .14332 .41505 Mdot .71984 .57852 Mdot .11487 .42674 Mdot .52563 .56582 Mdot .36029 .49324 Mdot .3674 .49414 Mdot .65352 .50426 Mdot .69818 .53688 Mdot .86073 .59124 Mdot .26822 .53128 Mdot .56633 .50932 Mdot .0882 .37913 Mdot .44682 .58068 Mdot .34456 .47313 Mdot .35577 .49363 Mdot .79511 .52806 Mdot .97167 .50076 Mdot .45152 .4969 Mdot .34196 .57006 Mdot .72953 .57392 Mdot .45742 .48994 Mdot .48512 .49006 Mdot .86727 .59212 Mdot .42072 .50878 Mdot .44426 .58767 Mdot .09998 .34907 Mdot .11231 .42976 Mdot .28105 .56274 Mdot .11683 .43601 Mdot .80572 .58782 Mdot .75106 .58904 Mdot .1 .33511 Mdot .31746 .55547 Mdot .59108 .50146 Mdot .42637 .58294 Mdot .19417 .51585 Mdot .9583 .60314 Mdot .118 .41699 Mdot .8698 .51764 Mdot .81314 .56451 Mdot .77678 .53949 Mdot .03122 .11227 Mdot .04953 .20994 Mdot .90741 .56988 Mdot .70826 .57174 Mdot .34014 .53355 Mdot .3057 .54796 Mdot .16979 .48979 Mdot .32359 .5546 Mdot .0756 .27752 Mdot .42998 .54904 Mdot .23865 .50379 Mdot .62939 .51751 Mdot .23124 .52441 Mdot .60367 .56666 Mdot .30002 .47369 Mdot .8716 .60332 Mdot .93606 .55368 Mdot .58971 .51464 Mdot .79009 .57922 Mdot .28992 .50614 Mdot .7383 .56167 Mdot .83613 .58235 Mdot .52345 .49788 Mdot .23056 .50444 Mdot .31602 .50331 Mdot .60308 .58563 Mdot .03982 .12069 Mdot .70767 .59071 Mdot .07994 .35282 Mdot .14177 .45755 Mdot .26605 .5598 Mdot .82804 .52663 Mdot .50394 .52642 Mdot .9681 .55083 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{1.42653, -0.00352407, 0, -0.00352407} :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 2.565137958227423, A -> 0.996475932308546, R -> 100} ;[o] {a -> 2.56514, A -> 0.996476, R -> 100} :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] angle fixed, fit for Hamaker and radius, the goal is A->1, R->100: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={a->1.05712}; sol=Sort[Join[NonlinearFit[dat2,ptmodel/.cond,d,{{R,50,150},{A,0.5,1.5}}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp2,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-100)/100,(A R-100)/100}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 299; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.557437 0.135087 [ [(2.5)] .05026 .55744 0 2 Msboxa [(7.5)] .31481 .55744 0 2 Msboxa [(10)] .44709 .55744 0 2 Msboxa [(12.5)] .57937 .55744 0 2 Msboxa [(15)] .71164 .55744 0 2 Msboxa [(17.5)] .84392 .55744 0 2 Msboxa [(20)] .97619 .55744 0 2 Msboxa [(-4)] .17004 .01709 1 0 Msboxa [(-3)] .17004 .15218 1 0 Msboxa [(-2)] .17004 .28726 1 0 Msboxa [(-1)] .17004 .42235 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .55744 m .05026 .56369 L s P [(2.5)] .05026 .55744 0 2 Mshowa p .002 w .31481 .55744 m .31481 .56369 L s P [(7.5)] .31481 .55744 0 2 Mshowa p .002 w .44709 .55744 m .44709 .56369 L s P [(10)] .44709 .55744 0 2 Mshowa p .002 w .57937 .55744 m .57937 .56369 L s P [(12.5)] .57937 .55744 0 2 Mshowa p .002 w .71164 .55744 m .71164 .56369 L s P [(15)] .71164 .55744 0 2 Mshowa p .002 w .84392 .55744 m .84392 .56369 L s P [(17.5)] .84392 .55744 0 2 Mshowa p .002 w .97619 .55744 m .97619 .56369 L s P [(20)] .97619 .55744 0 2 Mshowa p .001 w .07672 .55744 m .07672 .56119 L s P p .001 w .10317 .55744 m .10317 .56119 L s P p .001 w .12963 .55744 m .12963 .56119 L s P p .001 w .15608 .55744 m .15608 .56119 L s P p .001 w .20899 .55744 m .20899 .56119 L s P p .001 w .23545 .55744 m .23545 .56119 L s P p .001 w .2619 .55744 m .2619 .56119 L s P p .001 w .28836 .55744 m .28836 .56119 L s P p .001 w .34127 .55744 m .34127 .56119 L s P p .001 w .36772 .55744 m .36772 .56119 L s P p .001 w .39418 .55744 m .39418 .56119 L s P p .001 w .42063 .55744 m .42063 .56119 L s P p .001 w .47354 .55744 m .47354 .56119 L s P p .001 w .5 .55744 m .5 .56119 L s P p .001 w .52646 .55744 m .52646 .56119 L s P p .001 w .55291 .55744 m .55291 .56119 L s P p .001 w .60582 .55744 m .60582 .56119 L s P p .001 w .63228 .55744 m .63228 .56119 L s P p .001 w .65873 .55744 m .65873 .56119 L s P p .001 w .68519 .55744 m .68519 .56119 L s P p .001 w .7381 .55744 m .7381 .56119 L s P p .001 w .76455 .55744 m .76455 .56119 L s P p .001 w .79101 .55744 m .79101 .56119 L s P p .001 w .81746 .55744 m .81746 .56119 L s P p .001 w .87037 .55744 m .87037 .56119 L s P p .001 w .89683 .55744 m .89683 .56119 L s P p .001 w .92328 .55744 m .92328 .56119 L s P p .001 w .94974 .55744 m .94974 .56119 L s P p .001 w .02381 .55744 m .02381 .56119 L s P p .002 w 0 .55744 m 1 .55744 L s P p .002 w .18254 .01709 m .18879 .01709 L s P [(-4)] .17004 .01709 1 0 Mshowa p .002 w .18254 .15218 m .18879 .15218 L s P [(-3)] .17004 .15218 1 0 Mshowa p .002 w .18254 .28726 m .18879 .28726 L s P [(-2)] .17004 .28726 1 0 Mshowa p .002 w .18254 .42235 m .18879 .42235 L s P [(-1)] .17004 .42235 1 0 Mshowa p .001 w .18254 .04411 m .18629 .04411 L s P p .001 w .18254 .07112 m .18629 .07112 L s P p .001 w .18254 .09814 m .18629 .09814 L s P p .001 w .18254 .12516 m .18629 .12516 L s P p .001 w .18254 .17919 m .18629 .17919 L s P p .001 w .18254 .20621 m .18629 .20621 L s P p .001 w .18254 .23323 m .18629 .23323 L s P p .001 w .18254 .26025 m .18629 .26025 L s P p .001 w .18254 .31428 m .18629 .31428 L s P p .001 w .18254 .3413 m .18629 .3413 L s P p .001 w .18254 .36832 m .18629 .36832 L s P p .001 w .18254 .39533 m .18629 .39533 L s P p .001 w .18254 .44937 m .18629 .44937 L s P p .001 w .18254 .47638 m .18629 .47638 L s P p .001 w .18254 .5034 m .18629 .5034 L s P p .001 w .18254 .53042 m .18629 .53042 L s P p .001 w .18254 .58445 m .18629 .58445 L s P p .001 w .18254 .61147 m .18629 .61147 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .10494 L .04365 .17455 L .05357 .22936 L .06349 .27329 L .07341 .30902 L .08333 .33847 L .10317 .38371 L .12302 .41637 L .14286 .44069 L .1627 .45929 L .18254 .47382 L .20238 .48539 L .22222 .49474 L .24206 .5024 L .2619 .50876 L .30159 .5186 L .34127 .52577 L .38095 .53114 L .42063 .53527 L .46032 .53851 L .5 .5411 L .53968 .54319 L .57937 .54492 L .61905 .54635 L .65873 .54755 L .69841 .54857 L .7381 .54944 L .77778 .55018 L .81746 .55083 L .85714 .5514 L .89683 .5519 L .93651 .55233 L .97619 .55272 L s P P p 1 0 0 r p .02 w .23535 .5333 Mdot .66225 .60004 Mdot .61379 .54093 Mdot .33114 .48979 Mdot .41729 .56677 Mdot .05434 .22572 Mdot .04695 .16432 Mdot .62614 .57552 Mdot .87534 .5629 Mdot .58312 .57488 Mdot .25548 .46337 Mdot .36925 .57628 Mdot .04394 .14608 Mdot .6832 .59663 Mdot .40634 .56881 Mdot .37587 .55739 Mdot .96524 .50151 Mdot .34534 .49995 Mdot .9421 .57812 Mdot .6954 .54743 Mdot .09057 .31255 Mdot .13609 .45286 Mdot .81127 .49853 Mdot .74303 .57856 Mdot .79115 .57514 Mdot .08364 .36895 Mdot .40862 .57641 Mdot .68396 .50042 Mdot .41957 .57436 Mdot .36243 .55962 Mdot .45367 .49732 Mdot .64323 .52518 Mdot .38691 .58457 Mdot .53095 .5517 Mdot .55182 .59478 Mdot .76411 .58381 Mdot .73687 .52158 Mdot .70428 .60198 Mdot .35206 .50957 Mdot .04413 .21722 Mdot .90868 .54908 Mdot .65789 .50815 Mdot .47882 .52218 Mdot .03847 .16914 Mdot .11572 .38514 Mdot .48372 .51232 Mdot .54009 .52843 Mdot .6958 .54088 Mdot .77941 .50915 Mdot .96771 .54557 Mdot .45116 .56697 Mdot .94739 .55237 Mdot .51867 .52011 Mdot .31331 .50465 Mdot .06366 .31756 Mdot .29864 .51907 Mdot .92413 .56403 Mdot .79111 .52562 Mdot .40784 .50529 Mdot .11913 .4438 Mdot .60462 .50422 Mdot .12761 .40875 Mdot .17727 .4507 Mdot .15641 .49452 Mdot .63479 .49666 Mdot .8193 .55499 Mdot .59494 .49977 Mdot .54447 .49068 Mdot .647 .54352 Mdot .72954 .57527 Mdot .26297 .47796 Mdot .63423 .58999 Mdot .63453 .50399 Mdot .53043 .54209 Mdot .48108 .5679 Mdot .39782 .54257 Mdot .82248 .575 Mdot .55472 .49535 Mdot .25136 .51934 Mdot .03406 .04947 Mdot .58055 .50917 Mdot .28071 .49011 Mdot .34139 .57283 Mdot .62267 .53418 Mdot .68304 .58435 Mdot .11606 .4445 Mdot .22578 .44809 Mdot .69442 .52577 Mdot .37949 .51341 Mdot .16352 .42881 Mdot .15194 .4688 Mdot .15326 .41553 Mdot .54759 .54613 Mdot .84874 .59466 Mdot .23001 .50495 Mdot .24988 .5047 Mdot .52315 .56813 Mdot .15763 .46532 Mdot .32118 .53965 Mdot .4394 .5175 Mdot .91789 .53291 Mdot .2997 .47259 Mdot .78976 .56315 Mdot .17024 .50499 Mdot .26598 .46418 Mdot .29769 .56914 Mdot .05978 .25275 Mdot .07162 .29492 Mdot .51282 .56973 Mdot .89017 .58975 Mdot .21545 .44569 Mdot .47458 .54144 Mdot .27375 .5395 Mdot .1987 .47402 Mdot .46018 .49962 Mdot .05226 .22065 Mdot .21801 .44286 Mdot .73077 .55463 Mdot .18204 .4761 Mdot .68296 .50211 Mdot .64541 .57723 Mdot .76898 .51989 Mdot .45377 .5577 Mdot .31821 .54209 Mdot .20383 .52942 Mdot .14332 .4127 Mdot .71984 .57821 Mdot .11487 .42454 Mdot .52563 .56536 Mdot .36029 .49186 Mdot .3674 .49277 Mdot .65352 .50302 Mdot .69818 .53605 Mdot .86073 .59109 Mdot .26822 .53038 Mdot .56633 .50815 Mdot .0882 .37633 Mdot .44682 .5804 Mdot .34456 .47151 Mdot .35577 .49226 Mdot .79511 .52712 Mdot .97167 .49948 Mdot .45152 .49557 Mdot .34196 .56965 Mdot .72953 .57355 Mdot .45742 .48852 Mdot .48512 .48864 Mdot .86727 .59198 Mdot .42072 .5076 Mdot .44426 .58747 Mdot .09998 .34589 Mdot .11231 .42759 Mdot .28105 .56224 Mdot .11683 .43392 Mdot .80572 .58763 Mdot .75106 .58886 Mdot .1 .33177 Mdot .31746 .55487 Mdot .59108 .50019 Mdot .42637 .58268 Mdot .19417 .51476 Mdot .9583 .60314 Mdot .118 .41467 Mdot .8698 .51657 Mdot .81314 .56402 Mdot .77678 .53869 Mdot .03122 .10615 Mdot .04953 .20503 Mdot .90741 .56947 Mdot .70826 .57135 Mdot .34014 .53268 Mdot .3057 .54727 Mdot .16979 .48837 Mdot .32359 .55399 Mdot .0756 .27345 Mdot .42998 .54836 Mdot .23865 .50255 Mdot .62939 .51644 Mdot .23124 .52342 Mdot .60367 .56621 Mdot .30002 .47207 Mdot .8716 .60332 Mdot .93606 .55306 Mdot .58971 .51354 Mdot .79009 .57892 Mdot .28992 .50493 Mdot .7383 .56115 Mdot .83613 .58209 Mdot .52345 .49656 Mdot .23056 .5032 Mdot .31602 .50206 Mdot .60308 .58541 Mdot .03982 .11467 Mdot .70767 .59055 Mdot .07994 .34969 Mdot .14177 .45573 Mdot .26605 .55925 Mdot .82804 .52567 Mdot .50394 .52546 Mdot .9681 .55017 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{0., 0.581237, -0.371767, -0.00661456} :[font = output; inactive; formatted; output; nowordwrap; ] {a -> 1.05712, A -> 1.581237315822473, R -> 62.82329852043971} ;[o] {a -> 1.05712, A -> 1.58124, R -> 62.8233} :[font = input; startGroup; nowordwrap; ] A R/.% :[font = output; inactive; formatted; output; endGroup; endGroup; nowordwrap; ] 99.338543923574 ;[o] 99.3385 :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] now everything free: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={}; sol=Sort[Join[NonlinearFit[dat2,ptmodel/.cond,d,{a,{R,50,150},{A,0.5,1.5}}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp2,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-100)/100,(A R-100)/100}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 299; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.557409 0.135169 [ [(2.5)] .05026 .55741 0 2 Msboxa [(7.5)] .31481 .55741 0 2 Msboxa [(10)] .44709 .55741 0 2 Msboxa [(12.5)] .57937 .55741 0 2 Msboxa [(15)] .71164 .55741 0 2 Msboxa [(17.5)] .84392 .55741 0 2 Msboxa [(20)] .97619 .55741 0 2 Msboxa [(-4)] .17004 .01673 1 0 Msboxa [(-3)] .17004 .1519 1 0 Msboxa [(-2)] .17004 .28707 1 0 Msboxa [(-1)] .17004 .42224 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .55741 m .05026 .56366 L s P [(2.5)] .05026 .55741 0 2 Mshowa p .002 w .31481 .55741 m .31481 .56366 L s P [(7.5)] .31481 .55741 0 2 Mshowa p .002 w .44709 .55741 m .44709 .56366 L s P [(10)] .44709 .55741 0 2 Mshowa p .002 w .57937 .55741 m .57937 .56366 L s P [(12.5)] .57937 .55741 0 2 Mshowa p .002 w .71164 .55741 m .71164 .56366 L s P [(15)] .71164 .55741 0 2 Mshowa p .002 w .84392 .55741 m .84392 .56366 L s P [(17.5)] .84392 .55741 0 2 Mshowa p .002 w .97619 .55741 m .97619 .56366 L s P [(20)] .97619 .55741 0 2 Mshowa p .001 w .07672 .55741 m .07672 .56116 L s P p .001 w .10317 .55741 m .10317 .56116 L s P p .001 w .12963 .55741 m .12963 .56116 L s P p .001 w .15608 .55741 m .15608 .56116 L s P p .001 w .20899 .55741 m .20899 .56116 L s P p .001 w .23545 .55741 m .23545 .56116 L s P p .001 w .2619 .55741 m .2619 .56116 L s P p .001 w .28836 .55741 m .28836 .56116 L s P p .001 w .34127 .55741 m .34127 .56116 L s P p .001 w .36772 .55741 m .36772 .56116 L s P p .001 w .39418 .55741 m .39418 .56116 L s P p .001 w .42063 .55741 m .42063 .56116 L s P p .001 w .47354 .55741 m .47354 .56116 L s P p .001 w .5 .55741 m .5 .56116 L s P p .001 w .52646 .55741 m .52646 .56116 L s P p .001 w .55291 .55741 m .55291 .56116 L s P p .001 w .60582 .55741 m .60582 .56116 L s P p .001 w .63228 .55741 m .63228 .56116 L s P p .001 w .65873 .55741 m .65873 .56116 L s P p .001 w .68519 .55741 m .68519 .56116 L s P p .001 w .7381 .55741 m .7381 .56116 L s P p .001 w .76455 .55741 m .76455 .56116 L s P p .001 w .79101 .55741 m .79101 .56116 L s P p .001 w .81746 .55741 m .81746 .56116 L s P p .001 w .87037 .55741 m .87037 .56116 L s P p .001 w .89683 .55741 m .89683 .56116 L s P p .001 w .92328 .55741 m .92328 .56116 L s P p .001 w .94974 .55741 m .94974 .56116 L s P p .001 w .02381 .55741 m .02381 .56116 L s P p .002 w 0 .55741 m 1 .55741 L s P p .002 w .18254 .01673 m .18879 .01673 L s P [(-4)] .17004 .01673 1 0 Mshowa p .002 w .18254 .1519 m .18879 .1519 L s P [(-3)] .17004 .1519 1 0 Mshowa p .002 w .18254 .28707 m .18879 .28707 L s P [(-2)] .17004 .28707 1 0 Mshowa p .002 w .18254 .42224 m .18879 .42224 L s P [(-1)] .17004 .42224 1 0 Mshowa p .001 w .18254 .04377 m .18629 .04377 L s P p .001 w .18254 .0708 m .18629 .0708 L s P p .001 w .18254 .09783 m .18629 .09783 L s P p .001 w .18254 .12487 m .18629 .12487 L s P p .001 w .18254 .17894 m .18629 .17894 L s P p .001 w .18254 .20597 m .18629 .20597 L s P p .001 w .18254 .233 m .18629 .233 L s P p .001 w .18254 .26004 m .18629 .26004 L s P p .001 w .18254 .3141 m .18629 .3141 L s P p .001 w .18254 .34114 m .18629 .34114 L s P p .001 w .18254 .36817 m .18629 .36817 L s P p .001 w .18254 .39521 m .18629 .39521 L s P p .001 w .18254 .44927 m .18629 .44927 L s P p .001 w .18254 .47631 m .18629 .47631 L s P p .001 w .18254 .50334 m .18629 .50334 L s P p .001 w .18254 .53038 m .18629 .53038 L s P p .001 w .18254 .58444 m .18629 .58444 L s P p .001 w .18254 .61148 m .18629 .61148 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .10448 L .04365 .17378 L .05357 .22838 L .06349 .27217 L .07341 .30782 L .08333 .33722 L .10317 .38242 L .12302 .4151 L .14286 .43947 L .1627 .45813 L .18254 .47273 L .20238 .48436 L .22222 .49378 L .24206 .5015 L .2619 .50792 L .30159 .51787 L .34127 .52514 L .38095 .5306 L .42063 .5348 L .46032 .53811 L .5 .54075 L .53968 .54289 L .57937 .54466 L .61905 .54613 L .65873 .54737 L .69841 .54841 L .7381 .54931 L .77778 .55008 L .81746 .55076 L .85714 .55134 L .89683 .55186 L .93651 .55231 L .97619 .55271 L s P P p 1 0 0 r p .02 w .23535 .53326 Mdot .66225 .60004 Mdot .61379 .54089 Mdot .33114 .48972 Mdot .41729 .56675 Mdot .05434 .22549 Mdot .04695 .16406 Mdot .62614 .57551 Mdot .87534 .56287 Mdot .58312 .57486 Mdot .25548 .46329 Mdot .36925 .57627 Mdot .04394 .1458 Mdot .6832 .59663 Mdot .40634 .56879 Mdot .37587 .55736 Mdot .96524 .50144 Mdot .34534 .49989 Mdot .9421 .5781 Mdot .6954 .54739 Mdot .09057 .31238 Mdot .13609 .45277 Mdot .81127 .49847 Mdot .74303 .57855 Mdot .79115 .57512 Mdot .08364 .36881 Mdot .40862 .5764 Mdot .68396 .50035 Mdot .41957 .57434 Mdot .36243 .5596 Mdot .45367 .49725 Mdot .64323 .52513 Mdot .38691 .58456 Mdot .53095 .55167 Mdot .55182 .59477 Mdot .76411 .5838 Mdot .73687 .52153 Mdot .70428 .60198 Mdot .35206 .50952 Mdot .04413 .21698 Mdot .90868 .54905 Mdot .65789 .50809 Mdot .47882 .52213 Mdot .03847 .16888 Mdot .11572 .38501 Mdot .48372 .51227 Mdot .54009 .52839 Mdot .6958 .54084 Mdot .77941 .5091 Mdot .96771 .54553 Mdot .45116 .56695 Mdot .94739 .55234 Mdot .51867 .52006 Mdot .31331 .50459 Mdot .06366 .31738 Mdot .29864 .51902 Mdot .92413 .56401 Mdot .79111 .52557 Mdot .40784 .50523 Mdot .11913 .44371 Mdot .60462 .50416 Mdot .12761 .40863 Mdot .17727 .45061 Mdot .15641 .49446 Mdot .63479 .49659 Mdot .8193 .55496 Mdot .59494 .49971 Mdot .54447 .49061 Mdot .647 .54348 Mdot .72954 .57525 Mdot .26297 .47788 Mdot .63423 .58998 Mdot .63453 .50393 Mdot .53043 .54205 Mdot .48108 .56788 Mdot .39782 .54253 Mdot .82248 .57498 Mdot .55472 .49528 Mdot .25136 .51928 Mdot .03406 .04913 Mdot .58055 .50911 Mdot .28071 .49004 Mdot .34139 .57282 Mdot .62267 .53414 Mdot .68304 .58434 Mdot .11606 .44441 Mdot .22578 .448 Mdot .69442 .52572 Mdot .37949 .51336 Mdot .16352 .42871 Mdot .15194 .46872 Mdot .15326 .41541 Mdot .54759 .54609 Mdot .84874 .59466 Mdot .23001 .50489 Mdot .24988 .50464 Mdot .52315 .56811 Mdot .15763 .46523 Mdot .32118 .53961 Mdot .4394 .51745 Mdot .91789 .53287 Mdot .2997 .47251 Mdot .78976 .56313 Mdot .17024 .50493 Mdot .26598 .46409 Mdot .29769 .56912 Mdot .05978 .25254 Mdot .07162 .29473 Mdot .51282 .56971 Mdot .89017 .58974 Mdot .21545 .4456 Mdot .47458 .54141 Mdot .27375 .53947 Mdot .1987 .47394 Mdot .46018 .49956 Mdot .05226 .22041 Mdot .21801 .44276 Mdot .73077 .5546 Mdot .18204 .47602 Mdot .68296 .50205 Mdot .64541 .57722 Mdot .76898 .51983 Mdot .45377 .55767 Mdot .31821 .54205 Mdot .20383 .52938 Mdot .14332 .41258 Mdot .71984 .57819 Mdot .11487 .42443 Mdot .52563 .56533 Mdot .36029 .49179 Mdot .3674 .49271 Mdot .65352 .50296 Mdot .69818 .53601 Mdot .86073 .59108 Mdot .26822 .53034 Mdot .56633 .50809 Mdot .0882 .3762 Mdot .44682 .58038 Mdot .34456 .47143 Mdot .35577 .49219 Mdot .79511 .52707 Mdot .97167 .49941 Mdot .45152 .4955 Mdot .34196 .56963 Mdot .72953 .57353 Mdot .45742 .48845 Mdot .48512 .48857 Mdot .86727 .59197 Mdot .42072 .50754 Mdot .44426 .58746 Mdot .09998 .34574 Mdot .11231 .42749 Mdot .28105 .56221 Mdot .11683 .43382 Mdot .80572 .58762 Mdot .75106 .58885 Mdot .1 .3316 Mdot .31746 .55484 Mdot .59108 .50013 Mdot .42637 .58267 Mdot .19417 .51471 Mdot .9583 .60314 Mdot .118 .41455 Mdot .8698 .51652 Mdot .81314 .564 Mdot .77678 .53865 Mdot .03122 .10584 Mdot .04953 .20479 Mdot .90741 .56944 Mdot .70826 .57133 Mdot .34014 .53263 Mdot .3057 .54723 Mdot .16979 .4883 Mdot .32359 .55396 Mdot .0756 .27325 Mdot .42998 .54832 Mdot .23865 .50248 Mdot .62939 .51638 Mdot .23124 .52337 Mdot .60367 .56618 Mdot .30002 .47199 Mdot .8716 .60332 Mdot .93606 .55303 Mdot .58971 .51348 Mdot .79009 .57891 Mdot .28992 .50487 Mdot .7383 .56112 Mdot .83613 .58208 Mdot .52345 .4965 Mdot .23056 .50314 Mdot .31602 .502 Mdot .60308 .58539 Mdot .03982 .11437 Mdot .70767 .59054 Mdot .07994 .34954 Mdot .14177 .45564 Mdot .26605 .55923 Mdot .82804 .52562 Mdot .50394 .52541 Mdot .9681 .55014 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{0.792229, -0.163228, 0.171261, -0.0199214} :[font = output; inactive; formatted; output; nowordwrap; ] {a -> 1.894601451192916, A -> 0.836772393787923, R -> 117.1260716927685} ;[o] {a -> 1.8946, A -> 0.836772, R -> 117.126} :[font = input; startGroup; nowordwrap; ] R A/.% :[font = output; inactive; formatted; output; endGroup; endGroup; nowordwrap; ] 98.0078633853338 ;[o] 98.0079 :[font = subsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 12; ] dat3 is the smooth data, small radius: :[font = input; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] dat3=Table[Module[{temp}, temp=Random[Real,{2,20}]; {temp,noise*0+ForceParametricTip[1.05712,temp,1,5]}], {x,1,200}]; :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] lp3=ListPlot[dat3,PlotRange->All,PlotStyle->{Red,PointSize[0.02]}] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0890576 0.0537647 0.61689 3.83475 [ [(2.5)] .04535 .61689 0 2 Msboxa [(7.5)] .31418 .61689 0 2 Msboxa [(10)] .44859 .61689 0 2 Msboxa [(12.5)] .583 .61689 0 2 Msboxa [(15)] .71741 .61689 0 2 Msboxa [(17.5)] .85182 .61689 0 2 Msboxa [(20)] .98624 .61689 0 2 Msboxa [(-0.15)] .16727 .04168 1 0 Msboxa [(-0.125)] .16727 .13755 1 0 Msboxa [(-0.1)] .16727 .23342 1 0 Msboxa [(-0.075)] .16727 .32928 1 0 Msboxa [(-0.05)] .16727 .42515 1 0 Msboxa [(-0.025)] .16727 .52102 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .62414 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .04535 .61689 m .04535 .62314 L s P [(2.5)] .04535 .61689 0 2 Mshowa p .002 w .31418 .61689 m .31418 .62314 L s P [(7.5)] .31418 .61689 0 2 Mshowa p .002 w .44859 .61689 m .44859 .62314 L s P [(10)] .44859 .61689 0 2 Mshowa p .002 w .583 .61689 m .583 .62314 L s P [(12.5)] .583 .61689 0 2 Mshowa p .002 w .71741 .61689 m .71741 .62314 L s P [(15)] .71741 .61689 0 2 Mshowa p .002 w .85182 .61689 m .85182 .62314 L s P [(17.5)] .85182 .61689 0 2 Mshowa p .002 w .98624 .61689 m .98624 .62314 L s P [(20)] .98624 .61689 0 2 Mshowa p .001 w .07224 .61689 m .07224 .62064 L s P p .001 w .09912 .61689 m .09912 .62064 L s P p .001 w .126 .61689 m .126 .62064 L s P p .001 w .15288 .61689 m .15288 .62064 L s P p .001 w .20665 .61689 m .20665 .62064 L s P p .001 w .23353 .61689 m .23353 .62064 L s P p .001 w .26041 .61689 m .26041 .62064 L s P p .001 w .2873 .61689 m .2873 .62064 L s P p .001 w .34106 .61689 m .34106 .62064 L s P p .001 w .36794 .61689 m .36794 .62064 L s P p .001 w .39482 .61689 m .39482 .62064 L s P p .001 w .42171 .61689 m .42171 .62064 L s P p .001 w .47547 .61689 m .47547 .62064 L s P p .001 w .50235 .61689 m .50235 .62064 L s P p .001 w .52924 .61689 m .52924 .62064 L s P p .001 w .55612 .61689 m .55612 .62064 L s P p .001 w .60988 .61689 m .60988 .62064 L s P p .001 w .63677 .61689 m .63677 .62064 L s P p .001 w .66365 .61689 m .66365 .62064 L s P p .001 w .69053 .61689 m .69053 .62064 L s P p .001 w .7443 .61689 m .7443 .62064 L s P p .001 w .77118 .61689 m .77118 .62064 L s P p .001 w .79806 .61689 m .79806 .62064 L s P p .001 w .82494 .61689 m .82494 .62064 L s P p .001 w .87871 .61689 m .87871 .62064 L s P p .001 w .90559 .61689 m .90559 .62064 L s P p .001 w .93247 .61689 m .93247 .62064 L s P p .001 w .95935 .61689 m .95935 .62064 L s P p .001 w .01847 .61689 m .01847 .62064 L s P p .002 w 0 .61689 m 1 .61689 L s P p .002 w .17977 .04168 m .18602 .04168 L s P [(-0.15)] .16727 .04168 1 0 Mshowa p .002 w .17977 .13755 m .18602 .13755 L s P [(-0.125)] .16727 .13755 1 0 Mshowa p .002 w .17977 .23342 m .18602 .23342 L s P [(-0.1)] .16727 .23342 1 0 Mshowa p .002 w .17977 .32928 m .18602 .32928 L s P [(-0.075)] .16727 .32928 1 0 Mshowa p .002 w .17977 .42515 m .18602 .42515 L s P [(-0.05)] .16727 .42515 1 0 Mshowa p .002 w .17977 .52102 m .18602 .52102 L s P [(-0.025)] .16727 .52102 1 0 Mshowa p .001 w .17977 .06085 m .18352 .06085 L s P p .001 w .17977 .08003 m .18352 .08003 L s P p .001 w .17977 .0992 m .18352 .0992 L s P p .001 w .17977 .11837 m .18352 .11837 L s P p .001 w .17977 .15672 m .18352 .15672 L s P p .001 w .17977 .17589 m .18352 .17589 L s P p .001 w .17977 .19507 m .18352 .19507 L s P p .001 w .17977 .21424 m .18352 .21424 L s P p .001 w .17977 .25259 m .18352 .25259 L s P p .001 w .17977 .27176 m .18352 .27176 L s P p .001 w .17977 .29094 m .18352 .29094 L s P p .001 w .17977 .31011 m .18352 .31011 L s P p .001 w .17977 .34846 m .18352 .34846 L s P p .001 w .17977 .36763 m .18352 .36763 L s P p .001 w .17977 .38681 m .18352 .38681 L s P p .001 w .17977 .40598 m .18352 .40598 L s P p .001 w .17977 .44433 m .18352 .44433 L s P p .001 w .17977 .4635 m .18352 .4635 L s P p .001 w .17977 .48267 m .18352 .48267 L s P p .001 w .17977 .50185 m .18352 .50185 L s P p .001 w .17977 .5402 m .18352 .5402 L s P p .001 w .17977 .55937 m .18352 .55937 L s P p .001 w .17977 .57854 m .18352 .57854 L s P p .001 w .17977 .59772 m .18352 .59772 L s P p .001 w .17977 .0225 m .18352 .0225 L s P p .001 w .17977 .00333 m .18352 .00333 L s P p .002 w .17977 0 m .17977 .61803 L s P P p 1 0 0 r p .02 w .88022 .60155 Mdot .68302 .59605 Mdot .64887 .5947 Mdot .3924 .57513 Mdot .6534 .59489 Mdot .36612 .57126 Mdot .47355 .58401 Mdot .48776 .58521 Mdot .86081 .60114 Mdot .78688 .59935 Mdot .60114 .5925 Mdot .47142 .58382 Mdot .70715 .59692 Mdot .77464 .59901 Mdot .07675 .3384 Mdot .40071 .57623 Mdot .40959 .57735 Mdot .05305 .23884 Mdot .92227 .60238 Mdot .55153 .58974 Mdot .07994 .34861 Mdot .75089 .59832 Mdot .46504 .58324 Mdot .29793 .55751 Mdot .74412 .59812 Mdot .50953 .58691 Mdot .68584 .59616 Mdot .1273 .45007 Mdot .29473 .55668 Mdot .09271 .38412 Mdot .35869 .57005 Mdot .52742 .58818 Mdot .29722 .55732 Mdot .76277 .59867 Mdot .81842 .60016 Mdot .48331 .58484 Mdot .09277 .38427 Mdot .96002 .60305 Mdot .39317 .57524 Mdot .85119 .60092 Mdot .11691 .43379 Mdot .77695 .59908 Mdot .74446 .59813 Mdot .268 .54893 Mdot .46571 .5833 Mdot .49146 .58551 Mdot .63353 .59404 Mdot .02662 .04419 Mdot .55923 .59021 Mdot .05284 .23772 Mdot .18453 .50868 Mdot .18789 .51102 Mdot .08609 .36672 Mdot .39718 .57577 Mdot .32787 .56435 Mdot .14765 .47588 Mdot .8484 .60086 Mdot .64242 .59443 Mdot .23334 .53594 Mdot .63427 .59407 Mdot .66035 .59517 Mdot .5999 .59244 Mdot .49429 .58574 Mdot .43048 .57978 Mdot .42667 .57936 Mdot .05178 .23194 Mdot .11727 .4344 Mdot .89036 .60176 Mdot .25511 .54455 Mdot .26641 .54842 Mdot .04024 .1598 Mdot .61838 .59334 Mdot .36612 .57126 Mdot .03696 .13546 Mdot .85807 .60108 Mdot .59272 .59207 Mdot .44987 .5818 Mdot .55941 .59022 Mdot .35107 .56875 Mdot .65529 .59497 Mdot .11443 .42951 Mdot .40734 .57708 Mdot .09266 .384 Mdot .7752 .59903 Mdot .71278 .59711 Mdot .75671 .5985 Mdot .84095 .60069 Mdot .18246 .50719 Mdot .40955 .57735 Mdot .60929 .59291 Mdot .07696 .33909 Mdot .94024 .6027 Mdot .94877 .60286 Mdot .55137 .58973 Mdot .87458 .60143 Mdot .76241 .59866 Mdot .18027 .50557 Mdot .02417 .01864 Mdot .32556 .56387 Mdot .82794 .60039 Mdot .90224 .60199 Mdot .23713 .53756 Mdot .84375 .60076 Mdot .28312 .55352 Mdot .88122 .60157 Mdot .71799 .59729 Mdot .02511 .02866 Mdot .94182 .60273 Mdot .83108 .60046 Mdot .93612 .60263 Mdot .52399 .58794 Mdot .12665 .44912 Mdot .60799 .59284 Mdot .87575 .60146 Mdot .75047 .59831 Mdot .6111 .59299 Mdot .85548 .60102 Mdot .87935 .60153 Mdot .84884 .60087 Mdot .92376 .60241 Mdot .03624 .12983 Mdot .97388 .60328 Mdot .49848 .58607 Mdot .86571 .60124 Mdot .87673 .60148 Mdot .97619 .60332 Mdot .14474 .47261 Mdot .38356 .5739 Mdot .27549 .55127 Mdot .49045 .58543 Mdot .41288 .57775 Mdot .55292 .58983 Mdot .39512 .5755 Mdot .56528 .59056 Mdot .88287 .60161 Mdot .68581 .59616 Mdot .02461 .02333 Mdot .69585 .59652 Mdot .86611 .60125 Mdot .33076 .56493 Mdot .6091 .5929 Mdot .82655 .60035 Mdot .21469 .52705 Mdot .2921 .55599 Mdot .8058 .59984 Mdot .71305 .59712 Mdot .90917 .60213 Mdot .04572 .19645 Mdot .90303 .60201 Mdot .3361 .56599 Mdot .05539 .25095 Mdot .02381 .01472 Mdot .43253 .58 Mdot .1835 .50793 Mdot .23632 .53721 Mdot .87764 .6015 Mdot .41675 .57822 Mdot .18306 .50762 Mdot .49382 .5857 Mdot .15581 .48445 Mdot .57702 .59123 Mdot .80595 .59984 Mdot .5401 .58902 Mdot .80061 .59971 Mdot .12604 .44823 Mdot .63558 .59413 Mdot .87596 .60146 Mdot .74418 .59812 Mdot .47768 .58437 Mdot .57959 .59137 Mdot .9701 .60322 Mdot .44225 .58103 Mdot .41155 .57759 Mdot .62254 .59354 Mdot .85768 .60107 Mdot .80817 .5999 Mdot .75011 .5983 Mdot .19106 .51316 Mdot .86039 .60113 Mdot .43311 .58007 Mdot .40118 .57629 Mdot .83975 .60066 Mdot .41731 .57828 Mdot .41597 .57812 Mdot .02423 .01931 Mdot .77966 .59915 Mdot .15279 .48138 Mdot .95773 .60301 Mdot .38891 .57465 Mdot .78515 .5993 Mdot .51476 .58729 Mdot .3705 .57195 Mdot .13205 .45676 Mdot .51699 .58745 Mdot .70097 .5967 Mdot .11949 .43806 Mdot .69521 .5965 Mdot .32606 .56397 Mdot .5609 .59031 Mdot .35456 .56935 Mdot P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] No Input Form was saved for this expression. ;[o] -Graphics- :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] First test, radius fixed, fit for Hamaker and angle, the goal is A->1, a->1.05712: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={R->5}; sol=Sort[Join[NonlinearFit[dat3,ptmodel/.cond,d,{a,A}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp3,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-5)/5,(A R-5)/5}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.615437 3.46501 [ [(2.5)] .05026 .61544 0 2 Msboxa [(7.5)] .31481 .61544 0 2 Msboxa [(10)] .44709 .61544 0 2 Msboxa [(12.5)] .57937 .61544 0 2 Msboxa [(15)] .71164 .61544 0 2 Msboxa [(17.5)] .84392 .61544 0 2 Msboxa [(20)] .97619 .61544 0 2 Msboxa [(-0.175)] .17004 .00906 1 0 Msboxa [(-0.15)] .17004 .09569 1 0 Msboxa [(-0.125)] .17004 .18231 1 0 Msboxa [(-0.1)] .17004 .26894 1 0 Msboxa [(-0.075)] .17004 .35556 1 0 Msboxa [(-0.05)] .17004 .44219 1 0 Msboxa [(-0.025)] .17004 .52881 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .62269 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .61544 m .05026 .62169 L s P [(2.5)] .05026 .61544 0 2 Mshowa p .002 w .31481 .61544 m .31481 .62169 L s P [(7.5)] .31481 .61544 0 2 Mshowa p .002 w .44709 .61544 m .44709 .62169 L s P [(10)] .44709 .61544 0 2 Mshowa p .002 w .57937 .61544 m .57937 .62169 L s P [(12.5)] .57937 .61544 0 2 Mshowa p .002 w .71164 .61544 m .71164 .62169 L s P [(15)] .71164 .61544 0 2 Mshowa p .002 w .84392 .61544 m .84392 .62169 L s P [(17.5)] .84392 .61544 0 2 Mshowa p .002 w .97619 .61544 m .97619 .62169 L s P [(20)] .97619 .61544 0 2 Mshowa p .001 w .07672 .61544 m .07672 .61919 L s P p .001 w .10317 .61544 m .10317 .61919 L s P p .001 w .12963 .61544 m .12963 .61919 L s P p .001 w .15608 .61544 m .15608 .61919 L s P p .001 w .20899 .61544 m .20899 .61919 L s P p .001 w .23545 .61544 m .23545 .61919 L s P p .001 w .2619 .61544 m .2619 .61919 L s P p .001 w .28836 .61544 m .28836 .61919 L s P p .001 w .34127 .61544 m .34127 .61919 L s P p .001 w .36772 .61544 m .36772 .61919 L s P p .001 w .39418 .61544 m .39418 .61919 L s P p .001 w .42063 .61544 m .42063 .61919 L s P p .001 w .47354 .61544 m .47354 .61919 L s P p .001 w .5 .61544 m .5 .61919 L s P p .001 w .52646 .61544 m .52646 .61919 L s P p .001 w .55291 .61544 m .55291 .61919 L s P p .001 w .60582 .61544 m .60582 .61919 L s P p .001 w .63228 .61544 m .63228 .61919 L s P p .001 w .65873 .61544 m .65873 .61919 L s P p .001 w .68519 .61544 m .68519 .61919 L s P p .001 w .7381 .61544 m .7381 .61919 L s P p .001 w .76455 .61544 m .76455 .61919 L s P p .001 w .79101 .61544 m .79101 .61919 L s P p .001 w .81746 .61544 m .81746 .61919 L s P p .001 w .87037 .61544 m .87037 .61919 L s P p .001 w .89683 .61544 m .89683 .61919 L s P p .001 w .92328 .61544 m .92328 .61919 L s P p .001 w .94974 .61544 m .94974 .61919 L s P p .001 w .02381 .61544 m .02381 .61919 L s P p .002 w 0 .61544 m 1 .61544 L s P p .002 w .18254 .00906 m .18879 .00906 L s P [(-0.175)] .17004 .00906 1 0 Mshowa p .002 w .18254 .09569 m .18879 .09569 L s P [(-0.15)] .17004 .09569 1 0 Mshowa p .002 w .18254 .18231 m .18879 .18231 L s P [(-0.125)] .17004 .18231 1 0 Mshowa p .002 w .18254 .26894 m .18879 .26894 L s P [(-0.1)] .17004 .26894 1 0 Mshowa p .002 w .18254 .35556 m .18879 .35556 L s P [(-0.075)] .17004 .35556 1 0 Mshowa p .002 w .18254 .44219 m .18879 .44219 L s P [(-0.05)] .17004 .44219 1 0 Mshowa p .002 w .18254 .52881 m .18879 .52881 L s P [(-0.025)] .17004 .52881 1 0 Mshowa p .001 w .18254 .02639 m .18629 .02639 L s P p .001 w .18254 .04371 m .18629 .04371 L s P p .001 w .18254 .06104 m .18629 .06104 L s P p .001 w .18254 .07836 m .18629 .07836 L s P p .001 w .18254 .11301 m .18629 .11301 L s P p .001 w .18254 .13034 m .18629 .13034 L s P p .001 w .18254 .14766 m .18629 .14766 L s P p .001 w .18254 .16499 m .18629 .16499 L s P p .001 w .18254 .19964 m .18629 .19964 L s P p .001 w .18254 .21696 m .18629 .21696 L s P p .001 w .18254 .23429 m .18629 .23429 L s P p .001 w .18254 .25161 m .18629 .25161 L s P p .001 w .18254 .28626 m .18629 .28626 L s P p .001 w .18254 .30359 m .18629 .30359 L s P p .001 w .18254 .32091 m .18629 .32091 L s P p .001 w .18254 .33824 m .18629 .33824 L s P p .001 w .18254 .37289 m .18629 .37289 L s P p .001 w .18254 .39021 m .18629 .39021 L s P p .001 w .18254 .40754 m .18629 .40754 L s P p .001 w .18254 .42486 m .18629 .42486 L s P p .001 w .18254 .45951 m .18629 .45951 L s P p .001 w .18254 .47684 m .18629 .47684 L s P p .001 w .18254 .49416 m .18629 .49416 L s P p .001 w .18254 .51149 m .18629 .51149 L s P p .001 w .18254 .54614 m .18629 .54614 L s P p .001 w .18254 .56346 m .18629 .56346 L s P p .001 w .18254 .58079 m .18629 .58079 L s P p .001 w .18254 .59811 m .18629 .59811 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .11505 L .04365 .19187 L .05357 .25193 L .06349 .29977 L .07341 .33848 L .08333 .37023 L .10317 .41873 L .12302 .45354 L .14286 .47937 L .1627 .49909 L .18254 .5145 L .20238 .52677 L .22222 .53673 L .2619 .55174 L .30159 .56241 L .34127 .57029 L .38095 .5763 L .42063 .58101 L .46032 .58478 L .5 .58785 L .53968 .5904 L .57937 .59254 L .61905 .59437 L .65873 .59593 L .69841 .59729 L .7381 .59848 L .77778 .59953 L .81746 .60046 L .85714 .60129 L .89683 .60203 L .93651 .60271 L .97619 .60332 L s P P p 1 0 0 r p .02 w .87186 .60158 Mdot .6778 .59661 Mdot .64419 .59539 Mdot .3918 .5777 Mdot .64865 .59556 Mdot .36594 .57421 Mdot .47165 .58572 Mdot .48564 .58681 Mdot .85275 .6012 Mdot .78 .59958 Mdot .59721 .5934 Mdot .46956 .58555 Mdot .70154 .59739 Mdot .76796 .59928 Mdot .08116 .3638 Mdot .39997 .5787 Mdot .40871 .57971 Mdot .05784 .27383 Mdot .91324 .60232 Mdot .54839 .59091 Mdot .0843 .37302 Mdot .74458 .59866 Mdot .46328 .58503 Mdot .29883 .56178 Mdot .73792 .59848 Mdot .50706 .58835 Mdot .68057 .5967 Mdot .1309 .4647 Mdot .29567 .56103 Mdot .09687 .40511 Mdot .35862 .57311 Mdot .52467 .5895 Mdot .29813 .56161 Mdot .75628 .59898 Mdot .81104 .60032 Mdot .48126 .58648 Mdot .09693 .40524 Mdot .95039 .60293 Mdot .39255 .5778 Mdot .84329 .60101 Mdot .12068 .44999 Mdot .77023 .59934 Mdot .73826 .59849 Mdot .26937 .55403 Mdot .46394 .58509 Mdot .48928 .58709 Mdot .62909 .59479 Mdot .03183 .09795 Mdot .55597 .59133 Mdot .05764 .27282 Mdot .18723 .51766 Mdot .19053 .51978 Mdot .09036 .38938 Mdot .39649 .57828 Mdot .32829 .56796 Mdot .15093 .48803 Mdot .84054 .60095 Mdot .63784 .59514 Mdot .23527 .54229 Mdot .62982 .59482 Mdot .65549 .59581 Mdot .59599 .59334 Mdot .49207 .58729 Mdot .42927 .5819 Mdot .42552 .58152 Mdot .05658 .2676 Mdot .12104 .45054 Mdot .88184 .60176 Mdot .25669 .55007 Mdot .26781 .55357 Mdot .04523 .20242 Mdot .61418 .59416 Mdot .36593 .57421 Mdot .042 .18043 Mdot .85006 .60115 Mdot .58893 .59301 Mdot .44835 .58373 Mdot .55615 .59134 Mdot .35112 .57194 Mdot .6505 .59563 Mdot .11824 .44613 Mdot .4065 .57946 Mdot .09682 .405 Mdot .7685 .5993 Mdot .70708 .59757 Mdot .75031 .59882 Mdot .83321 .6008 Mdot .1852 .51631 Mdot .40868 .57971 Mdot .60523 .59377 Mdot .08137 .36442 Mdot .93092 .60262 Mdot .93932 .60276 Mdot .54824 .5909 Mdot .86631 .60147 Mdot .75592 .59897 Mdot .18304 .51485 Mdot .02942 .07487 Mdot .32602 .56753 Mdot .82041 .60052 Mdot .89353 .60198 Mdot .23899 .54375 Mdot .83597 .60086 Mdot .28426 .55817 Mdot .87284 .6016 Mdot .71221 .59772 Mdot .03035 .08392 Mdot .93248 .60264 Mdot .8235 .60059 Mdot .92687 .60255 Mdot .52129 .58928 Mdot .13026 .46384 Mdot .60395 .59371 Mdot .86746 .60149 Mdot .74417 .59865 Mdot .60702 .59385 Mdot .84752 .6011 Mdot .871 .60156 Mdot .84098 .60096 Mdot .91471 .60235 Mdot .04129 .17534 Mdot .96403 .60314 Mdot .49619 .58759 Mdot .85758 .6013 Mdot .86843 .60151 Mdot .9663 .60317 Mdot .14807 .48507 Mdot .3831 .57659 Mdot .27674 .55614 Mdot .48829 .58701 Mdot .41195 .58008 Mdot .54977 .59098 Mdot .39447 .57804 Mdot .56192 .59165 Mdot .87447 .60163 Mdot .68054 .5967 Mdot .02985 .07911 Mdot .69042 .59703 Mdot .85797 .60131 Mdot .33114 .56849 Mdot .60504 .59376 Mdot .81904 .60049 Mdot .2169 .53426 Mdot .29308 .56041 Mdot .79863 .60003 Mdot .70735 .59757 Mdot .90035 .6021 Mdot .05062 .23554 Mdot .89431 .60199 Mdot .33639 .56944 Mdot .06014 .28478 Mdot .02906 .07132 Mdot .43129 .58211 Mdot .18621 .51699 Mdot .23819 .54344 Mdot .86932 .60153 Mdot .41576 .58049 Mdot .18578 .5167 Mdot .4916 .58726 Mdot .15897 .49577 Mdot .57348 .59225 Mdot .79877 .60003 Mdot .53715 .59026 Mdot .79351 .59991 Mdot .12967 .46304 Mdot .63111 .59487 Mdot .86767 .6015 Mdot .73798 .59848 Mdot .47572 .58605 Mdot .57601 .59238 Mdot .96031 .60308 Mdot .44086 .58304 Mdot .41064 .57993 Mdot .61828 .59434 Mdot .84968 .60114 Mdot .80095 .60009 Mdot .74382 .59864 Mdot .19365 .52171 Mdot .85234 .60119 Mdot .43186 .58216 Mdot .40043 .57875 Mdot .83204 .60078 Mdot .41631 .58055 Mdot .41499 .58041 Mdot .02948 .07548 Mdot .7729 .59941 Mdot .15599 .49299 Mdot .94814 .6029 Mdot .38836 .57727 Mdot .7783 .59954 Mdot .51221 .58869 Mdot .37025 .57483 Mdot .13558 .47075 Mdot .5144 .58884 Mdot .69546 .5972 Mdot .12322 .45385 Mdot .68979 .59701 Mdot .32651 .56762 Mdot .55761 .59142 Mdot .35456 .57248 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] -6 -6 Errors:{-0.0000423732, 1.06084 10 , 0, 1.06084 10 } :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.057075206413038, A -> 1.000001060844471, R -> 5} ;[o] {a -> 1.05708, A -> 1., R -> 5} :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] angle fixed, fit for Hamaker and radius, the goal is A->1, R->5: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={a->1.05712}; sol=Sort[Join[NonlinearFit[dat3,ptmodel/.cond,d,{R,A}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp3,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-5)/5,(A R-5)/5}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.615435 3.46511 [ [(2.5)] .05026 .61544 0 2 Msboxa [(7.5)] .31481 .61544 0 2 Msboxa [(10)] .44709 .61544 0 2 Msboxa [(12.5)] .57937 .61544 0 2 Msboxa [(15)] .71164 .61544 0 2 Msboxa [(17.5)] .84392 .61544 0 2 Msboxa [(20)] .97619 .61544 0 2 Msboxa [(-0.175)] .17004 .00904 1 0 Msboxa [(-0.15)] .17004 .09567 1 0 Msboxa [(-0.125)] .17004 .1823 1 0 Msboxa [(-0.1)] .17004 .26892 1 0 Msboxa [(-0.075)] .17004 .35555 1 0 Msboxa [(-0.05)] .17004 .44218 1 0 Msboxa [(-0.025)] .17004 .52881 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .62269 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .61544 m .05026 .62169 L s P [(2.5)] .05026 .61544 0 2 Mshowa p .002 w .31481 .61544 m .31481 .62169 L s P [(7.5)] .31481 .61544 0 2 Mshowa p .002 w .44709 .61544 m .44709 .62169 L s P [(10)] .44709 .61544 0 2 Mshowa p .002 w .57937 .61544 m .57937 .62169 L s P [(12.5)] .57937 .61544 0 2 Mshowa p .002 w .71164 .61544 m .71164 .62169 L s P [(15)] .71164 .61544 0 2 Mshowa p .002 w .84392 .61544 m .84392 .62169 L s P [(17.5)] .84392 .61544 0 2 Mshowa p .002 w .97619 .61544 m .97619 .62169 L s P [(20)] .97619 .61544 0 2 Mshowa p .001 w .07672 .61544 m .07672 .61919 L s P p .001 w .10317 .61544 m .10317 .61919 L s P p .001 w .12963 .61544 m .12963 .61919 L s P p .001 w .15608 .61544 m .15608 .61919 L s P p .001 w .20899 .61544 m .20899 .61919 L s P p .001 w .23545 .61544 m .23545 .61919 L s P p .001 w .2619 .61544 m .2619 .61919 L s P p .001 w .28836 .61544 m .28836 .61919 L s P p .001 w .34127 .61544 m .34127 .61919 L s P p .001 w .36772 .61544 m .36772 .61919 L s P p .001 w .39418 .61544 m .39418 .61919 L s P p .001 w .42063 .61544 m .42063 .61919 L s P p .001 w .47354 .61544 m .47354 .61919 L s P p .001 w .5 .61544 m .5 .61919 L s P p .001 w .52646 .61544 m .52646 .61919 L s P p .001 w .55291 .61544 m .55291 .61919 L s P p .001 w .60582 .61544 m .60582 .61919 L s P p .001 w .63228 .61544 m .63228 .61919 L s P p .001 w .65873 .61544 m .65873 .61919 L s P p .001 w .68519 .61544 m .68519 .61919 L s P p .001 w .7381 .61544 m .7381 .61919 L s P p .001 w .76455 .61544 m .76455 .61919 L s P p .001 w .79101 .61544 m .79101 .61919 L s P p .001 w .81746 .61544 m .81746 .61919 L s P p .001 w .87037 .61544 m .87037 .61919 L s P p .001 w .89683 .61544 m .89683 .61919 L s P p .001 w .92328 .61544 m .92328 .61919 L s P p .001 w .94974 .61544 m .94974 .61919 L s P p .001 w .02381 .61544 m .02381 .61919 L s P p .002 w 0 .61544 m 1 .61544 L s P p .002 w .18254 .00904 m .18879 .00904 L s P [(-0.175)] .17004 .00904 1 0 Mshowa p .002 w .18254 .09567 m .18879 .09567 L s P [(-0.15)] .17004 .09567 1 0 Mshowa p .002 w .18254 .1823 m .18879 .1823 L s P [(-0.125)] .17004 .1823 1 0 Mshowa p .002 w .18254 .26892 m .18879 .26892 L s P [(-0.1)] .17004 .26892 1 0 Mshowa p .002 w .18254 .35555 m .18879 .35555 L s P [(-0.075)] .17004 .35555 1 0 Mshowa p .002 w .18254 .44218 m .18879 .44218 L s P [(-0.05)] .17004 .44218 1 0 Mshowa p .002 w .18254 .52881 m .18879 .52881 L s P [(-0.025)] .17004 .52881 1 0 Mshowa p .001 w .18254 .02637 m .18629 .02637 L s P p .001 w .18254 .04369 m .18629 .04369 L s P p .001 w .18254 .06102 m .18629 .06102 L s P p .001 w .18254 .07834 m .18629 .07834 L s P p .001 w .18254 .11299 m .18629 .11299 L s P p .001 w .18254 .13032 m .18629 .13032 L s P p .001 w .18254 .14765 m .18629 .14765 L s P p .001 w .18254 .16497 m .18629 .16497 L s P p .001 w .18254 .19962 m .18629 .19962 L s P p .001 w .18254 .21695 m .18629 .21695 L s P p .001 w .18254 .23427 m .18629 .23427 L s P p .001 w .18254 .2516 m .18629 .2516 L s P p .001 w .18254 .28625 m .18629 .28625 L s P p .001 w .18254 .30358 m .18629 .30358 L s P p .001 w .18254 .3209 m .18629 .3209 L s P p .001 w .18254 .33823 m .18629 .33823 L s P p .001 w .18254 .37288 m .18629 .37288 L s P p .001 w .18254 .3902 m .18629 .3902 L s P p .001 w .18254 .40753 m .18629 .40753 L s P p .001 w .18254 .42485 m .18629 .42485 L s P p .001 w .18254 .45951 m .18629 .45951 L s P p .001 w .18254 .47683 m .18629 .47683 L s P p .001 w .18254 .49416 m .18629 .49416 L s P p .001 w .18254 .51148 m .18629 .51148 L s P p .001 w .18254 .54613 m .18629 .54613 L s P p .001 w .18254 .56346 m .18629 .56346 L s P p .001 w .18254 .58078 m .18629 .58078 L s P p .001 w .18254 .59811 m .18629 .59811 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .11506 L .04365 .19187 L .05357 .25194 L .06349 .29977 L .07341 .33848 L .08333 .37023 L .10317 .41873 L .12302 .45354 L .14286 .47938 L .1627 .4991 L .18254 .5145 L .20238 .52678 L .22222 .53673 L .2619 .55175 L .30159 .56241 L .34127 .57029 L .38095 .5763 L .42063 .58101 L .46032 .58478 L .5 .58786 L .53968 .5904 L .57937 .59254 L .61905 .59437 L .65873 .59593 L .69841 .59729 L .7381 .59848 L .77778 .59953 L .81746 .60046 L .85714 .60129 L .89683 .60203 L .93651 .60271 L .97619 .60332 L s P P p 1 0 0 r p .02 w .87186 .60157 Mdot .6778 .59661 Mdot .64419 .59538 Mdot .3918 .5777 Mdot .64865 .59556 Mdot .36594 .57421 Mdot .47165 .58572 Mdot .48564 .58681 Mdot .85275 .6012 Mdot .78 .59958 Mdot .59721 .5934 Mdot .46956 .58555 Mdot .70154 .59739 Mdot .76796 .59928 Mdot .08116 .36379 Mdot .39997 .5787 Mdot .40871 .57971 Mdot .05784 .27382 Mdot .91324 .60232 Mdot .54839 .5909 Mdot .0843 .37301 Mdot .74458 .59866 Mdot .46328 .58503 Mdot .29883 .56178 Mdot .73792 .59847 Mdot .50706 .58834 Mdot .68057 .5967 Mdot .1309 .4647 Mdot .29567 .56103 Mdot .09687 .40511 Mdot .35862 .57311 Mdot .52467 .58949 Mdot .29813 .56161 Mdot .75628 .59898 Mdot .81104 .60031 Mdot .48126 .58648 Mdot .09693 .40523 Mdot .95039 .60293 Mdot .39255 .5778 Mdot .84329 .60101 Mdot .12068 .44999 Mdot .77023 .59934 Mdot .73826 .59848 Mdot .26937 .55403 Mdot .46394 .58509 Mdot .48928 .58708 Mdot .62909 .59478 Mdot .03183 .09794 Mdot .55597 .59133 Mdot .05764 .27281 Mdot .18723 .51765 Mdot .19053 .51977 Mdot .09036 .38938 Mdot .39649 .57828 Mdot .32829 .56796 Mdot .15093 .48802 Mdot .84054 .60095 Mdot .63784 .59514 Mdot .23527 .54229 Mdot .62982 .59481 Mdot .65549 .59581 Mdot .59599 .59334 Mdot .49207 .58729 Mdot .42927 .5819 Mdot .42552 .58152 Mdot .05658 .26759 Mdot .12104 .45053 Mdot .88184 .60176 Mdot .25669 .55007 Mdot .26781 .55356 Mdot .04523 .20241 Mdot .61418 .59416 Mdot .36593 .57421 Mdot .042 .18041 Mdot .85006 .60115 Mdot .58893 .59301 Mdot .44835 .58373 Mdot .55615 .59134 Mdot .35112 .57193 Mdot .6505 .59563 Mdot .11824 .44612 Mdot .4065 .57946 Mdot .09682 .40499 Mdot .7685 .59929 Mdot .70708 .59756 Mdot .75031 .59882 Mdot .83321 .6008 Mdot .1852 .51631 Mdot .40868 .57971 Mdot .60523 .59376 Mdot .08137 .36441 Mdot .93092 .60262 Mdot .93932 .60275 Mdot .54824 .5909 Mdot .86631 .60147 Mdot .75592 .59897 Mdot .18304 .51484 Mdot .02942 .07486 Mdot .32602 .56753 Mdot .82041 .60052 Mdot .89353 .60198 Mdot .23899 .54375 Mdot .83597 .60086 Mdot .28426 .55817 Mdot .87284 .60159 Mdot .71221 .59772 Mdot .03035 .0839 Mdot .93248 .60264 Mdot .8235 .60059 Mdot .92687 .60255 Mdot .52129 .58928 Mdot .13026 .46384 Mdot .60395 .59371 Mdot .86746 .60149 Mdot .74417 .59865 Mdot .60702 .59384 Mdot .84752 .6011 Mdot .871 .60156 Mdot .84098 .60096 Mdot .91471 .60235 Mdot .04129 .17532 Mdot .96403 .60314 Mdot .49619 .58759 Mdot .85758 .6013 Mdot .86843 .60151 Mdot .9663 .60317 Mdot .14807 .48506 Mdot .3831 .57659 Mdot .27674 .55614 Mdot .48829 .58701 Mdot .41195 .58007 Mdot .54977 .59098 Mdot .39447 .57803 Mdot .56192 .59165 Mdot .87447 .60162 Mdot .68054 .5967 Mdot .02985 .07909 Mdot .69042 .59703 Mdot .85797 .6013 Mdot .33114 .56849 Mdot .60504 .59376 Mdot .81904 .60049 Mdot .2169 .53426 Mdot .29308 .5604 Mdot .79863 .60003 Mdot .70735 .59757 Mdot .90035 .6021 Mdot .05062 .23553 Mdot .89431 .60199 Mdot .33639 .56944 Mdot .06014 .28477 Mdot .02906 .07131 Mdot .43129 .5821 Mdot .18621 .51698 Mdot .23819 .54344 Mdot .86932 .60153 Mdot .41576 .58049 Mdot .18578 .51669 Mdot .4916 .58725 Mdot .15897 .49576 Mdot .57348 .59225 Mdot .79877 .60003 Mdot .53715 .59025 Mdot .79351 .59991 Mdot .12967 .46304 Mdot .63111 .59487 Mdot .86767 .60149 Mdot .73798 .59848 Mdot .47572 .58605 Mdot .57601 .59238 Mdot .96031 .60308 Mdot .44086 .58303 Mdot .41064 .57993 Mdot .61828 .59433 Mdot .84968 .60114 Mdot .80095 .60008 Mdot .74382 .59864 Mdot .19365 .5217 Mdot .85234 .60119 Mdot .43186 .58216 Mdot .40043 .57875 Mdot .83204 .60077 Mdot .41631 .58055 Mdot .41499 .58041 Mdot .02948 .07546 Mdot .7729 .59941 Mdot .15599 .49299 Mdot .94814 .60289 Mdot .38836 .57727 Mdot .7783 .59954 Mdot .51221 .58869 Mdot .37025 .57483 Mdot .13558 .47074 Mdot .5144 .58883 Mdot .69546 .5972 Mdot .12322 .45384 Mdot .68979 .59701 Mdot .32651 .56762 Mdot .55761 .59142 Mdot .35456 .57248 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] -6 -6 -6 Errors:{0., 4.74334 10 , -8.63491 10 , -3.89161 10 } :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.05712, A -> 1.000004743341057, R -> 4.999956825452255} ;[o] {a -> 1.05712, A -> 1., R -> 4.99996} :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] the results are not very good, but the product R A is good: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] R A/.% :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] 4.999980541952747 ;[o] 4.99998 :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] fixing hamaker and fitting for the radius and angle is almost equivalent to fixing the radius, so it is not worth doing again. :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] now everything free: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={}; sol=Sort[Join[NonlinearFit[dat3,ptmodel/.cond,d,{a,R,A}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp3,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-5)/5,(A R-5)/5}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.614861 3.31197 [ [(2.5)] .05026 .61486 0 2 Msboxa [(7.5)] .31481 .61486 0 2 Msboxa [(10)] .44709 .61486 0 2 Msboxa [(12.5)] .57937 .61486 0 2 Msboxa [(15)] .71164 .61486 0 2 Msboxa [(17.5)] .84392 .61486 0 2 Msboxa [(20)] .97619 .61486 0 2 Msboxa [(-0.175)] .17004 .03527 1 0 Msboxa [(-0.15)] .17004 .11807 1 0 Msboxa [(-0.125)] .17004 .20086 1 0 Msboxa [(-0.1)] .17004 .28366 1 0 Msboxa [(-0.075)] .17004 .36646 1 0 Msboxa [(-0.05)] .17004 .44926 1 0 Msboxa [(-0.025)] .17004 .53206 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .62211 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .61486 m .05026 .62111 L s P [(2.5)] .05026 .61486 0 2 Mshowa p .002 w .31481 .61486 m .31481 .62111 L s P [(7.5)] .31481 .61486 0 2 Mshowa p .002 w .44709 .61486 m .44709 .62111 L s P [(10)] .44709 .61486 0 2 Mshowa p .002 w .57937 .61486 m .57937 .62111 L s P [(12.5)] .57937 .61486 0 2 Mshowa p .002 w .71164 .61486 m .71164 .62111 L s P [(15)] .71164 .61486 0 2 Mshowa p .002 w .84392 .61486 m .84392 .62111 L s P [(17.5)] .84392 .61486 0 2 Mshowa p .002 w .97619 .61486 m .97619 .62111 L s P [(20)] .97619 .61486 0 2 Mshowa p .001 w .07672 .61486 m .07672 .61861 L s P p .001 w .10317 .61486 m .10317 .61861 L s P p .001 w .12963 .61486 m .12963 .61861 L s P p .001 w .15608 .61486 m .15608 .61861 L s P p .001 w .20899 .61486 m .20899 .61861 L s P p .001 w .23545 .61486 m .23545 .61861 L s P p .001 w .2619 .61486 m .2619 .61861 L s P p .001 w .28836 .61486 m .28836 .61861 L s P p .001 w .34127 .61486 m .34127 .61861 L s P p .001 w .36772 .61486 m .36772 .61861 L s P p .001 w .39418 .61486 m .39418 .61861 L s P p .001 w .42063 .61486 m .42063 .61861 L s P p .001 w .47354 .61486 m .47354 .61861 L s P p .001 w .5 .61486 m .5 .61861 L s P p .001 w .52646 .61486 m .52646 .61861 L s P p .001 w .55291 .61486 m .55291 .61861 L s P p .001 w .60582 .61486 m .60582 .61861 L s P p .001 w .63228 .61486 m .63228 .61861 L s P p .001 w .65873 .61486 m .65873 .61861 L s P p .001 w .68519 .61486 m .68519 .61861 L s P p .001 w .7381 .61486 m .7381 .61861 L s P p .001 w .76455 .61486 m .76455 .61861 L s P p .001 w .79101 .61486 m .79101 .61861 L s P p .001 w .81746 .61486 m .81746 .61861 L s P p .001 w .87037 .61486 m .87037 .61861 L s P p .001 w .89683 .61486 m .89683 .61861 L s P p .001 w .92328 .61486 m .92328 .61861 L s P p .001 w .94974 .61486 m .94974 .61861 L s P p .001 w .02381 .61486 m .02381 .61861 L s P p .002 w 0 .61486 m 1 .61486 L s P p .002 w .18254 .03527 m .18879 .03527 L s P [(-0.175)] .17004 .03527 1 0 Mshowa p .002 w .18254 .11807 m .18879 .11807 L s P [(-0.15)] .17004 .11807 1 0 Mshowa p .002 w .18254 .20086 m .18879 .20086 L s P [(-0.125)] .17004 .20086 1 0 Mshowa p .002 w .18254 .28366 m .18879 .28366 L s P [(-0.1)] .17004 .28366 1 0 Mshowa p .002 w .18254 .36646 m .18879 .36646 L s P [(-0.075)] .17004 .36646 1 0 Mshowa p .002 w .18254 .44926 m .18879 .44926 L s P [(-0.05)] .17004 .44926 1 0 Mshowa p .002 w .18254 .53206 m .18879 .53206 L s P [(-0.025)] .17004 .53206 1 0 Mshowa p .001 w .18254 .05183 m .18629 .05183 L s P p .001 w .18254 .06839 m .18629 .06839 L s P p .001 w .18254 .08495 m .18629 .08495 L s P p .001 w .18254 .10151 m .18629 .10151 L s P p .001 w .18254 .13463 m .18629 .13463 L s P p .001 w .18254 .15119 m .18629 .15119 L s P p .001 w .18254 .16775 m .18629 .16775 L s P p .001 w .18254 .18431 m .18629 .18431 L s P p .001 w .18254 .21742 m .18629 .21742 L s P p .001 w .18254 .23398 m .18629 .23398 L s P p .001 w .18254 .25054 m .18629 .25054 L s P p .001 w .18254 .2671 m .18629 .2671 L s P p .001 w .18254 .30022 m .18629 .30022 L s P p .001 w .18254 .31678 m .18629 .31678 L s P p .001 w .18254 .33334 m .18629 .33334 L s P p .001 w .18254 .3499 m .18629 .3499 L s P p .001 w .18254 .38302 m .18629 .38302 L s P p .001 w .18254 .39958 m .18629 .39958 L s P p .001 w .18254 .41614 m .18629 .41614 L s P p .001 w .18254 .4327 m .18629 .4327 L s P p .001 w .18254 .46582 m .18629 .46582 L s P p .001 w .18254 .48238 m .18629 .48238 L s P p .001 w .18254 .49894 m .18629 .49894 L s P p .001 w .18254 .5155 m .18629 .5155 L s P p .001 w .18254 .54862 m .18629 .54862 L s P p .001 w .18254 .56518 m .18629 .56518 L s P p .001 w .18254 .58174 m .18629 .58174 L s P p .001 w .18254 .5983 m .18629 .5983 L s P p .001 w .18254 .01871 m .18629 .01871 L s P p .001 w .18254 .00215 m .18629 .00215 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .11532 L .04365 .19234 L .05357 .25256 L .06349 .30051 L .07341 .33931 L .08333 .37112 L .10317 .4197 L .12302 .45454 L .14286 .48038 L .1627 .50008 L .18254 .51546 L .20238 .5277 L .22222 .53762 L .2619 .55255 L .30159 .56313 L .34127 .57094 L .38095 .57688 L .42063 .58152 L .46032 .58523 L .5 .58825 L .53968 .59075 L .57937 .59284 L .61905 .59462 L .65873 .59615 L .69841 .59747 L .7381 .59863 L .77778 .59965 L .81746 .60055 L .85714 .60135 L .89683 .60208 L .93651 .60273 L .97619 .60332 L s P P p 1 0 0 r p .02 w .87186 .60161 Mdot .6778 .59687 Mdot .64419 .5957 Mdot .3918 .5788 Mdot .64865 .59586 Mdot .36594 .57546 Mdot .47165 .58646 Mdot .48564 .5875 Mdot .85275 .60126 Mdot .78 .59971 Mdot .59721 .5938 Mdot .46956 .5863 Mdot .70154 .59761 Mdot .76796 .59942 Mdot .08116 .37434 Mdot .39997 .57975 Mdot .40871 .58071 Mdot .05784 .28835 Mdot .91324 .60233 Mdot .54839 .59141 Mdot .0843 .38315 Mdot .74458 .59883 Mdot .46328 .5858 Mdot .29883 .56357 Mdot .73792 .59865 Mdot .50706 .58897 Mdot .68057 .59696 Mdot .1309 .47078 Mdot .29567 .56286 Mdot .09687 .41383 Mdot .35862 .57441 Mdot .52467 .59007 Mdot .29813 .56342 Mdot .75628 .59913 Mdot .81104 .60041 Mdot .48126 .58718 Mdot .09693 .41395 Mdot .95039 .60291 Mdot .39255 .57888 Mdot .84329 .60107 Mdot .12068 .45672 Mdot .77023 .59947 Mdot .73826 .59866 Mdot .26937 .55617 Mdot .46394 .58585 Mdot .48928 .58776 Mdot .62909 .59512 Mdot .03183 .12023 Mdot .55597 .59182 Mdot .05764 .28738 Mdot .18723 .5214 Mdot .19053 .52343 Mdot .09036 .39879 Mdot .39649 .57935 Mdot .32829 .56948 Mdot .15093 .49308 Mdot .84054 .60102 Mdot .63784 .59546 Mdot .23527 .54494 Mdot .62982 .59515 Mdot .65549 .5961 Mdot .59599 .59374 Mdot .49207 .58796 Mdot .42927 .58281 Mdot .42552 .58245 Mdot .05658 .28239 Mdot .12104 .45725 Mdot .88184 .60179 Mdot .25669 .55238 Mdot .26781 .55572 Mdot .04523 .22009 Mdot .61418 .59452 Mdot .36593 .57546 Mdot .042 .19906 Mdot .85006 .6012 Mdot .58893 .59343 Mdot .44835 .58455 Mdot .55615 .59183 Mdot .35112 .57328 Mdot .6505 .59593 Mdot .11824 .45303 Mdot .4065 .58047 Mdot .09682 .41372 Mdot .7685 .59943 Mdot .70708 .59778 Mdot .75031 .59898 Mdot .83321 .60087 Mdot .1852 .52011 Mdot .40868 .58071 Mdot .60523 .59415 Mdot .08137 .37493 Mdot .93092 .60261 Mdot .93932 .60274 Mdot .54824 .59141 Mdot .86631 .60151 Mdot .75592 .59912 Mdot .18304 .51872 Mdot .02942 .09817 Mdot .32602 .56907 Mdot .82041 .60061 Mdot .89353 .602 Mdot .23899 .54634 Mdot .83597 .60093 Mdot .28426 .56013 Mdot .87284 .60163 Mdot .71221 .59793 Mdot .03035 .10682 Mdot .93248 .60263 Mdot .8235 .60067 Mdot .92687 .60255 Mdot .52129 .58986 Mdot .13026 .46996 Mdot .60395 .59409 Mdot .86746 .60153 Mdot .74417 .59882 Mdot .60702 .59422 Mdot .84752 .60115 Mdot .871 .6016 Mdot .84098 .60103 Mdot .91471 .60235 Mdot .04129 .1942 Mdot .96403 .60311 Mdot .49619 .58824 Mdot .85758 .60135 Mdot .86843 .60155 Mdot .9663 .60314 Mdot .14807 .49025 Mdot .3831 .57773 Mdot .27674 .55819 Mdot .48829 .58769 Mdot .41195 .58106 Mdot .54977 .59149 Mdot .39447 .57911 Mdot .56192 .59212 Mdot .87447 .60166 Mdot .68054 .59695 Mdot .02985 .10222 Mdot .69042 .59727 Mdot .85797 .60136 Mdot .33114 .56999 Mdot .60504 .59414 Mdot .81904 .60058 Mdot .2169 .53727 Mdot .29308 .56226 Mdot .79863 .60014 Mdot .70735 .59779 Mdot .90035 .60211 Mdot .05062 .25174 Mdot .89431 .60201 Mdot .33639 .5709 Mdot .06014 .29881 Mdot .02906 .09478 Mdot .43129 .583 Mdot .18621 .52076 Mdot .23819 .54605 Mdot .86932 .60157 Mdot .41576 .58146 Mdot .18578 .52048 Mdot .4916 .58792 Mdot .15897 .50048 Mdot .57348 .5927 Mdot .79877 .60014 Mdot .53715 .59079 Mdot .79351 .60002 Mdot .12967 .4692 Mdot .63111 .5952 Mdot .86767 .60154 Mdot .73798 .59865 Mdot .47572 .58677 Mdot .57601 .59282 Mdot .96031 .60305 Mdot .44086 .58389 Mdot .41064 .58092 Mdot .61828 .59469 Mdot .84968 .6012 Mdot .80095 .60019 Mdot .74382 .59881 Mdot .19365 .52527 Mdot .85234 .60125 Mdot .43186 .58306 Mdot .40043 .5798 Mdot .83204 .60085 Mdot .41631 .58152 Mdot .41499 .58138 Mdot .02948 .09875 Mdot .7729 .59954 Mdot .15599 .49782 Mdot .94814 .60287 Mdot .38836 .57838 Mdot .7783 .59967 Mdot .51221 .5893 Mdot .37025 .57605 Mdot .13558 .47656 Mdot .5144 .58944 Mdot .69546 .59743 Mdot .12322 .46041 Mdot .68979 .59725 Mdot .32651 .56916 Mdot .55761 .5919 Mdot .35456 .5738 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{-0.000533105, -0.0239861, 0.0671604, 0.0415634} :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.056556444391034, A -> 0.976013911902572, R -> 5.335802019749945} ;[o] {a -> 1.05656, A -> 0.976014, R -> 5.3358} :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] and the product A R is very good: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] R A/.% :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] 5.207817002433787 ;[o] 5.20782 :[font = subsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 12; ] dat4, data with 10% noise, small radius : :[font = input; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] dat4=Table[Module[{temp}, temp=Random[Real,{2,20}]; {temp,noise*0.011+ForceParametricTip[1.05712,temp,1,5]}], {x,1,200}]; :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] lp4=ListPlot[dat4,PlotRange->All,PlotStyle->{Red,PointSize[0.02]}] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0905923 0.054118 0.580479 3.41696 [ [(2.5)] .0447 .58048 0 2 Msboxa [(7.5)] .31529 .58048 0 2 Msboxa [(10)] .45059 .58048 0 2 Msboxa [(12.5)] .58588 .58048 0 2 Msboxa [(15)] .72118 .58048 0 2 Msboxa [(17.5)] .85647 .58048 0 2 Msboxa [(20)] .99177 .58048 0 2 Msboxa [(-0.15)] .1675 .06793 1 0 Msboxa [(-0.125)] .1675 .15336 1 0 Msboxa [(-0.1)] .1675 .23878 1 0 Msboxa [(-0.075)] .1675 .32421 1 0 Msboxa [(-0.05)] .1675 .40963 1 0 Msboxa [(-0.025)] .1675 .49505 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .0447 .58048 m .0447 .58673 L s P [(2.5)] .0447 .58048 0 2 Mshowa p .002 w .31529 .58048 m .31529 .58673 L s P [(7.5)] .31529 .58048 0 2 Mshowa p .002 w .45059 .58048 m .45059 .58673 L s P [(10)] .45059 .58048 0 2 Mshowa p .002 w .58588 .58048 m .58588 .58673 L s P [(12.5)] .58588 .58048 0 2 Mshowa p .002 w .72118 .58048 m .72118 .58673 L s P [(15)] .72118 .58048 0 2 Mshowa p .002 w .85647 .58048 m .85647 .58673 L s P [(17.5)] .85647 .58048 0 2 Mshowa p .002 w .99177 .58048 m .99177 .58673 L s P [(20)] .99177 .58048 0 2 Mshowa p .001 w .07176 .58048 m .07176 .58423 L s P p .001 w .09882 .58048 m .09882 .58423 L s P p .001 w .12588 .58048 m .12588 .58423 L s P p .001 w .15294 .58048 m .15294 .58423 L s P p .001 w .20706 .58048 m .20706 .58423 L s P p .001 w .23412 .58048 m .23412 .58423 L s P p .001 w .26117 .58048 m .26117 .58423 L s P p .001 w .28823 .58048 m .28823 .58423 L s P p .001 w .34235 .58048 m .34235 .58423 L s P p .001 w .36941 .58048 m .36941 .58423 L s P p .001 w .39647 .58048 m .39647 .58423 L s P p .001 w .42353 .58048 m .42353 .58423 L s P p .001 w .47765 .58048 m .47765 .58423 L s P p .001 w .50471 .58048 m .50471 .58423 L s P p .001 w .53176 .58048 m .53176 .58423 L s P p .001 w .55882 .58048 m .55882 .58423 L s P p .001 w .61294 .58048 m .61294 .58423 L s P p .001 w .64 .58048 m .64 .58423 L s P p .001 w .66706 .58048 m .66706 .58423 L s P p .001 w .69412 .58048 m .69412 .58423 L s P p .001 w .74824 .58048 m .74824 .58423 L s P p .001 w .7753 .58048 m .7753 .58423 L s P p .001 w .80235 .58048 m .80235 .58423 L s P p .001 w .82941 .58048 m .82941 .58423 L s P p .001 w .88353 .58048 m .88353 .58423 L s P p .001 w .91059 .58048 m .91059 .58423 L s P p .001 w .93765 .58048 m .93765 .58423 L s P p .001 w .96471 .58048 m .96471 .58423 L s P p .001 w .01764 .58048 m .01764 .58423 L s P p .002 w 0 .58048 m 1 .58048 L s P p .002 w .18 .06793 m .18625 .06793 L s P [(-0.15)] .1675 .06793 1 0 Mshowa p .002 w .18 .15336 m .18625 .15336 L s P [(-0.125)] .1675 .15336 1 0 Mshowa p .002 w .18 .23878 m .18625 .23878 L s P [(-0.1)] .1675 .23878 1 0 Mshowa p .002 w .18 .32421 m .18625 .32421 L s P [(-0.075)] .1675 .32421 1 0 Mshowa p .002 w .18 .40963 m .18625 .40963 L s P [(-0.05)] .1675 .40963 1 0 Mshowa p .002 w .18 .49505 m .18625 .49505 L s P [(-0.025)] .1675 .49505 1 0 Mshowa p .001 w .18 .08502 m .18375 .08502 L s P p .001 w .18 .1021 m .18375 .1021 L s P p .001 w .18 .11919 m .18375 .11919 L s P p .001 w .18 .13627 m .18375 .13627 L s P p .001 w .18 .17044 m .18375 .17044 L s P p .001 w .18 .18753 m .18375 .18753 L s P p .001 w .18 .20461 m .18375 .20461 L s P p .001 w .18 .2217 m .18375 .2217 L s P p .001 w .18 .25587 m .18375 .25587 L s P p .001 w .18 .27295 m .18375 .27295 L s P p .001 w .18 .29004 m .18375 .29004 L s P p .001 w .18 .30712 m .18375 .30712 L s P p .001 w .18 .34129 m .18375 .34129 L s P p .001 w .18 .35838 m .18375 .35838 L s P p .001 w .18 .37546 m .18375 .37546 L s P p .001 w .18 .39255 m .18375 .39255 L s P p .001 w .18 .42672 m .18375 .42672 L s P p .001 w .18 .4438 m .18375 .4438 L s P p .001 w .18 .46089 m .18375 .46089 L s P p .001 w .18 .47797 m .18375 .47797 L s P p .001 w .18 .51214 m .18375 .51214 L s P p .001 w .18 .52922 m .18375 .52922 L s P p .001 w .18 .54631 m .18375 .54631 L s P p .001 w .18 .56339 m .18375 .56339 L s P p .001 w .18 .05085 m .18375 .05085 L s P p .001 w .18 .03376 m .18375 .03376 L s P p .001 w .18 .01668 m .18375 .01668 L s P p .001 w .18 .59756 m .18375 .59756 L s P p .001 w .18 .61465 m .18375 .61465 L s P p .002 w .18 0 m .18 .61803 L s P P p 1 0 0 r p .02 w .67386 .57158 Mdot .90492 .53957 Mdot .26931 .49074 Mdot .66432 .54509 Mdot .15767 .45803 Mdot .88863 .54128 Mdot .47565 .53517 Mdot .36043 .53932 Mdot .89619 .60332 Mdot .03957 .14576 Mdot .48802 .54944 Mdot .12147 .39939 Mdot .80594 .59002 Mdot .20831 .4624 Mdot .55427 .57286 Mdot .53576 .57412 Mdot .41425 .53054 Mdot .6389 .56708 Mdot .93036 .53098 Mdot .29611 .49568 Mdot .05181 .20523 Mdot .27418 .55534 Mdot .55556 .56128 Mdot .17036 .49529 Mdot .74139 .58125 Mdot .95381 .58627 Mdot .20476 .45892 Mdot .43569 .58475 Mdot .78228 .58078 Mdot .78857 .55734 Mdot .84369 .58118 Mdot .5101 .53985 Mdot .80952 .57819 Mdot .25356 .50516 Mdot .27161 .49113 Mdot .10085 .3968 Mdot .52199 .54473 Mdot .13881 .47129 Mdot .33487 .5239 Mdot .69488 .58839 Mdot .54436 .56639 Mdot .89808 .56311 Mdot .69244 .59517 Mdot .40563 .5169 Mdot .87469 .54917 Mdot .16971 .49368 Mdot .62073 .53402 Mdot .0865 .33555 Mdot .11638 .43862 Mdot .93946 .55574 Mdot .77328 .55943 Mdot .26222 .51737 Mdot .24656 .5349 Mdot .35591 .50384 Mdot .54588 .5823 Mdot .94205 .59882 Mdot .66296 .56743 Mdot .78998 .53828 Mdot .05988 .27012 Mdot .72112 .59799 Mdot .93526 .57764 Mdot .77343 .57383 Mdot .17963 .45806 Mdot .52885 .52728 Mdot .92484 .5605 Mdot .19059 .45673 Mdot .3966 .54768 Mdot .24031 .48643 Mdot .72541 .59848 Mdot .4421 .51298 Mdot .68317 .56538 Mdot .71274 .52998 Mdot .73968 .59473 Mdot .93108 .56302 Mdot .5777 .57376 Mdot .41987 .53154 Mdot .64463 .54609 Mdot .24692 .49173 Mdot .26567 .53817 Mdot .02426 .02102 Mdot .53203 .5764 Mdot .57393 .52047 Mdot .84063 .54615 Mdot .85296 .5997 Mdot .11859 .40666 Mdot .91365 .56782 Mdot .53267 .56298 Mdot .51142 .53127 Mdot .87981 .55138 Mdot .28215 .55999 Mdot .63178 .56302 Mdot .27553 .55529 Mdot .11739 .43762 Mdot .69338 .5949 Mdot .26853 .52137 Mdot .83219 .60261 Mdot .16758 .44943 Mdot .91031 .56602 Mdot .62668 .56673 Mdot .41653 .52944 Mdot .73864 .54826 Mdot .15203 .44373 Mdot .12451 .44757 Mdot .86826 .55534 Mdot .02476 .09785 Mdot .19253 .4824 Mdot .748 .56216 Mdot .35211 .53194 Mdot .59807 .54194 Mdot .43357 .58077 Mdot .96316 .58215 Mdot .03468 .15258 Mdot .24216 .50309 Mdot .87443 .55629 Mdot .1353 .45174 Mdot .02381 .01472 Mdot .12818 .44321 Mdot .82305 .57274 Mdot .37195 .51454 Mdot .48859 .52732 Mdot .76564 .59244 Mdot .07266 .30085 Mdot .79425 .54152 Mdot .05562 .26313 Mdot .56973 .57961 Mdot .17297 .45391 Mdot .45207 .52412 Mdot .1668 .44847 Mdot .34154 .53821 Mdot .33552 .54221 Mdot .96136 .55878 Mdot .8387 .56055 Mdot .21336 .49965 Mdot .78368 .5422 Mdot .41088 .53107 Mdot .7457 .57372 Mdot .83291 .56673 Mdot .59038 .55273 Mdot .39848 .53227 Mdot .44122 .52829 Mdot .07459 .35092 Mdot .12335 .43698 Mdot .105 .39353 Mdot .27642 .50001 Mdot .88341 .59971 Mdot .4845 .51471 Mdot .49018 .56174 Mdot .73057 .55369 Mdot .64903 .53152 Mdot .15783 .49756 Mdot .26819 .53979 Mdot .70838 .57838 Mdot .21124 .53255 Mdot .60267 .52666 Mdot .12388 .42847 Mdot .3439 .5638 Mdot .23224 .51333 Mdot .85116 .59329 Mdot .73383 .59712 Mdot .13823 .44415 Mdot .10244 .41869 Mdot .97217 .57156 Mdot .82602 .53324 Mdot .28143 .54847 Mdot .63242 .53117 Mdot .67053 .58093 Mdot .52618 .5625 Mdot .34428 .56585 Mdot .31159 .49489 Mdot .48489 .51673 Mdot .56953 .52553 Mdot .3643 .54267 Mdot .48473 .53265 Mdot .3839 .57959 Mdot .65048 .53624 Mdot .12012 .39713 Mdot .03571 .11454 Mdot .44135 .54115 Mdot .50129 .54959 Mdot .11472 .41517 Mdot .20735 .49104 Mdot .6216 .55797 Mdot .62959 .54815 Mdot .27494 .55437 Mdot .1625 .43803 Mdot .8828 .59952 Mdot .50378 .58516 Mdot .78033 .58489 Mdot .48572 .57951 Mdot .35662 .52805 Mdot .97619 .56176 Mdot .25955 .54369 Mdot .78648 .60143 Mdot .62972 .5503 Mdot .17454 .48799 Mdot .37243 .5358 Mdot .02968 .10595 Mdot .48139 .55175 Mdot .51767 .56302 Mdot .69283 .57954 Mdot P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] No Input Form was saved for this expression. ;[o] -Graphics- :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] First test, radius fixed, fit for Hamaker and angle, the goal is A->1, a->1.05712: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={R->5}; sol=Sort[Join[NonlinearFit[dat4,ptmodel/.cond,d,{{a,0.5,Pi/2},A}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp4,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-5)/5,(A R-5)/5}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.581116 3.27406 [ [(2.5)] .05026 .58112 0 2 Msboxa [(7.5)] .31481 .58112 0 2 Msboxa [(10)] .44709 .58112 0 2 Msboxa [(12.5)] .57937 .58112 0 2 Msboxa [(15)] .71164 .58112 0 2 Msboxa [(17.5)] .84392 .58112 0 2 Msboxa [(20)] .97619 .58112 0 2 Msboxa [(-0.175)] .17004 .00816 1 0 Msboxa [(-0.15)] .17004 .09001 1 0 Msboxa [(-0.125)] .17004 .17186 1 0 Msboxa [(-0.1)] .17004 .25371 1 0 Msboxa [(-0.075)] .17004 .33556 1 0 Msboxa [(-0.05)] .17004 .41741 1 0 Msboxa [(-0.025)] .17004 .49926 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .58112 m .05026 .58737 L s P [(2.5)] .05026 .58112 0 2 Mshowa p .002 w .31481 .58112 m .31481 .58737 L s P [(7.5)] .31481 .58112 0 2 Mshowa p .002 w .44709 .58112 m .44709 .58737 L s P [(10)] .44709 .58112 0 2 Mshowa p .002 w .57937 .58112 m .57937 .58737 L s P [(12.5)] .57937 .58112 0 2 Mshowa p .002 w .71164 .58112 m .71164 .58737 L s P [(15)] .71164 .58112 0 2 Mshowa p .002 w .84392 .58112 m .84392 .58737 L s P [(17.5)] .84392 .58112 0 2 Mshowa p .002 w .97619 .58112 m .97619 .58737 L s P [(20)] .97619 .58112 0 2 Mshowa p .001 w .07672 .58112 m .07672 .58487 L s P p .001 w .10317 .58112 m .10317 .58487 L s P p .001 w .12963 .58112 m .12963 .58487 L s P p .001 w .15608 .58112 m .15608 .58487 L s P p .001 w .20899 .58112 m .20899 .58487 L s P p .001 w .23545 .58112 m .23545 .58487 L s P p .001 w .2619 .58112 m .2619 .58487 L s P p .001 w .28836 .58112 m .28836 .58487 L s P p .001 w .34127 .58112 m .34127 .58487 L s P p .001 w .36772 .58112 m .36772 .58487 L s P p .001 w .39418 .58112 m .39418 .58487 L s P p .001 w .42063 .58112 m .42063 .58487 L s P p .001 w .47354 .58112 m .47354 .58487 L s P p .001 w .5 .58112 m .5 .58487 L s P p .001 w .52646 .58112 m .52646 .58487 L s P p .001 w .55291 .58112 m .55291 .58487 L s P p .001 w .60582 .58112 m .60582 .58487 L s P p .001 w .63228 .58112 m .63228 .58487 L s P p .001 w .65873 .58112 m .65873 .58487 L s P p .001 w .68519 .58112 m .68519 .58487 L s P p .001 w .7381 .58112 m .7381 .58487 L s P p .001 w .76455 .58112 m .76455 .58487 L s P p .001 w .79101 .58112 m .79101 .58487 L s P p .001 w .81746 .58112 m .81746 .58487 L s P p .001 w .87037 .58112 m .87037 .58487 L s P p .001 w .89683 .58112 m .89683 .58487 L s P p .001 w .92328 .58112 m .92328 .58487 L s P p .001 w .94974 .58112 m .94974 .58487 L s P p .001 w .02381 .58112 m .02381 .58487 L s P p .002 w 0 .58112 m 1 .58112 L s P p .002 w .18254 .00816 m .18879 .00816 L s P [(-0.175)] .17004 .00816 1 0 Mshowa p .002 w .18254 .09001 m .18879 .09001 L s P [(-0.15)] .17004 .09001 1 0 Mshowa p .002 w .18254 .17186 m .18879 .17186 L s P [(-0.125)] .17004 .17186 1 0 Mshowa p .002 w .18254 .25371 m .18879 .25371 L s P [(-0.1)] .17004 .25371 1 0 Mshowa p .002 w .18254 .33556 m .18879 .33556 L s P [(-0.075)] .17004 .33556 1 0 Mshowa p .002 w .18254 .41741 m .18879 .41741 L s P [(-0.05)] .17004 .41741 1 0 Mshowa p .002 w .18254 .49926 m .18879 .49926 L s P [(-0.025)] .17004 .49926 1 0 Mshowa p .001 w .18254 .02453 m .18629 .02453 L s P p .001 w .18254 .0409 m .18629 .0409 L s P p .001 w .18254 .05727 m .18629 .05727 L s P p .001 w .18254 .07364 m .18629 .07364 L s P p .001 w .18254 .10638 m .18629 .10638 L s P p .001 w .18254 .12275 m .18629 .12275 L s P p .001 w .18254 .13912 m .18629 .13912 L s P p .001 w .18254 .15549 m .18629 .15549 L s P p .001 w .18254 .18823 m .18629 .18823 L s P p .001 w .18254 .2046 m .18629 .2046 L s P p .001 w .18254 .22097 m .18629 .22097 L s P p .001 w .18254 .23734 m .18629 .23734 L s P p .001 w .18254 .27008 m .18629 .27008 L s P p .001 w .18254 .28645 m .18629 .28645 L s P p .001 w .18254 .30282 m .18629 .30282 L s P p .001 w .18254 .31919 m .18629 .31919 L s P p .001 w .18254 .35193 m .18629 .35193 L s P p .001 w .18254 .3683 m .18629 .3683 L s P p .001 w .18254 .38467 m .18629 .38467 L s P p .001 w .18254 .40104 m .18629 .40104 L s P p .001 w .18254 .43378 m .18629 .43378 L s P p .001 w .18254 .45015 m .18629 .45015 L s P p .001 w .18254 .46652 m .18629 .46652 L s P p .001 w .18254 .48289 m .18629 .48289 L s P p .001 w .18254 .51564 m .18629 .51564 L s P p .001 w .18254 .53201 m .18629 .53201 L s P p .001 w .18254 .54838 m .18629 .54838 L s P p .001 w .18254 .56475 m .18629 .56475 L s P p .001 w .18254 .59749 m .18629 .59749 L s P p .001 w .18254 .61386 m .18629 .61386 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .1088 L .04365 .18086 L .05357 .23725 L .06349 .28218 L .07341 .31855 L .08333 .3484 L .10317 .39405 L .12302 .42684 L .14286 .45122 L .1627 .46984 L .18254 .48441 L .20238 .49604 L .22222 .50547 L .2619 .51973 L .30159 .52988 L .34127 .5374 L .38095 .54315 L .42063 .54765 L .46032 .55127 L .5 .55422 L .53968 .55668 L .57937 .55874 L .61905 .5605 L .65873 .56201 L .69841 .56333 L .7381 .56448 L .77778 .5655 L .81746 .5664 L .85714 .56721 L .89683 .56793 L .93651 .56859 L .97619 .56918 L s P P p 1 0 0 r p .02 w .32794 .54794 Mdot .09732 .37045 Mdot .44367 .52208 Mdot .04113 .13965 Mdot .18735 .46595 Mdot .03043 .11169 Mdot .78315 .54774 Mdot .03662 .10173 Mdot .9523 .60332 Mdot .53788 .54946 Mdot .42551 .54833 Mdot .23962 .49016 Mdot .12137 .44932 Mdot .16611 .49942 Mdot .65389 .5787 Mdot .14879 .4726 Mdot .49035 .55742 Mdot .14217 .46976 Mdot .68339 .54898 Mdot .12936 .44836 Mdot .70728 .55174 Mdot .56767 .54245 Mdot .30559 .55613 Mdot .35185 .51108 Mdot .20803 .46494 Mdot .20955 .48247 Mdot .53033 .54148 Mdot .08457 .35635 Mdot .06379 .30162 Mdot .91859 .59045 Mdot .35659 .53689 Mdot .81304 .54086 Mdot .6255 .53401 Mdot .26918 .54952 Mdot .34372 .52315 Mdot .89352 .5879 Mdot .1595 .45207 Mdot .70778 .56623 Mdot .60536 .59568 Mdot .64702 .59671 Mdot .56538 .54028 Mdot .70462 .54094 Mdot .2326 .53241 Mdot .86777 .53472 Mdot .5833 .57255 Mdot .6224 .53766 Mdot .26339 .51394 Mdot .70507 .557 Mdot .1277 .40753 Mdot .97348 .59649 Mdot .49853 .52893 Mdot .35027 .56832 Mdot .90933 .59717 Mdot .62184 .57697 Mdot .70054 .5351 Mdot .73026 .57774 Mdot .14106 .44482 Mdot .13167 .44017 Mdot .85386 .53332 Mdot .40279 .5184 Mdot .74997 .59406 Mdot .40551 .52785 Mdot .27526 .54338 Mdot .07905 .32558 Mdot .34211 .56478 Mdot .4334 .55859 Mdot .61776 .58026 Mdot .67934 .59531 Mdot .50051 .54397 Mdot .57148 .55451 Mdot .62284 .5486 Mdot .19249 .47958 Mdot .84906 .56219 Mdot .76318 .53659 Mdot .59761 .57224 Mdot .70794 .58272 Mdot .27931 .54752 Mdot .29834 .50399 Mdot .63774 .52858 Mdot .59519 .53797 Mdot .16104 .4836 Mdot .04753 .22307 Mdot .51439 .54651 Mdot .83122 .56192 Mdot .64152 .59403 Mdot .09186 .3589 Mdot .06772 .28337 Mdot .36011 .54669 Mdot .7646 .56403 Mdot .08558 .35112 Mdot .15066 .45607 Mdot .46658 .53451 Mdot .96581 .58735 Mdot .44286 .55017 Mdot .47523 .54459 Mdot .58783 .52868 Mdot .80989 .56215 Mdot .51978 .57333 Mdot .76598 .54112 Mdot .18348 .46201 Mdot .0252 .0415 Mdot .12171 .4423 Mdot .85073 .54908 Mdot .63132 .55984 Mdot .86111 .56746 Mdot .21227 .53694 Mdot .40969 .5198 Mdot .60063 .5903 Mdot .57599 .57234 Mdot .10466 .37492 Mdot .7862 .56796 Mdot .89737 .53361 Mdot .78481 .59088 Mdot .79947 .55172 Mdot .91027 .57584 Mdot .19196 .51458 Mdot .07298 .28947 Mdot .95588 .59373 Mdot .63947 .59718 Mdot .37906 .57442 Mdot .08729 .34607 Mdot .29821 .54796 Mdot .27728 .54573 Mdot .37703 .56016 Mdot .46866 .58469 Mdot .55374 .55359 Mdot .53458 .54589 Mdot .38559 .55309 Mdot .48541 .53524 Mdot .4059 .56709 Mdot .82213 .54908 Mdot .05065 .24785 Mdot .75865 .5975 Mdot .72864 .53668 Mdot .50518 .53076 Mdot .37542 .53477 Mdot .06032 .24897 Mdot .79787 .59853 Mdot .50193 .58199 Mdot .43609 .53767 Mdot .04033 .16822 Mdot .05399 .24337 Mdot .19439 .48454 Mdot .02715 .06251 Mdot .41192 .54341 Mdot .27471 .53076 Mdot .88294 .55255 Mdot .87548 .54694 Mdot .84642 .53636 Mdot .10142 .37469 Mdot .3683 .52051 Mdot .64152 .59337 Mdot .35178 .54455 Mdot .61134 .5526 Mdot .18121 .46173 Mdot .608 .57442 Mdot .74547 .58245 Mdot .3571 .51235 Mdot .83872 .56421 Mdot .45781 .58056 Mdot .96849 .56155 Mdot .38021 .55376 Mdot .624 .53935 Mdot .71487 .57888 Mdot .29603 .53931 Mdot .12734 .42078 Mdot .13863 .44457 Mdot .49554 .56381 Mdot .36935 .55813 Mdot .16224 .47241 Mdot .50682 .53818 Mdot .68062 .57242 Mdot .51451 .58269 Mdot .32422 .57082 Mdot .86671 .58067 Mdot .58554 .54464 Mdot .59449 .52738 Mdot .48201 .58799 Mdot .47967 .56021 Mdot .96266 .55859 Mdot .13413 .46974 Mdot .82423 .59064 Mdot .6035 .56896 Mdot .16742 .45349 Mdot .1128 .4302 Mdot .81939 .54624 Mdot .22228 .47573 Mdot .25765 .55003 Mdot .60398 .52632 Mdot .75184 .59767 Mdot .14813 .45316 Mdot .76536 .57289 Mdot .03781 .14453 Mdot .91732 .58837 Mdot .4105 .57928 Mdot .77371 .57104 Mdot .32151 .51302 Mdot .93052 .5588 Mdot .12304 .40121 Mdot .69667 .55936 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{-0.0134171, -0.0104952, 0, -0.0104952} :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.042936547557704, A -> 0.989504750389851, R -> 5} ;[o] {a -> 1.04294, A -> 0.989505, R -> 5} :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] angle fixed, fit for Hamaker and radius, the goal is A->1, R->5: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={a->1.05712}; sol=Sort[Join[NonlinearFit[dat4,ptmodel/.cond,d,{R,A}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp4,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-5)/5,(A R-5)/5}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.581089 3.27813 [ [(2.5)] .05026 .58109 0 2 Msboxa [(7.5)] .31481 .58109 0 2 Msboxa [(10)] .44709 .58109 0 2 Msboxa [(12.5)] .57937 .58109 0 2 Msboxa [(15)] .71164 .58109 0 2 Msboxa [(17.5)] .84392 .58109 0 2 Msboxa [(20)] .97619 .58109 0 2 Msboxa [(-0.175)] .17004 .00742 1 0 Msboxa [(-0.15)] .17004 .08937 1 0 Msboxa [(-0.125)] .17004 .17132 1 0 Msboxa [(-0.1)] .17004 .25328 1 0 Msboxa [(-0.075)] .17004 .33523 1 0 Msboxa [(-0.05)] .17004 .41718 1 0 Msboxa [(-0.025)] .17004 .49914 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .58109 m .05026 .58734 L s P [(2.5)] .05026 .58109 0 2 Mshowa p .002 w .31481 .58109 m .31481 .58734 L s P [(7.5)] .31481 .58109 0 2 Mshowa p .002 w .44709 .58109 m .44709 .58734 L s P [(10)] .44709 .58109 0 2 Mshowa p .002 w .57937 .58109 m .57937 .58734 L s P [(12.5)] .57937 .58109 0 2 Mshowa p .002 w .71164 .58109 m .71164 .58734 L s P [(15)] .71164 .58109 0 2 Mshowa p .002 w .84392 .58109 m .84392 .58734 L s P [(17.5)] .84392 .58109 0 2 Mshowa p .002 w .97619 .58109 m .97619 .58734 L s P [(20)] .97619 .58109 0 2 Mshowa p .001 w .07672 .58109 m .07672 .58484 L s P p .001 w .10317 .58109 m .10317 .58484 L s P p .001 w .12963 .58109 m .12963 .58484 L s P p .001 w .15608 .58109 m .15608 .58484 L s P p .001 w .20899 .58109 m .20899 .58484 L s P p .001 w .23545 .58109 m .23545 .58484 L s P p .001 w .2619 .58109 m .2619 .58484 L s P p .001 w .28836 .58109 m .28836 .58484 L s P p .001 w .34127 .58109 m .34127 .58484 L s P p .001 w .36772 .58109 m .36772 .58484 L s P p .001 w .39418 .58109 m .39418 .58484 L s P p .001 w .42063 .58109 m .42063 .58484 L s P p .001 w .47354 .58109 m .47354 .58484 L s P p .001 w .5 .58109 m .5 .58484 L s P p .001 w .52646 .58109 m .52646 .58484 L s P p .001 w .55291 .58109 m .55291 .58484 L s P p .001 w .60582 .58109 m .60582 .58484 L s P p .001 w .63228 .58109 m .63228 .58484 L s P p .001 w .65873 .58109 m .65873 .58484 L s P p .001 w .68519 .58109 m .68519 .58484 L s P p .001 w .7381 .58109 m .7381 .58484 L s P p .001 w .76455 .58109 m .76455 .58484 L s P p .001 w .79101 .58109 m .79101 .58484 L s P p .001 w .81746 .58109 m .81746 .58484 L s P p .001 w .87037 .58109 m .87037 .58484 L s P p .001 w .89683 .58109 m .89683 .58484 L s P p .001 w .92328 .58109 m .92328 .58484 L s P p .001 w .94974 .58109 m .94974 .58484 L s P p .001 w .02381 .58109 m .02381 .58484 L s P p .002 w 0 .58109 m 1 .58109 L s P p .002 w .18254 .00742 m .18879 .00742 L s P [(-0.175)] .17004 .00742 1 0 Mshowa p .002 w .18254 .08937 m .18879 .08937 L s P [(-0.15)] .17004 .08937 1 0 Mshowa p .002 w .18254 .17132 m .18879 .17132 L s P [(-0.125)] .17004 .17132 1 0 Mshowa p .002 w .18254 .25328 m .18879 .25328 L s P [(-0.1)] .17004 .25328 1 0 Mshowa p .002 w .18254 .33523 m .18879 .33523 L s P [(-0.075)] .17004 .33523 1 0 Mshowa p .002 w .18254 .41718 m .18879 .41718 L s P [(-0.05)] .17004 .41718 1 0 Mshowa p .002 w .18254 .49914 m .18879 .49914 L s P [(-0.025)] .17004 .49914 1 0 Mshowa p .001 w .18254 .02381 m .18629 .02381 L s P p .001 w .18254 .0402 m .18629 .0402 L s P p .001 w .18254 .05659 m .18629 .05659 L s P p .001 w .18254 .07298 m .18629 .07298 L s P p .001 w .18254 .10576 m .18629 .10576 L s P p .001 w .18254 .12215 m .18629 .12215 L s P p .001 w .18254 .13854 m .18629 .13854 L s P p .001 w .18254 .15493 m .18629 .15493 L s P p .001 w .18254 .18771 m .18629 .18771 L s P p .001 w .18254 .2041 m .18629 .2041 L s P p .001 w .18254 .22049 m .18629 .22049 L s P p .001 w .18254 .23689 m .18629 .23689 L s P p .001 w .18254 .26967 m .18629 .26967 L s P p .001 w .18254 .28606 m .18629 .28606 L s P p .001 w .18254 .30245 m .18629 .30245 L s P p .001 w .18254 .31884 m .18629 .31884 L s P p .001 w .18254 .35162 m .18629 .35162 L s P p .001 w .18254 .36801 m .18629 .36801 L s P p .001 w .18254 .3844 m .18629 .3844 L s P p .001 w .18254 .40079 m .18629 .40079 L s P p .001 w .18254 .43357 m .18629 .43357 L s P p .001 w .18254 .44996 m .18629 .44996 L s P p .001 w .18254 .46635 m .18629 .46635 L s P p .001 w .18254 .48275 m .18629 .48275 L s P p .001 w .18254 .51553 m .18629 .51553 L s P p .001 w .18254 .53192 m .18629 .53192 L s P p .001 w .18254 .54831 m .18629 .54831 L s P p .001 w .18254 .5647 m .18629 .5647 L s P p .001 w .18254 .59748 m .18629 .59748 L s P p .001 w .18254 .61387 m .18629 .61387 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .10912 L .04365 .18138 L .05357 .2379 L .06349 .28292 L .07341 .31934 L .08333 .34923 L .10317 .39489 L .12302 .42768 L .14286 .45203 L .1627 .47062 L .18254 .48515 L .20238 .49674 L .22222 .50614 L .2619 .52034 L .30159 .53044 L .34127 .5379 L .38095 .54361 L .42063 .54808 L .46032 .55166 L .5 .55459 L .53968 .55702 L .57937 .55906 L .61905 .5608 L .65873 .56229 L .69841 .56359 L .7381 .56473 L .77778 .56573 L .81746 .56663 L .85714 .56742 L .89683 .56814 L .93651 .56878 L .97619 .56937 L s P P p 1 0 0 r p .02 w .32794 .54787 Mdot .09732 .37016 Mdot .44367 .52198 Mdot .04113 .13907 Mdot .18735 .46578 Mdot .03043 .11108 Mdot .78315 .54767 Mdot .03662 .1011 Mdot .9523 .60332 Mdot .53788 .5494 Mdot .42551 .54826 Mdot .23962 .49002 Mdot .12137 .44913 Mdot .16611 .49929 Mdot .65389 .57867 Mdot .14879 .47244 Mdot .49035 .55736 Mdot .14217 .4696 Mdot .68339 .54892 Mdot .12936 .44816 Mdot .70728 .55167 Mdot .56767 .54238 Mdot .30559 .55607 Mdot .35185 .51097 Mdot .20803 .46476 Mdot .20955 .48232 Mdot .53033 .54141 Mdot .08457 .35604 Mdot .06379 .30124 Mdot .91859 .59044 Mdot .35659 .5368 Mdot .81304 .54078 Mdot .6255 .53392 Mdot .26918 .54945 Mdot .34372 .52305 Mdot .89352 .58788 Mdot .1595 .45188 Mdot .70778 .56619 Mdot .60536 .59567 Mdot .64702 .5967 Mdot .56538 .54021 Mdot .70462 .54086 Mdot .2326 .53232 Mdot .86777 .53464 Mdot .5833 .57251 Mdot .6224 .53758 Mdot .26339 .51383 Mdot .70507 .55694 Mdot .1277 .40729 Mdot .97348 .59648 Mdot .49853 .52883 Mdot .35027 .56828 Mdot .90933 .59717 Mdot .62184 .57693 Mdot .70054 .53502 Mdot .73026 .57771 Mdot .14106 .44463 Mdot .13167 .43997 Mdot .85386 .53323 Mdot .40279 .51829 Mdot .74997 .59405 Mdot .40551 .52775 Mdot .27526 .5433 Mdot .07905 .32524 Mdot .34211 .56473 Mdot .4334 .55854 Mdot .61776 .58023 Mdot .67934 .5953 Mdot .50051 .5439 Mdot .57148 .55445 Mdot .62284 .54853 Mdot .19249 .47943 Mdot .84906 .56214 Mdot .76318 .53651 Mdot .59761 .5722 Mdot .70794 .5827 Mdot .27931 .54745 Mdot .29834 .50387 Mdot .63774 .52849 Mdot .59519 .53789 Mdot .16104 .48345 Mdot .04753 .2226 Mdot .51439 .54644 Mdot .83122 .56187 Mdot .64152 .59402 Mdot .09186 .35859 Mdot .06772 .28297 Mdot .36011 .54662 Mdot .7646 .56398 Mdot .08558 .35081 Mdot .15066 .45588 Mdot .46658 .53443 Mdot .96581 .58733 Mdot .44286 .5501 Mdot .47523 .54452 Mdot .58783 .52859 Mdot .80989 .5621 Mdot .51978 .57329 Mdot .76598 .54104 Mdot .18348 .46184 Mdot .0252 .0408 Mdot .12171 .4421 Mdot .85073 .54901 Mdot .63132 .55978 Mdot .86111 .56742 Mdot .21227 .53686 Mdot .40969 .5197 Mdot .60063 .59028 Mdot .57599 .5723 Mdot .10466 .37464 Mdot .7862 .56792 Mdot .89737 .53353 Mdot .78481 .59086 Mdot .79947 .55165 Mdot .91027 .5758 Mdot .19196 .51447 Mdot .07298 .28908 Mdot .95588 .59372 Mdot .63947 .59717 Mdot .37906 .57438 Mdot .08729 .34575 Mdot .29821 .54789 Mdot .27728 .54566 Mdot .37703 .56011 Mdot .46866 .58467 Mdot .55374 .55352 Mdot .53458 .54582 Mdot .38559 .55303 Mdot .48541 .53516 Mdot .4059 .56705 Mdot .82213 .54901 Mdot .05065 .24741 Mdot .75865 .59749 Mdot .72864 .5366 Mdot .50518 .53067 Mdot .37542 .53469 Mdot .06032 .24853 Mdot .79787 .59852 Mdot .50193 .58196 Mdot .43609 .53758 Mdot .04033 .16768 Mdot .05399 .24293 Mdot .19439 .4844 Mdot .02715 .06184 Mdot .41192 .54334 Mdot .27471 .53067 Mdot .88294 .55248 Mdot .87548 .54687 Mdot .84642 .53628 Mdot .10142 .37441 Mdot .3683 .52041 Mdot .64152 .59336 Mdot .35178 .54448 Mdot .61134 .55254 Mdot .18121 .46155 Mdot .608 .57439 Mdot .74547 .58242 Mdot .3571 .51224 Mdot .83872 .56416 Mdot .45781 .58053 Mdot .96849 .56149 Mdot .38021 .5537 Mdot .624 .53927 Mdot .71487 .57885 Mdot .29603 .53923 Mdot .12734 .42056 Mdot .13863 .44437 Mdot .49554 .56376 Mdot .36935 .55808 Mdot .16224 .47224 Mdot .50682 .53809 Mdot .68062 .57238 Mdot .51451 .58266 Mdot .32422 .57078 Mdot .86671 .58064 Mdot .58554 .54456 Mdot .59449 .52729 Mdot .48201 .58797 Mdot .47967 .56016 Mdot .96266 .55853 Mdot .13413 .46958 Mdot .82423 .59062 Mdot .6035 .56892 Mdot .16742 .45331 Mdot .1128 .42999 Mdot .81939 .54616 Mdot .22228 .47557 Mdot .25765 .54997 Mdot .60398 .52622 Mdot .75184 .59767 Mdot .14813 .45298 Mdot .76536 .57285 Mdot .03781 .14396 Mdot .91732 .58836 Mdot .4105 .57925 Mdot .77371 .571 Mdot .32151 .51291 Mdot .93052 .55875 Mdot .12304 .40096 Mdot .69667 .5593 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{0., 0.031833, -0.032752, -0.00196161} :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.05712, A -> 1.031832964279692, R -> 4.836240078073561} ;[o] {a -> 1.05712, A -> 1.03183, R -> 4.83624} :[font = subsubsection; inactive; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; bold; fontName = "Arial"; fontSize = 10; ] now everything free: :[font = input; startGroup; nowordwrap; backColorRed = 65280; backColorGreen = 65280; backColorBlue = 65280; fontColorRed = 0; fontColorGreen = 0; fontColorBlue = 0; plain; fontName = "Courier New"; fontSize = 10; ] Module[{cond,sol}, cond={}; sol=Sort[Join[NonlinearFit[dat4,ptmodel/.cond,d,{a,{R,5},A}],cond]]; Show[Plot[ptmodel/.sol,{d,2,20},DisplayFunction->Identity], lp4,PlotRange->All,DisplayFunction->$DisplayFunction]; Print["Errors:", {(a-1.05712)/1.05712,(A-1),(R-5)/5,(A R-5)/5}/.sol]; sol] :[font = postscript; inactive; output; BITMAP; PostScript; pictureLeft = 100; pictureTop = 0; pictureWidth = 300; pictureHeight = 184; nowordwrap; ] %! %%Creator: Mathematica %%AspectRatio: .61803 MathPictureStart %% Graphics /Courier findfont 10 scalefont setfont % Scaling calculations -0.0820106 0.0529101 0.581912 3.20264 [ [(2.5)] .05026 .58191 0 2 Msboxa [(7.5)] .31481 .58191 0 2 Msboxa [(10)] .44709 .58191 0 2 Msboxa [(12.5)] .57937 .58191 0 2 Msboxa [(15)] .71164 .58191 0 2 Msboxa [(17.5)] .84392 .58191 0 2 Msboxa [(20)] .97619 .58191 0 2 Msboxa [(-0.175)] .17004 .02145 1 0 Msboxa [(-0.15)] .17004 .10152 1 0 Msboxa [(-0.125)] .17004 .18158 1 0 Msboxa [(-0.1)] .17004 .26165 1 0 Msboxa [(-0.075)] .17004 .34171 1 0 Msboxa [(-0.05)] .17004 .42178 1 0 Msboxa [(-0.025)] .17004 .50185 1 0 Msboxa [ -0.001 -0.001 0 0 ] [ 1.001 .61903 0 0 ] ] MathScale % Start of Graphics 1 setlinecap 1 setlinejoin newpath [ ] 0 setdash 0 g p p .002 w .05026 .58191 m .05026 .58816 L s P [(2.5)] .05026 .58191 0 2 Mshowa p .002 w .31481 .58191 m .31481 .58816 L s P [(7.5)] .31481 .58191 0 2 Mshowa p .002 w .44709 .58191 m .44709 .58816 L s P [(10)] .44709 .58191 0 2 Mshowa p .002 w .57937 .58191 m .57937 .58816 L s P [(12.5)] .57937 .58191 0 2 Mshowa p .002 w .71164 .58191 m .71164 .58816 L s P [(15)] .71164 .58191 0 2 Mshowa p .002 w .84392 .58191 m .84392 .58816 L s P [(17.5)] .84392 .58191 0 2 Mshowa p .002 w .97619 .58191 m .97619 .58816 L s P [(20)] .97619 .58191 0 2 Mshowa p .001 w .07672 .58191 m .07672 .58566 L s P p .001 w .10317 .58191 m .10317 .58566 L s P p .001 w .12963 .58191 m .12963 .58566 L s P p .001 w .15608 .58191 m .15608 .58566 L s P p .001 w .20899 .58191 m .20899 .58566 L s P p .001 w .23545 .58191 m .23545 .58566 L s P p .001 w .2619 .58191 m .2619 .58566 L s P p .001 w .28836 .58191 m .28836 .58566 L s P p .001 w .34127 .58191 m .34127 .58566 L s P p .001 w .36772 .58191 m .36772 .58566 L s P p .001 w .39418 .58191 m .39418 .58566 L s P p .001 w .42063 .58191 m .42063 .58566 L s P p .001 w .47354 .58191 m .47354 .58566 L s P p .001 w .5 .58191 m .5 .58566 L s P p .001 w .52646 .58191 m .52646 .58566 L s P p .001 w .55291 .58191 m .55291 .58566 L s P p .001 w .60582 .58191 m .60582 .58566 L s P p .001 w .63228 .58191 m .63228 .58566 L s P p .001 w .65873 .58191 m .65873 .58566 L s P p .001 w .68519 .58191 m .68519 .58566 L s P p .001 w .7381 .58191 m .7381 .58566 L s P p .001 w .76455 .58191 m .76455 .58566 L s P p .001 w .79101 .58191 m .79101 .58566 L s P p .001 w .81746 .58191 m .81746 .58566 L s P p .001 w .87037 .58191 m .87037 .58566 L s P p .001 w .89683 .58191 m .89683 .58566 L s P p .001 w .92328 .58191 m .92328 .58566 L s P p .001 w .94974 .58191 m .94974 .58566 L s P p .001 w .02381 .58191 m .02381 .58566 L s P p .002 w 0 .58191 m 1 .58191 L s P p .002 w .18254 .02145 m .18879 .02145 L s P [(-0.175)] .17004 .02145 1 0 Mshowa p .002 w .18254 .10152 m .18879 .10152 L s P [(-0.15)] .17004 .10152 1 0 Mshowa p .002 w .18254 .18158 m .18879 .18158 L s P [(-0.125)] .17004 .18158 1 0 Mshowa p .002 w .18254 .26165 m .18879 .26165 L s P [(-0.1)] .17004 .26165 1 0 Mshowa p .002 w .18254 .34171 m .18879 .34171 L s P [(-0.075)] .17004 .34171 1 0 Mshowa p .002 w .18254 .42178 m .18879 .42178 L s P [(-0.05)] .17004 .42178 1 0 Mshowa p .002 w .18254 .50185 m .18879 .50185 L s P [(-0.025)] .17004 .50185 1 0 Mshowa p .001 w .18254 .03746 m .18629 .03746 L s P p .001 w .18254 .05348 m .18629 .05348 L s P p .001 w .18254 .06949 m .18629 .06949 L s P p .001 w .18254 .0855 m .18629 .0855 L s P p .001 w .18254 .11753 m .18629 .11753 L s P p .001 w .18254 .13354 m .18629 .13354 L s P p .001 w .18254 .14955 m .18629 .14955 L s P p .001 w .18254 .16557 m .18629 .16557 L s P p .001 w .18254 .19759 m .18629 .19759 L s P p .001 w .18254 .21361 m .18629 .21361 L s P p .001 w .18254 .22962 m .18629 .22962 L s P p .001 w .18254 .24563 m .18629 .24563 L s P p .001 w .18254 .27766 m .18629 .27766 L s P p .001 w .18254 .29367 m .18629 .29367 L s P p .001 w .18254 .30969 m .18629 .30969 L s P p .001 w .18254 .3257 m .18629 .3257 L s P p .001 w .18254 .35773 m .18629 .35773 L s P p .001 w .18254 .37374 m .18629 .37374 L s P p .001 w .18254 .38975 m .18629 .38975 L s P p .001 w .18254 .40577 m .18629 .40577 L s P p .001 w .18254 .43779 m .18629 .43779 L s P p .001 w .18254 .45381 m .18629 .45381 L s P p .001 w .18254 .46982 m .18629 .46982 L s P p .001 w .18254 .48583 m .18629 .48583 L s P p .001 w .18254 .51786 m .18629 .51786 L s P p .001 w .18254 .53387 m .18629 .53387 L s P p .001 w .18254 .54989 m .18629 .54989 L s P p .001 w .18254 .5659 m .18629 .5659 L s P p .001 w .18254 .00544 m .18629 .00544 L s P p .001 w .18254 .59792 m .18629 .59792 L s P p .001 w .18254 .61394 m .18629 .61394 L s P p .002 w .18254 0 m .18254 .61803 L s P P p p p .004 w .02381 .01472 m .03373 .10987 L .04365 .18269 L .05357 .2396 L .06349 .2849 L .07341 .32153 L .08333 .35157 L .10317 .39743 L .12302 .4303 L .14286 .45468 L .1627 .47326 L .18254 .48777 L .20238 .49932 L .22222 .50867 L .2619 .52277 L .30159 .53277 L .34127 .54014 L .38095 .54575 L .42063 .55015 L .46032 .55366 L .5 .55652 L .53968 .55889 L .57937 .56087 L .61905 .56256 L .65873 .56401 L .69841 .56527 L .7381 .56637 L .77778 .56734 L .81746 .5682 L .85714 .56897 L .89683 .56965 L .93651 .57028 L .97619 .57084 L s P P p 1 0 0 r p .02 w .66537 .57357 Mdot .89128 .54357 Mdot .26986 .4978 Mdot .65605 .54874 Mdot .16071 .46714 Mdot .87536 .54517 Mdot .4716 .53945 Mdot .35894 .54334 Mdot .88275 .60332 Mdot .04525 .17446 Mdot .48369 .55282 Mdot .12532 .41218 Mdot .79451 .59085 Mdot .21022 .47124 Mdot .54846 .57477 Mdot .53036 .57595 Mdot .41156 .5351 Mdot .6312 .56935 Mdot .91615 .53552 Mdot .29606 .50243 Mdot .05721 .2302 Mdot .27462 .55835 Mdot .54971 .56392 Mdot .17311 .50206 Mdot .7314 .58263 Mdot .93908 .58734 Mdot .20675 .46797 Mdot .43253 .58591 Mdot .77138 .58219 Mdot .77752 .56022 Mdot .83141 .58257 Mdot .50527 .54383 Mdot .79801 .57977 Mdot .25446 .51132 Mdot .27211 .49817 Mdot .10516 .40976 Mdot .5169 .54841 Mdot .14227 .47957 Mdot .33396 .52888 Mdot .68593 .58933 Mdot .53877 .56871 Mdot .8846 .56564 Mdot .68355 .59568 Mdot .40313 .52232 Mdot .86173 .55257 Mdot .17248 .50056 Mdot .61343 .53836 Mdot .09113 .35234 Mdot .12034 .44895 Mdot .92505 .55873 Mdot .76258 .56218 Mdot .26293 .52276 Mdot .24762 .53919 Mdot .35452 .51008 Mdot .54026 .58362 Mdot .92758 .5991 Mdot .65472 .56968 Mdot .77891 .54236 Mdot .0651 .29102 Mdot .71159 .59832 Mdot .92095 .57925 Mdot .76273 .57568 Mdot .18218 .46718 Mdot .52361 .53205 Mdot .91075 .56319 Mdot .19289 .46592 Mdot .3943 .55117 Mdot .24151 .49376 Mdot .71577 .59878 Mdot .43879 .51864 Mdot .67448 .56776 Mdot .70339 .53458 Mdot .72973 .59527 Mdot .91685 .56555 Mdot .57136 .57562 Mdot .41706 .53604 Mdot .6368 .54968 Mdot .24797 .49873 Mdot .2663 .54225 Mdot .03027 .05754 Mdot .52672 .57809 Mdot .56768 .52567 Mdot .82843 .54974 Mdot .84048 .59992 Mdot .12251 .419 Mdot .89981 .57005 Mdot .52734 .56551 Mdot .50657 .53579 Mdot .86673 .55464 Mdot .28241 .56271 Mdot .62424 .56555 Mdot .27594 .5583 Mdot .12133 .44802 Mdot .68446 .59543 Mdot .26909 .52651 Mdot .82017 .60266 Mdot .1704 .45908 Mdot .89655 .56836 Mdot .61925 .56902 Mdot .41379 .53407 Mdot .72871 .55171 Mdot .15519 .45374 Mdot .12829 .45734 Mdot .85544 .55835 Mdot .03077 .12955 Mdot .19479 .48998 Mdot .73787 .56474 Mdot .35081 .53642 Mdot .59128 .54579 Mdot .43045 .58218 Mdot .94822 .58348 Mdot .04047 .18085 Mdot .24332 .50938 Mdot .86147 .55924 Mdot .13884 .46125 Mdot .02984 .05163 Mdot .13188 .45325 Mdot .81124 .57466 Mdot .3702 .5201 Mdot .48424 .53209 Mdot .75511 .59312 Mdot .0776 .31982 Mdot .78308 .5454 Mdot .06094 .28447 Mdot .56357 .5811 Mdot .17566 .46328 Mdot .44854 .52909 Mdot .16964 .45819 Mdot .34047 .5423 Mdot .33459 .54604 Mdot .94646 .56157 Mdot .82654 .56323 Mdot .21516 .50615 Mdot .77275 .54603 Mdot .40827 .5356 Mdot .73562 .57557 Mdot .82088 .56903 Mdot .58376 .5559 Mdot .39615 .53673 Mdot .43794 .53299 Mdot .07949 .36675 Mdot .12716 .44741 Mdot .10922 .40669 Mdot .27681 .50649 Mdot .87025 .59993 Mdot .48025 .52027 Mdot .4858 .56434 Mdot .72082 .5568 Mdot .6411 .53603 Mdot .16087 .50419 Mdot .26876 .54378 Mdot .69912 .57995 Mdot .21309 .53699 Mdot .59578 .53147 Mdot .12768 .43943 Mdot .34278 .56628 Mdot .23361 .51897 Mdot .83872 .59392 Mdot .72401 .59751 Mdot .14171 .45413 Mdot .10671 .43027 Mdot .95703 .57356 Mdot .81414 .53763 Mdot .28171 .55191 Mdot .62486 .5357 Mdot .66212 .58234 Mdot .521 .56506 Mdot .34315 .5682 Mdot .31119 .50169 Mdot .48063 .52216 Mdot .56338 .53041 Mdot .36273 .54647 Mdot .48047 .53708 Mdot .38189 .58108 Mdot .64252 .54045 Mdot .12399 .41006 Mdot .04147 .1452 Mdot .43806 .54505 Mdot .49666 .55296 Mdot .11872 .42697 Mdot .20928 .49808 Mdot .61428 .56081 Mdot .6221 .55161 Mdot .27536 .55744 Mdot .16543 .44839 Mdot .86966 .59976 Mdot .4991 .5863 Mdot .76947 .58604 Mdot .48143 .58101 Mdot .35522 .53277 Mdot .96096 .56436 Mdot .26031 .54743 Mdot .77549 .60155 Mdot .62222 .55362 Mdot .1772 .49522 Mdot .37067 .54004 Mdot .03558 .13715 Mdot .47721 .55499 Mdot .51268 .56554 Mdot .68393 .58104 Mdot P P P 0 0 m 1 0 L 1 .61803 L 0 .61803 L closepath clip newpath % End of Graphics MathPictureEnd :[font = print; inactive; formatted; output; nowordwrap; ] No Input Form Generated ;[o] Errors:{0.0130051, 0.0444899, -0.0134947, 0.0303948} :[font = output; inactive; formatted; output; endGroup; nowordwrap; ] {a -> 1.070867915525963, A -> 1.044489906992488, R -> 4.93252643345754} ;[o] {a -> 1.07087, A -> 1.04449, R -> 4.93253} ^*)