t=1000:20:3000;
lambda=100e-9:10e-9:10000e-9;
[Lambda,T]=meshgrid(lambda,t);
Me=3.741832e-16 ./Lambda .^5 ./(exp(1.438786e-2 ./Lambda ./T)-1);
mesh(lambda*1e6,t,Me)
xlabel('Våglängd i µm  ')
ylabel('Temperatur i K ')
