from graftarta2 import Graf
import turtle

def equi(T1,m1,cs1,T2,m2,cs2):
    fin=Graf("scambi di calore",600,400)
    fin.coordinate(-2,-10,8,140,True)
    #turtle.color('black')
    #sca.assi()
    t=0
    k=50
    dt=0.02
    while abs(T2-T1)>0.01:
        dQ=k*(T2-T1)*dt
        T1=T1+dQ/(m1*cs1)
        T2=T2-dQ/(m2*cs2)
        fin.plot(t,T1,'red')
        fin.plot(t,T2,'green')
        t=t+dt
    return T1
        
print equi(80,200,1,40,150,1.2)