clear; clc; close all;

sigma = 5.67051e-8;
x = linspace(0,2000);
y = sigma * x.^4 /1e3;

f = figure('Position',[50,50,600,450]); 
hold on;
plot(x,y,'LineWidth',2);
xlabel('Temperatur / Kelvin')
ylabel('Wärmestrom / kW/m²')
grid on;

ax2 = axes('Position',[0.2,0.4,0.5,0.5]); hold on;
plot(x,y,'LineWidth',2);
xlabel('Temperatur / Kelvin')
ylabel('Wärmestrom / kW/m²')
xlim([0,800])
grid on;

ax3 = axes('Position',[0.3,0.58,0.25,0.3]); hold on;
plot(x,y,'LineWidth',2);
xlabel('Temperatur / Kelvin')
ylabel('Wärmestrom / kW/m²')
xlim([0,150])
grid on;

print(f,'Thoch4.png','-dpng','-r150');


T = 20;
epsilon = 1;
A = 2;
Phi = epsilon*sigma*A*(T+273.15)^4
sprintf('%.3f kcal / Tag',Phi*60*60*24*0.000239006)