Als je geen stapgrootte dt aangeeft, dan zal er niets werken:

t := t + dt

Daarvoor is een dt waarde (1 seconde, 30 seconde, wat ook) wel nodig.
In die k-waarde zit van alles verborgen zoals massa van het water en de soortelijke warmte?

Ik heb nu een dt toegevoegd van 0,1 maar er komt nogsteeds geen grafiek uit het model

Het is me nu gelukt om een grafiek te maken, maar het verband klopt helemaal niet. het is nu een recht evenredig verband door de oorsprong, maar dan koelt het water helemaal niet af en stijgt het juist. Dit is niet de bedoeling...

De k-waarde is een evenredigheidsconstante in de volgende differentiaalvergelijking: 
Deze heb ik samen met een vriend van mij berekend door een aantal randwaarden in te voeren. Toen kwamen we uit op een k van 0,000402567

Ohhh... de geleidingsformule, in plat Nederlands P =ΔQ/Δt = λ A/d  ΔT (energie/seconde geleid door een materiaal met oppervlak A en dikte d en temperatuursverschil tussen uiteinden ΔT en als geleidingscoefficient materiaalspecifiek λ)

Die energie wordt verloren door het warmere water: P = ΔQ/Δt = m.c_w.ΔT  (energieverlies (per seconde) van m kg water bij daling in temperatuur van ΔT graden) is. 

Let op dat de beide ΔT's niet hetzelfde zijn.
Enerzijds is het de afkoeling van het water, anderzijds het verschil in temperatuur tussen de omgeving (van de kop) en het water.
Daarnaast heb je met een kopje natuurlijk te maken met wanden maar ook een open wateroppervlak aan de bovenkant, direct aan de lucht. De geleiding naar lucht is meer verdamping.

niet raar:

dT := -k * (T-Tomg) * dt

dT wordt dus onvermijdelijk een negatief getal (want k definieer je positief, en T-Tomg zal ook altijd positief zijn)

T := T - dT

een negatief getal ergens aftrekken levert een toenemend getal

groet, Jan

ik denk dat het een rechtevenredig verband LIJKT, maar als je je model lang genoeg laat lopen móet hij wegkrommen vanwege je steeds grotere (in dit foute model) temperatuurverschil. 

Allemaal heel erg bedankt! het is gelukt

Je kunt ook rekening houden met stralingseffecten: naast warmtegeleiding is er ook sprake van warmtestraling. Je hebt dan de wet van Stefan-Boltzmann nodig. Mbv Coach krijg je dan ook een mooie afkoelingscurve.