>> format long e
>> param = [1.2; 1.4; 1.8; 2.0];
>> t = [0.1; 0.3; 0.5; 0.7];
>> y = [2.1; 2.6; 2.1; 2.7];
>> [f,gr,hess,g,jac_g] = fit_damped_osc(param,t,y)

f =

    1.232691643666291e+000


gr =

   -2.690430401693836e+000
    6.529320048130283e+000
    9.815170148375403e+000
   -2.278804231368617e+000


hess =

  Columns 1 through 3 

    4.291687428671427e+000   -6.403403965969554e+000   -1.249709395757857e+001
   -6.403403965969554e+000    2.496241766728911e+001    3.943549860688516e+001
   -1.249709395757857e+001    3.943549860688516e+001    6.648905071413478e+001
    4.919860774073265e+000   -6.248546978789287e+000   -1.423024901806703e+001

  Column 4 

    4.919860774073265e+000
   -6.248546978789287e+000
   -1.423024901806703e+001
    8.539849946798238e+000


g =

    3.040984091372767e-003
   -3.165185560417383e-001
    8.097469056367368e-002
   -1.061103905963698e+000


jac_g =

  Columns 1 through 3 

    2.103040984091373e-001   -8.137660723559378e-002   -8.137660723559378e-001
    6.850444331874785e-001   -5.199135316479518e-001   -1.733045105493173e+000
    1.090487345281837e+000   -1.459778843700894e+000   -2.919557687401788e+000
    1.147227265825411e+000   -3.033209406543804e+000   -4.333156295062578e+000

  Column 4 

    1.051520492045686e+000
    1.141740721979131e+000
    1.090487345281837e+000
    8.194480470181510e-001

>>