function u = g2(tkk,uk,ukk,Dt)

u = (uk + Dt*sin(ukk))/(1.0+4*Dt/(1+tkk^2));