// input image: map of bandgap from fit in [eV]
image bandgap := getfrontimage()

number egan, einn, bow

egan = 3.4
einn = 0.69
bow = 1.36

image content = -((einn - egan - bow) + sqrt((einn-egan-bow)**2 - 4*(egan- bandgap )*bow))/(2*bow)

content.showimage()