# Is there the function of "DEBYE" in MATLAB ?

43 views (last 30 days)
ADNAN KIRAL on 22 Jan 2021
Answered: Jens on 5 Feb 2022
Hi
I need the "DEBYE model". how can I plot that in matlab ? any suggestion ?
ADNAN KIRAL on 23 Jan 2021
many thanks. Appriciate that.

Jens on 5 Feb 2022
Hello,
you mean the Debye-function, specific heat in 3D?
you can easily program it yourself. Either by numerical integration of the definition equatin, or a rational approximant:
R=8.314462618; % molar gas constant in J/mol/K CODATA 2018, exact.
xn=Theta./T;
I=@(x) x.^4.*exp(-x)./(exp(-x)-1).^2;
DEBYE=R*9*(1./xn).^3.*integral(I,0,xn,'reltol',1e-9,'abstol',1e-11);
or
N=[0.0279254827 -0.0358437761 1.05224663 17.2857105 64.7520511 226.68446];
D=[0.000119451046 -0.000210411972 0.00744002583 0.0346090463 0.337538084 ...
1.4272431 9.53997783 21.5840170 75.5614867]; % R. J. Goetsch+ (2011)
res1=R*polyval(N,Tn)./polyval(D,Tn);
res=reshape(res1,size(T));

### Categories

Find more on Logical in Help Center and File Exchange

R2020a

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!