248
21 Piecewise Interpolation Polynomials
Fig. 21.4 Example for the global Hamiltonian matrix for five finite elements. Each of the blocks
is given by the local matrix
A realization of the equations above is given by the SPECFUNPHYS-program E
= Hatomeig(n, m, rmax). The structure of this program is as follows:
function E = Hatomeig(n, m, rmax)
%% Finite element ...
% first argument n: number of nodes
% 2nd argument number of finite elements
% rmax integration limit
% example: n = 6; m = 1000; rmax = 100000;
%
E = Hatomeig(n, m, rmax);
% testing inputs
.....
%% Lagrange interpolation polynomial
[obj, C] = ippolynom(1,n);
... % some matrix gymnastic
%% auxiliary integrals