XProfan

Window Title upper$("   density of luftfreiem and of luftgesättigtem Wasser")
' Q: https://de.wikipedia.org/wiki/Anomalie_des_Wassers
' (CL) Copyleft 2013-10 P.Specht; No Gewähr!
Windowstyle 24:Window (%maxx-640)/2,(%maxy-480)/2 - 640,480:Font 2

proc DLF'Wassserdichte luftfrei, [gram/liters]

    parameters T!'[°C]
    'details see de.wikipedia.org/wiki/Virialgleichung

    if T!<0:print " T<0-Fehler":else

        return (999.83952 + 16.952577*T! - 7.9905127*10^(-3)*sqr(T!) \
        - 4.6241757*10^(-5)*T!*sqr(T!)+ 1.0584601*10^(-7)*sqr(sqr(T!)) \
        - 2.8103006*10^(-10)*T!*sqr(T!)*sqr(T!))/ (1 + 0.0168872*T!)

    endif

endproc

proc DLG'Wasserdichte, luftgesaettigt [gram/liters = kg/m³]

    parameters T!
    'for the Berechnung the density of luftgesättigtem
    'waters korrigiert one whom obigen worth to DLG = \
    return DLF(T!)-0.004612+0.000106*T!

endproc

proc Ueberschrift

    Cls
    Print "\n DICHTE [g/ltr] v.";tab(20);"LUFTFREIEM  |";;tab(40-5);"LUFTGESÄTTIGTEM WASSER , Difference:"
    print "\n with T [°C]";tab(20);"[gram/liters]";tab(40);"[gram/liters]";tab(60);"    [%]"
    print mkstr$("=",60+10)

endproc

var n&=0
Ueberschrift
var T! = -0.10

While T!<99.9

    case T!<0:goto "hop"
    print " ";T!;tab(20);DLF(T!);tab(40);DLG(T!);tab(60);if(DLF(T!)>0,100*(DLG(T!)-DLF(T!))/DLF(T!),"")
    inc n& :casenot n& mod 5:print

    if %csrlin > 29:waitinput:Ueberschrift:endif

        hop:
        T!=T!+0.1'°C

    endwhile

    Waitinput
    beep
    waitinput 4000
    End