PHY 531 Homework 3 Problem 5, Part C

clc;
x = linspace(-5,-1,100);
rho = 1e9;
T = 1e8;
k = 1.38e-23;
h = 6.626e-34;
c = 3e8;
a=(-(1./1.32e56));
b=(2*T.^5)*(k.^6);
d=(h.^5)*(c.^2);
p=(rho.^2)*(x.^7);
q=d*p;
z=exp(-2*x);
r=z.*q;
func = a*b./r;
plot(x,func)
hold on;
grid on;
xlabel("x = $-\frac{h \nu}{\kappa T}$","interpreter", "latex")
ylabel('$\frac{\partial B_{\nu}}{\partial T} (\rho \kappa_{\nu})$',"interpreter", "latex")
title('Problem 5,Part c')