XProfan

Deutsch

p.specht

Eine Formel für den durchschnittlichen Blutzucker-Verlauf nach Kohlehydrat-Zufuhr war im gesamten Internet nicht zu finden, bloß Bilder und Tabellen. Die Fachbücher strotzen dagegen vor Doppelintegralen und komplexen, hochtheoretischen physiologischen Modellen, mit denen wir Normalos genau nix anfangen können.

Aber da gab es doch mal ein Buch mit Funktionsgraphen und zugehörigen Formeln: Abramowitz-Stegun: Handbook of Mathematical Functions, und so findet man schließlich (Tataaaa!) den Integralcosinus in seinem ersten positiven Abschnitt. Der modelliert dann - laienhaft ausgedrückt - sogar die anfängliche Insulinverzögerung. Nachstehend die nackte, unskalierte Mathe-Funktion, wie immer ohne jede Gewähr (Die Reihenentwicklung sollte nicht überbeansprucht werden: Bei x=24 ist Schluss)!

Kompilieren Markieren Separieren

WindowTitle "Integralcosinus Ci() aus Reihenentwicklung"
Windowstyle 24:Window 0,0-%maxx,%maxy
var xh&=width(%hwnd)\2:var yh&=height(%hwnd)\2
Declare x!,y!,f!:f!=1/1000:set("decimals",17)
Main:
cls:line 0,yh& - 2*xh&,yh&: line 300,0 - 300,2*yh&:usepen 0,1,255

whileloop 0,24300,50:x!=&loop*f!

    y!=Ci(x!)
    print format$("##0.#####",x!),tab(12);y!
    lineto 300+x!*20,yh&-y!*500
    case %csrlin>50:locate 1,1

endwhile

waitinput:cls:goto "Main"

Proc Ci :parameters x!:case x!<=0:return -999999999

    case x!>24.3:print " **Err: x > 24.3** ";
    var EulerMascheroni!=0.57721566490153286:var lm2!=0
    var sum!=EulerMascheroni!+ln(x!)

    whileloop 30:lm2!=&Loop*2

        sum!=sum! + if(&loop mod 2,-1,1)*x!^Lm2!/(faku(Lm2!)*Lm2!)

    endwhile

    return sum!

endproc

proc faku :parameters x!:x!=abs(x!)

    if x!>169:print " ***Error in Faku: Overflow!***"

        waitinput:End:endif :var prod!=1

        whileloop int(x!):prod!=prod!*&loop

            endwhile:return prod!

        endproc