2 %
this routine computes the stiffness matrix
for a constant c=1.
3 %K_one = int_omega (grad phi_j) . grad phi_i dx
5 c_ones=ones(1,length(c));
7 [K_one,dummy1,dummy2]=assema(p,t,c_ones,a,f);
function K_one = H10_scalar_product_matrix(p, t, c, a, f)
this routine computes the stiffness matrix for a constant c=1. K_one = int_omega (grad phi_j) ...