You're calling your Newton Raphson function like:
[root(j), ep(j), n(j)] = NewtonRaphson(@(x) func(x, e, l, i, w(j)),@(x) dfunc(x, e, l, i, w(j)), xs, epf, e, l, i, w(j));
How does this function call the function you pass into it as the first and second inputs?
xs_new = xs - func(xs, varargin{:}) / dfunc(xs, varargin{:});
What does varargin contain in this call to NewtonRaphson? Looking at the definition:
function[root, ep, n] = NewtonRaphson(func, dfunc, xs, epf, varargin)
func is @(x) func(x, e, l, i, w(j)).
dfunc is @(x) dfunc(x, e, l, i, w(j)).
xs is xs.
epf is epf.
Everything else is packed into varargin. That means it contains {e, l, i, w(j)}.
So your line of code starting with xs_new is trying to call func with five inputs. The first is xs, the second through fifth are the elements of varargin. But the way you called NewtonRaphson, you told MATLAB that the functions whose handles you passed as the first and second input arguments only accept one input argument.
Either just pass the anonymous functions func and dfunc into NewtonRaphson (don't try to wrap them again in an anonymous function) or (the way I'd prefer) eliminate the trailing inputs in your call to NewtonRaphson.
[root(j), ep(j), n(j)] = NewtonRaphson(func, dfunc, xs, epf, e, l, i, w(j)); % Approach 1 or
[root(j), ep(j), n(j)] = NewtonRaphson(@(x) func(x, e, l, i, w(j)), ...
@(x) dfunc(x, e, l, i, w(j)), xs, epf); % Approach 2
