MATLAB: Plot for specific vibrational heat

chart;plotplotting

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')

Best Answer

Use element-wise operations:
C_vv=((R*theta_v^2*(exp(theta_v./T)))./(((exp(theta_v./T))-1).^2.*T.^2));