T=input('Temperatur i Kelvin?');
lambda=100e-9:10e-9:15000e-9;
Me=3.741832e-16 ./lambda .^5 ./(exp(1.438786e-2 ./lambda ./T)-1);
plot(lambda*1e6,Me)
xlabel('Våglängd i µm  ')
