Интегрирование по постоянной функции

Я пытаюсь интегрировать постоянную функцию в MATLAB 2017a, но я застрял. Прежде всего, когда я интегрирую с помощью следующего скрипта, я получаю правильный результат. Итак, скрипт работает для x0, который зависит от t.

function E=sol(n,k)
x0 = @(t)  t^(2);
j = 0;
E = zeros(n,1);
 while j < n+1 ;
    K = matlabFunction(subs(po(j,k))) ; 
    eval(sprintf('x%d = integral(K,0,1);',j+1)) ;
    E(j+1,1) = subs(sprintf('x%d',j+1))
    j = j+1;
 end
end

Где функция po(j,k) выглядит следующим образом:

function A_j = po(j,k)          % Adomian polynomials
 if j >0
  x =  sym('x',[1 j]);
    syms p;                            % Assinging a symbolic variable for p
    syms x0;
    S = x0+ sum(p.^(1:j) .* x) ;     % Sum of p*x up to order j
    Q =f(S,k);                        % Taking the k-th power of S, i.e. 
    A_nc = diff(Q,p,j)/factorial(j);  % Taking the j-th order derivative 
    A_j = subs(A_nc,p,0) ;            % Filling in p=0
 else
    syms x0;
    S = x0; 
    A_j =f(S,k);                      % Taking the k-th power of S,
  end
 end 

А где f(x,k),

function F = f(x,k) % Nonlinear function of k power
 F = x^k ;
end

Теперь, когда я вызываю sol(n,k), он работает. Но когда я пытаюсь изменить свою функцию x0 в sol(n,k) в постоянной функции, например,

function E=solcon(n,k)
x0 = @(t) 2.*ones(size(t));
j = 0;
E = zeros(n,1);
 while j < n+1 ;
    K = matlabFunction(subs(po(j,k))) ; 
    eval(sprintf('x%d = integral(K,0,1);',j+1)) ;
    E(j+1,1) = subs(sprintf('x%d',j+1))
    j = j+1;
 end
end

Это не работает. Как видите, я добавил *ones(size(t)); только для того, чтобы сделать его функцией t. Но, к сожалению, он все еще не работает, когда я звоню,

K = matlabFunction(subs(po(j,k))) ;

Я получил,

@()4.0

А то при звонке получаю ошибку

eval(sprintf('x%d = integral(K,0,1);',j+1)) 

Может ли кто-нибудь помочь мне интегрировать константу?

Ошибка, которую я получаю, когда звоню solcon(10,2),

Error using symengine>@()4.0
Too many input arguments.

Error in integralCalc/iterateScalarValued (line 314)
            fx = FUN(t);

Error in integralCalc/vadapt (line 132)
        [q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

Error in integralCalc (line 75)
    [q,errbnd] = vadapt(@AtoBInvTransform,interval);

Error in integral (line 88)
Q = integralCalc(fun,a,b,opstruct);

Error in solcon1 (line 7)
eval(sprintf('x%d = integral(K,0,1);',j+1)) ;

РЕДАКТИРОВАТЬ 2 Я использовал следующий скрипт,

function E=solcon(n,k)
x0 = @(t) 2.*ones(size(t));
j = 0;
E = zeros(n,1);
 while j < n+1 ;
    K = matlabFunction(subs(po(j,k))) ; 
    fstr= func2str(K)
    if fstr(3) == ')';
       x{j+1} = K*(1-0)
    else x{j+1} = integral(K,0,1)
    end
 E(j+1,1) = subs(x{j+1},1);
 j = j+1
 end
end

Но возникает следующая ошибка,

Undefined operator '*' for input arguments of type 'function_handle'.
Error in solcone1 (line 9)
x{j+1} = K*(1-0);