Hi Pooneh,
An example code for you to use the function would be:
sgvalues=[0; 5.953452709; 3.819761106; 0.168685691; 2.300182045; 0.015065368];
ldist=[0; 1.5; 3; 4.5; 6; 7.5];
[a,c,n,S]=variogramfit(ldist,sgvalues);
The model is set as spherical by default.
The function also accepts three more optional arguments, initial value for range, initial value for sill variance and the number of observations per lag distance. The initial values will be taken as their specified defaults if not passed as arguments. The outputs are
a : Range
c : Sill
n : Nugget (empty if nugget variance is not applied)
S : structure array with additional information
You can get the parameters of the theoretical model from the output structure S.
You can download the variogramfit.m script from the functions tab of the link you provided. It is quite well commented and you can refer to it for further information.
