# -*- coding: utf-8 -

import numpy as np
import math
import matplotlib.pyplot as plt

f = np.linspace(0, 1, 100) # x-axel (f andel solstrålning som stannar kvar)

Fs = 1370 # Solstrålning från solen

A = 0.28 # Jordens Albedo

sigma = 5.670367*10**(-8)

y = []

nominator = Fs * (1 - A)

def denominator(fVal):
    return 4.0 * sigma * (1.0 - fVal / 2.0)

def temp(nom, denom):
    return (nom/denom)**(1.0/4.0)

for i in range(0, 100):
    y.append(temp(nominator, denominator(f[i])))

print (temp(nominator, denominator(0.77)))

plt.plot(f, y)
plt.show()

7b3317bdb00ac692ec8c7608d5c338cb86c6758f
<mask token>


def denominator(fVal):
    return 4.0 * sigma * (1.0 - fVal / 2.0)


def temp(nom, denom):
    return (nom / denom) ** (1.0 / 4.0)


<mask token>

<mask token>


def denominator(fVal):
    return 4.0 * sigma * (1.0 - fVal / 2.0)


def temp(nom, denom):
    return (nom / denom) ** (1.0 / 4.0)


for i in range(0, 100):
    y.append(temp(nominator, denominator(f[i])))
print(temp(nominator, denominator(0.77)))
plt.plot(f, y)
plt.show()

<mask token>
f = np.linspace(0, 1, 100)
Fs = 1370
A = 0.28
sigma = 5.670367 * 10 ** -8
y = []
nominator = Fs * (1 - A)


def denominator(fVal):
    return 4.0 * sigma * (1.0 - fVal / 2.0)


def temp(nom, denom):
    return (nom / denom) ** (1.0 / 4.0)


for i in range(0, 100):
    y.append(temp(nominator, denominator(f[i])))
print(temp(nominator, denominator(0.77)))
plt.plot(f, y)
plt.show()

import numpy as np
import math
import matplotlib.pyplot as plt
f = np.linspace(0, 1, 100)
Fs = 1370
A = 0.28
sigma = 5.670367 * 10 ** -8
y = []
nominator = Fs * (1 - A)


def denominator(fVal):
    return 4.0 * sigma * (1.0 - fVal / 2.0)


def temp(nom, denom):
    return (nom / denom) ** (1.0 / 4.0)


for i in range(0, 100):
    y.append(temp(nominator, denominator(f[i])))
print(temp(nominator, denominator(0.77)))
plt.plot(f, y)
plt.show()

# -*- coding: utf-8 -

import numpy as np
import math
import matplotlib.pyplot as plt

f = np.linspace(0, 1, 100) # x-axel (f andel solstrålning som stannar kvar)

Fs = 1370 # Solstrålning från solen

A = 0.28 # Jordens Albedo

sigma = 5.670367*10**(-8)

y = []

nominator = Fs * (1 - A)

def denominator(fVal):
    return 4.0 * sigma * (1.0 - fVal / 2.0)

def temp(nom, denom):
    return (nom/denom)**(1.0/4.0)

for i in range(0, 100):
    y.append(temp(nominator, denominator(f[i])))

print (temp(nominator, denominator(0.77)))

plt.plot(f, y)
plt.show()
