Modify the solution program accordingly. To this end, declare the function u_sol to be of type TrialFunction, generate the bilinear form ‘a’ and the corresponding matrix A, using the techniques in...

Modify the solution program accordingly. To this end, declare the

function u_sol to be of type TrialFunction, generate the bilinear form ‘a’ and the

corresponding matrix A, using the techniques in Remark 14.1 in Sect. 14.2. Ten

redeclare u_sol as a Function object. Tis is done outside of the while loop.

Inside the loop, construct the vector ‘b’, work the boundary condition into the

equation, pass it to solve, where the solution is inserted into the function u_sol in

vector form u_sol.vector().