“…The solution domain a ≤ x ≤ b is uniformly partitioned by knots x j into M subintervals [ x j , x j + 1 ] of equal length h , j = 0,1,2, … , M − 1, where a = x 0 < x 1 < … < x n − 1 < x M = b . Our scheme for solving () requires approximate solution U ( x , t ) to the exact solution u ( x , t ) in the following form: 25‐29 where d i ( t ) are unknown functions to be determined and 20,30,31 are twice differentiable cubic trigonometric basis functions given by …”