Versão 1.0

Objetivo:
Calcula a contribuição da rede para o calor específico (modelo de Debye).

Síntaxe:
c=crede(t,theta,zero).
't' é um vetor contendo as temperaturas, 'theta' é a temperatura de Debye
e 'zero' o limite inferior da integral de Debye, sendo opcional (caso
não seja atribuido valor, é considerado 1e-10).

Exemplo:
t=[1:300]';theta=100;
c=crede(t,theta,1e-5);

%%%%%%%%%%%%%%%%%%
Rotina otimizada para MATLAB 7.1

Rotina elaborada por Mario Reis (marior@if.uff.br).
Favor comunicar possíveis erros e/ou 'bugs'.

Histórico de alterações:
Versão 1.0 (abril 2010) *ORIGINAL*