function ukk = funz(tk,uk)

ukk = -4*uk/(1+tk) + sin(uk);