% Example:
%
 function [x, s] = edo_s(n, c, f, alpha, beta, varargin)
 x = 1/n*(1:n-1);
 cx = feval(c, x, varargin {:});
 b = feval(f, x, varargin {:});
 b(1) = b(1) + alpha*n^2;
 b(n - 1) = b(n - 1) + beta*n^2;
 A = diag(-n^2*ones(n-2,1),1)+diag(-n^2*ones(n-2,1),-1) + diag(2*n^2*ones(n-1, 1), 0) + diag(cx, 0);
 s = A\b';
 x = [0, x, 1];
 s = [alpha, s', beta];
 return