function feul = feul(t,x)
    feul = .061*(70.+t/5 - x);