Question: Fitting Debye model of resistivity

hello,

 

I am trying to fit my data (in the form of two vectors "temp" and "resistance") on  resitivity as a function of the temperature "T" using a Dabye model.

The fit, intended to give as a result the two parameters "alpha" and "theta", doesn' seem to give any results.

Here my sample code:

 

>with(plots):
with(Statistics):

>rho[0]=1.33;

>T[0]=6;

 

> NonLinearFit(rho[0]+alpha*(T-T[0])*((T-T[0])/theta)^5*(int(x^5/(exp(x)-1)(1-exp(-x)), x = 0 .. theta/(T-T[0]))),temp, resistance, T,output=[parametervalues, residuals],parameterranges=[alpha=0.001..0.01, theta=10..100],initialvalues=[alpha=0.0027,theta=50]);

I am looking forward for your suggestions.

 

Thanks in advance

Please Wait...