Plot for specific vibrational heat

Hello! I'm trying to plot the outputs of cv for the formula given in the picture to the range of temperatures between 300 to 3000. But it still gives me error due to the matrix size doesn't match.

dT=100;
T=300:dT:3000;
theta_v=3400;
R=287;
C_vv=((R*theta_v^2*(exp(theta_v/T)))/(((exp(theta_v/T))-1)^2*T^2)); 
figure(1)
plot(T, C_vv)
grid
xlabel('Temperature')
ylabel('Cv')
title('Vibrational Specific Heat vs Temperature for Nitrogent')       

Answer 1

[Star Strider]

Use element-wise operations:

C_vv=((R*theta_v^2*(exp(theta_v./T)))./(((exp(theta_v./T))-1).^2.*T.^2)); 
                               ↑     ↑               ↑       ↑  ↑  ↑