def Euler(dx,dy,t0,C0, T0,xf,h):
    tn = t0
    Cn = C0
    Tn = T0
    Cn += h*dC(C0,T0)
    Tn += h*dT(C0, T0)
    tn += h
#    return Cn #usar apenas 1º iteracao
    print(tn,Cn,Tn)
#    while (tn<xf):
    for i in range(round((xf-t0)/h)):
        C0 = Cn
        T0 = Tn
        Cn += h * dC(C0, T0)
        Tn += h * dT(C0,T0)
        tn += h
#        print(tn,Cn, Tn)
    return Cn