XProfan

Título de la ventana upper$("   Dichte de luftfreiem y de luftgesättigtem Wasser")
' Q: https://de.wikipedia.org/wiki/Anomalie_des_Wassers
' (CL) Copyleft 2013-10 P.Specht; Keine Gewähr!
Windowstyle 24:Ventana (%maxx-640)/2,(%maxy-480)/2 - 640,480:Font 2

proc DLF'Wassserdichte luftfrei, [Gramm/Liter]

    parámetros T!'[°C]
    'Details siehe de.wikipedia.org/wiki/Virialgleichung

    if T!<0:imprimir " T<0-Fehler":más

        volver (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 [Gramm/Liter = kg/m³]

    parámetros T!
    'Für el Berechnung el Dichte de luftgesättigtem
    'Wasser korrigiert uno el obigen Valor después de DLG = \
    volver DLF(T!)-0.004612+0.000106*T!

ENDPROC

proc Ueberschrift

    Cls
    Imprimir "\n DICHTE [g/ltr] v.";tab(20);"LUFTFREIEM  |";;tab(40-5);"LUFTGESÄTTIGTEM WASSER , Diferencia:"
    imprimir "\n en T [°C]";tab(20);"[Gramm/Liter]";tab(40);"[Gramm/Liter]";tab(60);"    [%]"
    imprimir mkstr$("=",60+10)

ENDPROC

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

Mientras que T!<99.9

    caso T!<0:goto "hop"
    imprimir " ";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:imprimir

    if %csrlin > 29:waitinput:Ueberschrift:endif

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

    endwhile

    Waitinput
    beep
    waitinput 4000
    End