% Sam Lazarus (110121560)

% Part 1

x = [.01:.01:1000];
b = [0:.01:1000];
b = x.^3 ./ (exp(x)-1);
int_ = sum(b) * .01;

k = 1.3807 .* 10.^(-23);
c = 2.9979 .* 10.^8;
hbar = 1.0546 .* 10.^(-34);

int = k.^4 ./ (4 .* pi.^2 .* c.^2 .* hbar.^3) .* int_;


% Part 2

