From: Maxime Debon <maxime.debon@et...>  20080227 19:48:16

Hi, Using spatial coordinates (q_point) in matrix Ke almost as it's presented in the example 14 for matrix Fe, I may misunderstand something important. Here is the code. //  // const std::vector<Point>& qrule_points = fe>get_xyz(); for (unsigned int qp=0; qp<qrule.n_points(); qp++) { Real S = qrule_points[qp](0); // Real S = 10.0; for (unsigned int i=0; i<phi.size(); i++) { for (unsigned int j=0; j<phi.size(); j++) { Ke(i,j) += JxW[qp]*(dphi[i][qp]*dphi[j][qp]) * S; } } } //  // When the Real S (line 4) is a constant, the equation solution is the same as Matlab's results whereas it's different when S depends on the first spatial factor (line 3). What's the good way ? Best regards, Maxime 