The Power of the Atom

logo:=disk([6,4],1,thickness=2,color=yellow):

logo2:=plot([[5.75,3.5],[6.25,3.5],[6,4],[5.75,4.5],[5.25,4.5],[6,4],[6.3,4.5],[6.75,4.5],[6,4],[5.75,3.5]],color=black,filled=true):

mainlogo:=display([logo2,logo]):

cplant:=display(polygon([[2.5,0],[3,15],[4,15],[4.5,0]],color="firebrick",filled=true,scaling=constrained)):

cplant2:=display(rectangle([3.5,0],[14,6],color="dimgray")):

cplantwindow1:=plot([[6,3],[6,5],[8,5],[8,3],[6,3]],color="DarkTurquoise",filled=true, transparency=0):

coalplant:=display([cplantwindow1,cplant,cplant2]):

cplantwindow2:=plot([[10,3],[10,5],[12,5],[12,3],[10,3]],color="DarkTurquoise",filled=true, transparency=0):

coalplant:=display([cplantwindow1, cplantwindow2, cplant, cplant2]):

smoke1:=disk([3.5, 15], 0.5, thickness=2, color=gray):

smoke2:=disk([3.5, 15.3], 0.5, thickness=2, color=gray):

smoke3:=disk([3.25, 15], 0.5, thickness=2, color=gray):

smoke4:=disk([3.8, 15], 0.5, thickness=2, color=gray):

mainsmoke:=display(smoke1, smoke2, smoke3, smoke4):

p1:=display([coalplant, mainsmoke]):

p2:=display([mainsmoke, coalplant, translate(mainsmoke, 1, 1)]):

p3:=display([mainsmoke, coalplant, translate(mainsmoke, 2, 2)]):

p4:=display([mainsmoke, coalplant, translate(mainsmoke, 3, 3)]):

p5:=display([coalplant, textplot([6, 11, "This is a coal power plant"], 'align'={'above', 'right'}, font=[Times,14])]):

p17:=display([coalplant, textplot([6, 11, "Everything seems to be ok, right?"], 'color'="black", 'align'={'above', 'right'}, font=[Times,14])]):

p33:=display([coalplant, textplot([6, 11, "WRONG!!!"], 'color'="red",'font'=["times", "roman", 20], 'align'={'above', 'right'}, font=[Times,14])]):

p48:=display([coalplant, textplot([5, 10, "Coal power plants produce millions of tons of CO2"], 'color'="black", 'align'={'above', 'right'}, font=[Times,14])]):

badsmoke1:=display([smoke1,smoke2,smoke3,smoke4, textplot([3.8,15,"CO2"], 'color'="black", 'align'='right', font=[Times,14])]):

smoketype2:=ellipse([8, 13], 2,1, thickness=2, color =gray, filled=true):

badsmoke2:=display([smoketype2, textplot([10,13,"CO2"], 'color'="black",'align'='right', font=[Times,16])]):

coolingtower:=display(polygon([[14,0],[16.5,13],[18.5,13],[21,0],[14,0]],color="Gainsboro",filled=true)):

nplant:=display(rectangle([-2,0],[3.4,7],color="darkgray")):

nplant2:=display(sector([0.5,7],2.6,0..Pi),color=darkgray):

coolingtower2:=display(polygon([[18,0],[21,11],[24,11],[27,0],[19.5,0]],color="Gainsboro",filled=true)):

nplantmain:=display([mainlogo,cplantwindow2,cplant2,nplant,nplant2,coolingtower,coolingtower2]):

smoketype3:=ellipse([12, 10], 2,1, thickness=2, color =gray, filled=true):

badsmoke3:=display([smoketype3, textplot([14,10,"CO2"], 'color'="black", 'align'='right', font=[Times,16])]):

p49:=display([translate(badsmoke1, 1, 1), coalplant, translate(badsmoke1, 2, 2)]):

p50:=display([translate(badsmoke1, 1, 1), coalplant, translate(badsmoke1, 3, 3)]):

p51:=display([translate(badsmoke1, 1, 1), coalplant, translate(badsmoke1, 4, 4)]):

p52:=display([translate(badsmoke1, 1, 1), coalplant, badsmoke2]):

p53:=display([translate(badsmoke1, 1, 1), coalplant, badsmoke2,badsmoke3]):

thermo1:=disk([8,6],1,thickness=2,color=gray):

thermo2:=display(rectangle([7.5,6.5],[8.5,10.5],color=gray)):

mainthermo:=display(thermo1,thermo2, textplot([9,6,"(thermometer)"], 'color'="black",'font'=["times", "roman", 12], 'align'='right')):

thermomilk:=disk([8,6],0.5,thickness=2,color=red):

thermomilk2:=display(rectangle([7.75,5.75],[8.25,8.25],color=red)):

warming:=textplot([2,15,"CO2 causes global warming"], 'color'="black",'align'='right', font=[Times,16]):

p54:=display([thermomilk,mainthermo,warming]):

p56:=display([thermomilk,thermomilk2,translate(thermomilk2, 0, 0.5),mainthermo,warming]):

p57:=display([thermomilk,thermomilk2,translate(thermomilk2, 0, 1),mainthermo,warming]):

p58 := display([thermomilk, thermomilk2, translate(thermomilk2, 0, 1.5), mainthermo,warming]):

p59 := display([thermomilk, thermomilk2, translate(thermomilk2, 0, 2), mainthermo,warming]):

ice := display(polygon([[0, 0], [0, 6], [7, 6], [7, 0], [0, 0]], color="powderblue\134n", filled=true, thickness=0.1)):

icecy02:=display(polygon([[0.5, 0.5], [0.5, 5.5], [6, 5.5],[5.5,5.5],[6,4],[5,3],[5.5,2], [6, 0.5], [0.5, 0.5]], color="paleturquoise", filled=true, thickness=1)):

berg := display(polygon([[7, 0], [7, 6], [11, 6], [9, 3], [9, 2], [7, 0]], color="powderblue\134n", filled=true, outlines=none, thickness=0.5)):

icecy := display(polygon([[0.5, 0.5], [0.5, 5.5], [6, 5.5], [6, 0.5], [0.5, 0.5]], color="paleturquoise\134n", filled=true, thickness=1)):

wateryboi := display(rectangle([4, 0], [20, 3], color=deepskyblue)):

wateryboiye := display(rectangle([0, -4], [20, 1.5], color=deepskyblue)):

warming2:=textplot([2, 15, "CO2 causes icecaps to melt"], 'color'="black",'align'='right', font=[Times,16]):

p60:=display([wateryboiye, ice, berg, wateryboi,warming2]):

p61:=display([wateryboiye, icecy02, rotate(translate(berg, 1, -1), -Pi/6, [8, 4]), wateryboi,warming2]):

p62:=display([wateryboiye, icecy02, rotate(translate(berg, 2, -2), -Pi/6, [8, 4]), wateryboi,warming2]):

p63:=display([wateryboiye, icecy02, rotate(translate(berg, 2, -3), -Pi/6, [8, 4]), wateryboi,warming2]):

p64:=display([textplot([1, 9, "Its time to change to a new energy source..."], 'color'="black", 'align'='right', font=[Times,20])]):

p65:=display([textplot([3, 9, "Nuclear Power!"], 'color'="black", 'align'='right', font=[Times,24])]):

p66:=display([nplantmain, textplot([5, 13, "Nuclear power has 0 carbon emmisions"], 'color'="black", 'align'='right', font=[Times,20])]):

atom:=disk([12,10],1,thickness=2,color=red):

atomfull:=display([atom, textplot([14, 10, "atom"], 'color'="black", 'align'='right', font=[Times,24])]):

neutron:=disk([6,6],0.5,thickness=2,color=darkgray):

fullneutron:=display([neutron, textplot([7, 6, "Neutron"], 'color'="black", 'align'='right', font=[Times,24])]):

p67:=display([atomfull,fullneutron, textplot([2, 15, "Nuclear power works by splitting the atom"], 'color'="black", 'align'='right', font=[Times,14])]):

p68:=display([atomfull,translate(fullneutron,2,1)]):

p69:=display([atomfull,translate(fullneutron,4,2)]):

p70:=display([atomfull,translate(fullneutron,6,3)]):

p71:=display([textplot([2, 10, "Thats how Nuclear power works!"], 'color'="black", 'align'='right', font=[Times,16])]):

display([p1, p2, p3, p4, p5,p5,p5,p5,p5, p17,p17,p17,p17,p17, p33,p33,p33,p33,p33,p48,p48,p48,p48,p48,#p48,p48,p48,p48,p48, p49,p49,p50,p50,p51,p51,p52,p52,p52,p53,p53,p53, p54,p54,p54,p56,p56,p56,p57,p57,p57,p58,p58,p58, p59,p59,p59,p60,p60,p60,p61,p61,p61,p62,p62,p62, p63,p63,p63, p64,p64,p64,p64,p64,p65,p65,p65,p65,p65, p67,p67,p67,p67,p67,p68,p68,p68, p69,p69,p69,p70,p70,p70,p71,p71,p71], insequence, scaling=constrained, axes=none, size=[500,500]); # speed 4