Problem in using FindFit












4














I have the following set of data



data={{0,0,0},{0,2,1},{0,4,2.247},{0,6,3.627},{0,8,5.031},{1,0,3.346}};


where the values are {n, L,$varepsilon$} and satisfy the following equations



$E(n,L) = 2n+1 + sqrt{L(L+1)-frac{3}{4}(L)^2 + 1 + beta_0^4}$



e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]


$varepsilon = frac{E(n,L)-E(0,0)}{E(0,2)-E(0,0)}$,



where $beta_0$ should be determined. I don't know how I can use FindFit command of Mathematica to find the best value of $beta_0$ to have the best fit for $varepsilon$.






fitting







share|improve this question





























    4














    I have the following set of data



    data={{0,0,0},{0,2,1},{0,4,2.247},{0,6,3.627},{0,8,5.031},{1,0,3.346}};


    where the values are {n, L,$varepsilon$} and satisfy the following equations



    $E(n,L) = 2n+1 + sqrt{L(L+1)-frac{3}{4}(L)^2 + 1 + beta_0^4}$



    e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]


    $varepsilon = frac{E(n,L)-E(0,0)}{E(0,2)-E(0,0)}$,



    where $beta_0$ should be determined. I don't know how I can use FindFit command of Mathematica to find the best value of $beta_0$ to have the best fit for $varepsilon$.






    fitting







    share|improve this question



























      4












      4








      4







      I have the following set of data



      data={{0,0,0},{0,2,1},{0,4,2.247},{0,6,3.627},{0,8,5.031},{1,0,3.346}};


      where the values are {n, L,$varepsilon$} and satisfy the following equations



      $E(n,L) = 2n+1 + sqrt{L(L+1)-frac{3}{4}(L)^2 + 1 + beta_0^4}$



      e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]


      $varepsilon = frac{E(n,L)-E(0,0)}{E(0,2)-E(0,0)}$,



      where $beta_0$ should be determined. I don't know how I can use FindFit command of Mathematica to find the best value of $beta_0$ to have the best fit for $varepsilon$.






      fitting







      share|improve this question















      I have the following set of data



      data={{0,0,0},{0,2,1},{0,4,2.247},{0,6,3.627},{0,8,5.031},{1,0,3.346}};


      where the values are {n, L,$varepsilon$} and satisfy the following equations



      $E(n,L) = 2n+1 + sqrt{L(L+1)-frac{3}{4}(L)^2 + 1 + beta_0^4}$



      e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]


      $varepsilon = frac{E(n,L)-E(0,0)}{E(0,2)-E(0,0)}$,



      where $beta_0$ should be determined. I don't know how I can use FindFit command of Mathematica to find the best value of $beta_0$ to have the best fit for $varepsilon$.






      fitting




      fitting






      share|improve this question















      share|improve this question













      share|improve this question




      share|improve this question








      edited Nov 12 '18 at 7:19









      Coolwater

      14.7k32553




      14.7k32553










      asked Nov 12 '18 at 6:13









      Hadi Sobhani

      37327




      37327






















          2 Answers
          2






          active

          oldest

          votes


















          6














          e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]
          

          
FindFit[data, (e[n, L] - e[0, 0])/(e[0, 2] - e[0, 0]), b0, {n, L}]



          {b0 -> 1.3514967}




          Which seems reasonable in view of the residuals:



          Plot[Evaluate[(e[#, #2] - e[0, 0])/(e[0, 2] - e[0, 0]) - #3 & @@@ data], {b0, 0, 3}]




          The brown and purple residual has bigger slope around the roots in the plots. Hence for Mathematica to minimize the sum of squares in the y-dimension, the mean of the 2 data points that correspond to the big slopes are cared more about than the others. It is purpose specific whether this is appropriate. If it isn't you can add the NormFunction-option to FindFit.






          share|improve this answer























          • Thank you dear @Coolwater. How does Mathematica recognize {n,L} for each data?
            – Hadi Sobhani
            Nov 12 '18 at 7:25










          • Also, how does it recognize that the expr is the 3rd value of every data element?
            – J42161217
            Nov 12 '18 at 7:34






          • 1




            FindFit assumes its first argument has the form {{var1, var2, ..., varN, expr}, ... , {var1, var2, ..., varN, expr}} where {var1, var2, ..., varN} is the 4th argument of FindFit
            – Coolwater
            Nov 12 '18 at 7:36










          • Ok! Given b -> {1.27225, 1.29505, 1.28573, 1.40411} having Mean=1.31428 and Medean=1.29039 do you think Mathematica did a good job? Anyway +1 from me
            – J42161217
            Nov 12 '18 at 7:52










          • @J42161217 It uses least squares, see edit
            – Coolwater
            Nov 12 '18 at 8:05



















          2














          You can also use NMinimize. First we need to write cost function, i.e. residual.



          data = {{0, 0, 0}, {0, 2, 1}, {0, 4, 2.247}, {0, 6, 3.627}, {0, 8, 
          
    5.031}, {1, 0, 3.346}};
          
e[n_, L_] := 2 n + 1 + Sqrt[L (L + 1) - 3/4 L^2 + 1 + b0^4]
          
cost[b0_] =Sum[(e @@data[[i, 1 ;; 2]] - (data[[i, 3]] (e[0, 2] - e[0, 0]) + 
          
      e[0, 0]))^2, {i, 6}];
          
(*or Total[(e[#1, #2] - (#3 (e[0, 2] - e[0, 0]) + e[0, 0]))^2 & @@@ data]*)
          

          
fit = NMinimize[cost[b0] , b0]



          {0.0196376, {b0 -> 1.35462}}




          Since your cost function has only one variable you can also use grid search.



          Ordering[val,1] gives position of min value.



          b0Val = Range[0, 10, 0.0001];
          
val = cost[b0Val];
          
b0Val[[Ordering[val, 1]]]



          {1.3546}




          Note that there is another min at b0=-1.3546



          b0Val = Range[-1000, 1000, 0.001];    
          
val = cost[b0Val];    
          
b0Val[[Ordering[val, 2]]]



          {-1.3546, 1.3546}




          We can plot cost function



          $text{cost}(b0)=left(-5.031 left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+25}right)^2\+left(-3.627
          left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+
          sqrt{text{b0}^4+16}right)^2\+left(2-3.346
          left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)right)^2+left(-2.247
          left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+9}right)^2$



          Plot[cost[b0], {b0, -10, 10}]


          enter image description here






          share|improve this answer























            Your Answer





            StackExchange.ifUsing("editor", function () {
            return StackExchange.using("mathjaxEditing", function () {
            StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix) {
            StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
            });
            });
            }, "mathjax-editing");

            StackExchange.ready(function() {
            var channelOptions = {
            tags: "".split(" "),
            id: "387"
            };
            initTagRenderer("".split(" "), "".split(" "), channelOptions);

            StackExchange.using("externalEditor", function() {
            // Have to fire editor after snippets, if snippets enabled
            if (StackExchange.settings.snippets.snippetsEnabled) {
            StackExchange.using("snippets", function() {
            createEditor();
            });
            }
            else {
            createEditor();
            }
            });

            function createEditor() {
            StackExchange.prepareEditor({
            heartbeatType: 'answer',
            autoActivateHeartbeat: false,
            convertImagesToLinks: false,
            noModals: true,
            showLowRepImageUploadWarning: true,
            reputationToPostImages: null,
            bindNavPrevention: true,
            postfix: "",
            imageUploader: {
            brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
            contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
            allowUrls: true
            },
            onDemand: true,
            discardSelector: ".discard-answer"
            ,immediatelyShowMarkdownHelp:true
            });


            }
            });














            draft saved

            draft discarded


















            StackExchange.ready(
            function () {
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f185829%2fproblem-in-using-findfit%23new-answer', 'question_page');
            }
            );

            Post as a guest















            Required, but never shown

























            2 Answers
            2






            active

            oldest

            votes








            2 Answers
            2






            active

            oldest

            votes









            active

            oldest

            votes






            active

            oldest

            votes









            6














            e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]
            

            
FindFit[data, (e[n, L] - e[0, 0])/(e[0, 2] - e[0, 0]), b0, {n, L}]



            {b0 -> 1.3514967}




            Which seems reasonable in view of the residuals:



            Plot[Evaluate[(e[#, #2] - e[0, 0])/(e[0, 2] - e[0, 0]) - #3 & @@@ data], {b0, 0, 3}]




            The brown and purple residual has bigger slope around the roots in the plots. Hence for Mathematica to minimize the sum of squares in the y-dimension, the mean of the 2 data points that correspond to the big slopes are cared more about than the others. It is purpose specific whether this is appropriate. If it isn't you can add the NormFunction-option to FindFit.






            share|improve this answer























            • Thank you dear @Coolwater. How does Mathematica recognize {n,L} for each data?
              – Hadi Sobhani
              Nov 12 '18 at 7:25










            • Also, how does it recognize that the expr is the 3rd value of every data element?
              – J42161217
              Nov 12 '18 at 7:34






            • 1




              FindFit assumes its first argument has the form {{var1, var2, ..., varN, expr}, ... , {var1, var2, ..., varN, expr}} where {var1, var2, ..., varN} is the 4th argument of FindFit
              – Coolwater
              Nov 12 '18 at 7:36










            • Ok! Given b -> {1.27225, 1.29505, 1.28573, 1.40411} having Mean=1.31428 and Medean=1.29039 do you think Mathematica did a good job? Anyway +1 from me
              – J42161217
              Nov 12 '18 at 7:52










            • @J42161217 It uses least squares, see edit
              – Coolwater
              Nov 12 '18 at 8:05
















            6














            e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]
            

            
FindFit[data, (e[n, L] - e[0, 0])/(e[0, 2] - e[0, 0]), b0, {n, L}]



            {b0 -> 1.3514967}




            Which seems reasonable in view of the residuals:



            Plot[Evaluate[(e[#, #2] - e[0, 0])/(e[0, 2] - e[0, 0]) - #3 & @@@ data], {b0, 0, 3}]




            The brown and purple residual has bigger slope around the roots in the plots. Hence for Mathematica to minimize the sum of squares in the y-dimension, the mean of the 2 data points that correspond to the big slopes are cared more about than the others. It is purpose specific whether this is appropriate. If it isn't you can add the NormFunction-option to FindFit.






            share|improve this answer























            • Thank you dear @Coolwater. How does Mathematica recognize {n,L} for each data?
              – Hadi Sobhani
              Nov 12 '18 at 7:25










            • Also, how does it recognize that the expr is the 3rd value of every data element?
              – J42161217
              Nov 12 '18 at 7:34






            • 1




              FindFit assumes its first argument has the form {{var1, var2, ..., varN, expr}, ... , {var1, var2, ..., varN, expr}} where {var1, var2, ..., varN} is the 4th argument of FindFit
              – Coolwater
              Nov 12 '18 at 7:36










            • Ok! Given b -> {1.27225, 1.29505, 1.28573, 1.40411} having Mean=1.31428 and Medean=1.29039 do you think Mathematica did a good job? Anyway +1 from me
              – J42161217
              Nov 12 '18 at 7:52










            • @J42161217 It uses least squares, see edit
              – Coolwater
              Nov 12 '18 at 8:05














            6












            6








            6






            e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]
            

            
FindFit[data, (e[n, L] - e[0, 0])/(e[0, 2] - e[0, 0]), b0, {n, L}]



            {b0 -> 1.3514967}




            Which seems reasonable in view of the residuals:



            Plot[Evaluate[(e[#, #2] - e[0, 0])/(e[0, 2] - e[0, 0]) - #3 & @@@ data], {b0, 0, 3}]




            The brown and purple residual has bigger slope around the roots in the plots. Hence for Mathematica to minimize the sum of squares in the y-dimension, the mean of the 2 data points that correspond to the big slopes are cared more about than the others. It is purpose specific whether this is appropriate. If it isn't you can add the NormFunction-option to FindFit.






            share|improve this answer














            e[n_, L_] = 2n + 1 + Sqrt[L(L + 1) - 3/4 L^2 + 1 + b0^4]
            

            
FindFit[data, (e[n, L] - e[0, 0])/(e[0, 2] - e[0, 0]), b0, {n, L}]



            {b0 -> 1.3514967}




            Which seems reasonable in view of the residuals:



            Plot[Evaluate[(e[#, #2] - e[0, 0])/(e[0, 2] - e[0, 0]) - #3 & @@@ data], {b0, 0, 3}]




            The brown and purple residual has bigger slope around the roots in the plots. Hence for Mathematica to minimize the sum of squares in the y-dimension, the mean of the 2 data points that correspond to the big slopes are cared more about than the others. It is purpose specific whether this is appropriate. If it isn't you can add the NormFunction-option to FindFit.







            share|improve this answer














            share|improve this answer



            share|improve this answer








            edited Nov 12 '18 at 8:02

























            answered Nov 12 '18 at 7:22









            Coolwater

            14.7k32553




            14.7k32553












            • Thank you dear @Coolwater. How does Mathematica recognize {n,L} for each data?
              – Hadi Sobhani
              Nov 12 '18 at 7:25










            • Also, how does it recognize that the expr is the 3rd value of every data element?
              – J42161217
              Nov 12 '18 at 7:34






            • 1




              FindFit assumes its first argument has the form {{var1, var2, ..., varN, expr}, ... , {var1, var2, ..., varN, expr}} where {var1, var2, ..., varN} is the 4th argument of FindFit
              – Coolwater
              Nov 12 '18 at 7:36










            • Ok! Given b -> {1.27225, 1.29505, 1.28573, 1.40411} having Mean=1.31428 and Medean=1.29039 do you think Mathematica did a good job? Anyway +1 from me
              – J42161217
              Nov 12 '18 at 7:52










            • @J42161217 It uses least squares, see edit
              – Coolwater
              Nov 12 '18 at 8:05


















            • Thank you dear @Coolwater. How does Mathematica recognize {n,L} for each data?
              – Hadi Sobhani
              Nov 12 '18 at 7:25










            • Also, how does it recognize that the expr is the 3rd value of every data element?
              – J42161217
              Nov 12 '18 at 7:34






            • 1




              FindFit assumes its first argument has the form {{var1, var2, ..., varN, expr}, ... , {var1, var2, ..., varN, expr}} where {var1, var2, ..., varN} is the 4th argument of FindFit
              – Coolwater
              Nov 12 '18 at 7:36










            • Ok! Given b -> {1.27225, 1.29505, 1.28573, 1.40411} having Mean=1.31428 and Medean=1.29039 do you think Mathematica did a good job? Anyway +1 from me
              – J42161217
              Nov 12 '18 at 7:52










            • @J42161217 It uses least squares, see edit
              – Coolwater
              Nov 12 '18 at 8:05
















            Thank you dear @Coolwater. How does Mathematica recognize {n,L} for each data?
            – Hadi Sobhani
            Nov 12 '18 at 7:25




            Thank you dear @Coolwater. How does Mathematica recognize {n,L} for each data?
            – Hadi Sobhani
            Nov 12 '18 at 7:25












            Also, how does it recognize that the expr is the 3rd value of every data element?
            – J42161217
            Nov 12 '18 at 7:34




            Also, how does it recognize that the expr is the 3rd value of every data element?
            – J42161217
            Nov 12 '18 at 7:34




            1




            1




            FindFit assumes its first argument has the form {{var1, var2, ..., varN, expr}, ... , {var1, var2, ..., varN, expr}} where {var1, var2, ..., varN} is the 4th argument of FindFit
            – Coolwater
            Nov 12 '18 at 7:36




            FindFit assumes its first argument has the form {{var1, var2, ..., varN, expr}, ... , {var1, var2, ..., varN, expr}} where {var1, var2, ..., varN} is the 4th argument of FindFit
            – Coolwater
            Nov 12 '18 at 7:36












            Ok! Given b -> {1.27225, 1.29505, 1.28573, 1.40411} having Mean=1.31428 and Medean=1.29039 do you think Mathematica did a good job? Anyway +1 from me
            – J42161217
            Nov 12 '18 at 7:52




            Ok! Given b -> {1.27225, 1.29505, 1.28573, 1.40411} having Mean=1.31428 and Medean=1.29039 do you think Mathematica did a good job? Anyway +1 from me
            – J42161217
            Nov 12 '18 at 7:52












            @J42161217 It uses least squares, see edit
            – Coolwater
            Nov 12 '18 at 8:05




            @J42161217 It uses least squares, see edit
            – Coolwater
            Nov 12 '18 at 8:05











            2














            You can also use NMinimize. First we need to write cost function, i.e. residual.



            data = {{0, 0, 0}, {0, 2, 1}, {0, 4, 2.247}, {0, 6, 3.627}, {0, 8, 
            
    5.031}, {1, 0, 3.346}};
            
e[n_, L_] := 2 n + 1 + Sqrt[L (L + 1) - 3/4 L^2 + 1 + b0^4]
            
cost[b0_] =Sum[(e @@data[[i, 1 ;; 2]] - (data[[i, 3]] (e[0, 2] - e[0, 0]) + 
            
      e[0, 0]))^2, {i, 6}];
            
(*or Total[(e[#1, #2] - (#3 (e[0, 2] - e[0, 0]) + e[0, 0]))^2 & @@@ data]*)
            

            
fit = NMinimize[cost[b0] , b0]



            {0.0196376, {b0 -> 1.35462}}




            Since your cost function has only one variable you can also use grid search.



            Ordering[val,1] gives position of min value.



            b0Val = Range[0, 10, 0.0001];
            
val = cost[b0Val];
            
b0Val[[Ordering[val, 1]]]



            {1.3546}




            Note that there is another min at b0=-1.3546



            b0Val = Range[-1000, 1000, 0.001];    
            
val = cost[b0Val];    
            
b0Val[[Ordering[val, 2]]]



            {-1.3546, 1.3546}




            We can plot cost function



            $text{cost}(b0)=left(-5.031 left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+25}right)^2\+left(-3.627
            left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+
            sqrt{text{b0}^4+16}right)^2\+left(2-3.346
            left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)right)^2+left(-2.247
            left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+9}right)^2$



            Plot[cost[b0], {b0, -10, 10}]


            enter image description here






            share|improve this answer




























              2














              You can also use NMinimize. First we need to write cost function, i.e. residual.



              data = {{0, 0, 0}, {0, 2, 1}, {0, 4, 2.247}, {0, 6, 3.627}, {0, 8, 
              
    5.031}, {1, 0, 3.346}};
              
e[n_, L_] := 2 n + 1 + Sqrt[L (L + 1) - 3/4 L^2 + 1 + b0^4]
              
cost[b0_] =Sum[(e @@data[[i, 1 ;; 2]] - (data[[i, 3]] (e[0, 2] - e[0, 0]) + 
              
      e[0, 0]))^2, {i, 6}];
              
(*or Total[(e[#1, #2] - (#3 (e[0, 2] - e[0, 0]) + e[0, 0]))^2 & @@@ data]*)
              

              
fit = NMinimize[cost[b0] , b0]



              {0.0196376, {b0 -> 1.35462}}




              Since your cost function has only one variable you can also use grid search.



              Ordering[val,1] gives position of min value.



              b0Val = Range[0, 10, 0.0001];
              
val = cost[b0Val];
              
b0Val[[Ordering[val, 1]]]



              {1.3546}




              Note that there is another min at b0=-1.3546



              b0Val = Range[-1000, 1000, 0.001];    
              
val = cost[b0Val];    
              
b0Val[[Ordering[val, 2]]]



              {-1.3546, 1.3546}




              We can plot cost function



              $text{cost}(b0)=left(-5.031 left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+25}right)^2\+left(-3.627
              left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+
              sqrt{text{b0}^4+16}right)^2\+left(2-3.346
              left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)right)^2+left(-2.247
              left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+9}right)^2$



              Plot[cost[b0], {b0, -10, 10}]


              enter image description here






              share|improve this answer


























                2












                2








                2






                You can also use NMinimize. First we need to write cost function, i.e. residual.



                data = {{0, 0, 0}, {0, 2, 1}, {0, 4, 2.247}, {0, 6, 3.627}, {0, 8, 
                
    5.031}, {1, 0, 3.346}};
                
e[n_, L_] := 2 n + 1 + Sqrt[L (L + 1) - 3/4 L^2 + 1 + b0^4]
                
cost[b0_] =Sum[(e @@data[[i, 1 ;; 2]] - (data[[i, 3]] (e[0, 2] - e[0, 0]) + 
                
      e[0, 0]))^2, {i, 6}];
                
(*or Total[(e[#1, #2] - (#3 (e[0, 2] - e[0, 0]) + e[0, 0]))^2 & @@@ data]*)
                

                
fit = NMinimize[cost[b0] , b0]



                {0.0196376, {b0 -> 1.35462}}




                Since your cost function has only one variable you can also use grid search.



                Ordering[val,1] gives position of min value.



                b0Val = Range[0, 10, 0.0001];
                
val = cost[b0Val];
                
b0Val[[Ordering[val, 1]]]



                {1.3546}




                Note that there is another min at b0=-1.3546



                b0Val = Range[-1000, 1000, 0.001];    
                
val = cost[b0Val];    
                
b0Val[[Ordering[val, 2]]]



                {-1.3546, 1.3546}




                We can plot cost function



                $text{cost}(b0)=left(-5.031 left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+25}right)^2\+left(-3.627
                left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+
                sqrt{text{b0}^4+16}right)^2\+left(2-3.346
                left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)right)^2+left(-2.247
                left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+9}right)^2$



                Plot[cost[b0], {b0, -10, 10}]


                enter image description here






                share|improve this answer














                You can also use NMinimize. First we need to write cost function, i.e. residual.



                data = {{0, 0, 0}, {0, 2, 1}, {0, 4, 2.247}, {0, 6, 3.627}, {0, 8, 
                
    5.031}, {1, 0, 3.346}};
                
e[n_, L_] := 2 n + 1 + Sqrt[L (L + 1) - 3/4 L^2 + 1 + b0^4]
                
cost[b0_] =Sum[(e @@data[[i, 1 ;; 2]] - (data[[i, 3]] (e[0, 2] - e[0, 0]) + 
                
      e[0, 0]))^2, {i, 6}];
                
(*or Total[(e[#1, #2] - (#3 (e[0, 2] - e[0, 0]) + e[0, 0]))^2 & @@@ data]*)
                

                
fit = NMinimize[cost[b0] , b0]



                {0.0196376, {b0 -> 1.35462}}




                Since your cost function has only one variable you can also use grid search.



                Ordering[val,1] gives position of min value.



                b0Val = Range[0, 10, 0.0001];
                
val = cost[b0Val];
                
b0Val[[Ordering[val, 1]]]



                {1.3546}




                Note that there is another min at b0=-1.3546



                b0Val = Range[-1000, 1000, 0.001];    
                
val = cost[b0Val];    
                
b0Val[[Ordering[val, 2]]]



                {-1.3546, 1.3546}




                We can plot cost function



                $text{cost}(b0)=left(-5.031 left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+25}right)^2\+left(-3.627
                left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+
                sqrt{text{b0}^4+16}right)^2\+left(2-3.346
                left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)right)^2+left(-2.247
                left(sqrt{text{b0}^4+4}-sqrt{text{b0}^4+1}right)-sqrt{text{b0}^4+1}+sqrt{text{b0}^4+9}right)^2$



                Plot[cost[b0], {b0, -10, 10}]


                enter image description here







                share|improve this answer














                share|improve this answer



                share|improve this answer








                edited Nov 12 '18 at 15:38

























                answered Nov 12 '18 at 12:05









                Okkes Dulgerci

                4,1851816




                4,1851816






























                    draft saved

                    draft discarded




















































                    Thanks for contributing an answer to Mathematica Stack Exchange!


                    • Please be sure to answer the question. Provide details and share your research!

                    But avoid …



                    • Asking for help, clarification, or responding to other answers.

                    • Making statements based on opinion; back them up with references or personal experience.


                    Use MathJax to format equations. MathJax reference.


                    To learn more, see our tips on writing great answers.





                    Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


                    Please pay close attention to the following guidance:


                    • Please be sure to answer the question. Provide details and share your research!

                    But avoid …



                    • Asking for help, clarification, or responding to other answers.

                    • Making statements based on opinion; back them up with references or personal experience.


                    To learn more, see our tips on writing great answers.




                    draft saved


                    draft discarded














                    StackExchange.ready(
                    function () {
                    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fmathematica.stackexchange.com%2fquestions%2f185829%2fproblem-in-using-findfit%23new-answer', 'question_page');
                    }
                    );

                    Post as a guest















                    Required, but never shown





















































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown

































                    Required, but never shown














                    Required, but never shown












                    Required, but never shown







                    Required, but never shown







                    -

                    Popular posts from this blog

                    Full-time equivalent

                    Image
                    Full-time equivalent From Wikipedia, the free encyclopedia Jump to navigation Jump to search Full-time equivalent ( FTE ) or whole time equivalent ( WTE ) is a unit that indicates the workload of an employed person (or student) in a way that makes workloads or class loads comparable [1] across various contexts. FTE is often used to measure a worker's or student's involvement in a project, or to track cost reductions in an organization. An FTE of 1.0 is equivalent to a full-time worker or student, while an FTE of 0.5 signals half of a full work or school load. [2] Contents 1 U.S. Federal Government 2 In education 2.1 Example 3 Notes 4 References U.S. Federal Government [ edit ] In the U.S. Federal Government, FTE is defined by the Government Accountability Office (GAO) as the number of total hours worked divided by the maximum number of compensable hours in a full-time schedule as
                    Read more

                    Bicuculline

                    Image
                    Bicuculline From Wikipedia, the free encyclopedia Jump to navigation Jump to search Bicuculline Clinical data ATC code none Identifiers IUPAC name (6 R )-​6-​[(5 S )-​6-​methyl-​5,​6,​7,​8-​tetrahydro​[1,3]dioxolo​[4,5- g ]​isoquinolin-​5-​yl]​furo[3,4- e ]​[1,3]​benzodioxol-​8(6 H )-​one CAS Number 485-49-4   Y PubChem CID 10237 IUPHAR/BPS 2312 ChemSpider 9820   Y UNII Y37615DVKC ChEBI CHEBI:3092   N ChEMBL ChEMBL417990   N ECHA InfoCard 100.006.927 Chemical and physical data Formula C 20 H 17 N O 6 Molar mass 367.352 g/mol 3D model (JSmol) Interactive image Melting point 215 °C (419 °F) SMILES O=C1O[C@H](c3c1c2OCOc2cc3)[C@@H]5c4cc6OCOc6cc4CCN5C InChI InChI=1S/C20H17NO6/c1-21-5-4-10-6-14-15(25-8-24-14)7-12(10)17(21)18-11-2-3-13-19(26-9-23-13)16(11)20(22)27-18/h2-3,6-7,17-18H,4-5,8-9H2,1H3/t17-,18+/m0/s1   Y Key:IYGYMKDQCDOMRE-ZWKOTPCHSA-N
                    Read more

                    さくらももこ

                    Image
                    さくら ももこ 本名 非公開 [1] 生誕 ( 1965-05-08 ) 1965年5月8日 日本・静岡県清水市 （現・静岡市清水区） 死没 ( 2018-08-15 ) 2018年8月15日（53歳没） 日本・東京都目黒区 国籍 日本 職業 漫画家 エッセイスト 作詞家 活動期間 1984年 - 2018年 ジャンル 少女漫画 青年漫画 児童漫画 代表作 『ちびまる子ちゃん』 『コジコジ』 受賞 1989年：第13回講談社漫画賞少女部門 （『ちびまる子ちゃん』） 1990年：第32回日本レコード大賞 （『おどるポンポコリン』作詞） 公式サイト さくらプロダクション テンプレートを表示 さくら ももこ （1965年5月8日 - 2018年8月15日 [2] ）は、日本の漫画家、エッセイスト、作詞家、脚本家。また、自身の少女時代をモデルとした代表作のコミック『ちびまる子ちゃん』の主人公の名前でもある。静岡県清水市（現・静岡市清水区）出身。血液型はA型。身長159cm。一男の母親。 代表作のコミック『ちびまる子ちゃん』の単行本の売上は累計3000万部を超える。また、エッセイストとしても独特の視点と語り口で人気が高く、初期エッセイ集三部作『もものかんづめ』『さるのこしかけ』『たいのおかしら』はいずれもミリオンセラーを記録。 目次 1 経歴 2 人物・交遊関係 3 主な作品 3.1 漫画 3.2 エッセイ 3.3 作詞 3.4 詩集 3.5 雑誌・ムック本 3.6 翻訳 3.7 絵本、ドラマ脚本他 3.8 ラジオ出演 3.9 その他 4 関係人物 4.1 アシスタント 4.2 さくらももこを演じた人物 5 脚注 5.1 注釈 5.2 出典 6 外部リンク 経歴 年譜形式の経歴は推奨されていません 。人物の伝記は流れのあるまとまった文章で記述し、年譜は補助的な使用にとどめてください。 （ 2018年8月 ） 清水市立（当時）入江小学校
                    Read more