I am trying to solve wave equation in FreeFem using standard finite element. Due to the source term the approximate solution become unbounded. Can anyone tell me why is this happening? I am attaching the code here. Thanks in advance.
wave 2.edp (2.5 KB)
Debug:
line 74 u1=uexact(0) → u1=uexact(dt)
line 89 u=uexact(0) → u=uexact(t)
line 84 missing dt^2
wave2.edp (2.5 KB)