I am applying FEM on Boussinesq equations. i write code for the same in matlab but while executing after few iterations it shows NaNa. Can anyone help me?
The system that i am solving is as follows
Md'' +Cd'+Kd=F+G* (1)
Where d' and d'' represents the derivatives w.r.t 't' and denotes the inner product. M, C, K are mxm matrices while F is force vector containing incident and reflected wave system and G is the vector corresponding to nonlinear terms.