XProfan

Titre de la fenêtre upper$("   Dichte de luftfreiem et de luftgesättigtem Wasser")
' Q: https://de.wikipedia.org/wiki/Anomalie_des_Wassers
' (CL) Copyleft 2013-10 P.Specht; aucun Gewähr!
Windowstyle 24:Fenêtre (%maxx-640)/2,(%maxy-480)/2 - 640,480:Font 2

proc DLF'Wassserdichte luftfrei, [Gramm/litre]

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

    si T!<0:imprimer " T<0-Fehler":d'autre

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

    parameters T!
    'Pour qui Berechnung qui Dichte de luftgesättigtem
    'l'eau korrigiert on den obigen Wert pour DLG = \
    return DLF(T!)-0.004612+0.000106*T!

endproc

proc Ueberschrift

    Cls
    Imprimer "\n DICHTE [g/ltr] v.";tab(20);"LUFTFREIEM  |";;tab(40-5);"LUFTGESÄTTIGTEM WASSER , Unterschied:"
    imprimer "\n chez T [°C]";tab(20);"[Gramm/litre]";tab(40);"[Gramm/litre]";tab(60);"    [%]"
    imprimer mkstr$("=",60+10)

endproc

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

Tandis que T!<99.9

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

    si %csrlin > 29:waitinput:Ueberschrift:endif

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

    endwhile

    Waitinput
    beep
    waitinput 4000
    Fin