cubic_spline.asv

efg code with matlab

% *************************************************************************
%                 TWO DIMENSIONAL ELEMENT FREE GALERKIN CODE
%                            Nguyen Vinh Phu
%                        LTDS, ENISE, Juillet 2006
% *************************************************************************
function [w,dwdr] = cubic_spline(r)
% Compute cubic spline function

if r <= 0.5
    w    = 2/3 - 4*r*r + 4*r*r*r ;
    dwdr = -8*r + 
end