I am getting error "not enough input arguments" , I don't know how to resolve this.

Question

KAJAL KUMARI 2021-6-30

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/868633-i-am-getting-error-not-enough-input-arguments-i-don-t-know-how-to-resolve-this

评论： KAJAL KUMARI 2021-6-30

function [dn]=numparticles(t, n)
global SECT funcall No m
funcall=funcall+1;
dn=zeros(SECT,1);
sums=zeros(3,1);
    disp(['Function call ', num2str(funcall),' t value ', num2str(t) ,' of ',num2str(120),'.'])
for i=1:SECT   
   if (i-2 >= 1)
        for j=1:i-2
          sums(1)=sums(1)+beta(i-1,j)*2^(j+1-i)*n(j);
        end
   end
   if (i-1 >= 1)
         for j=1:i-1
           sums(2)=sums(2)+beta(i,j)*2^(j-i)*n(j);
         end
         dn(i)=dn(i)+ n(i-1)*sum(1)+ 0.5*beta(i-1,i-1)*n(i-1)^2-n(i)*sums(2);
   end 
   if (i+1 <=SECT)
%        	dn(i)=dn(i)+2*gamma(i+1)*n(i+1);
        for j=1:SECT
             sums(3)=sums(3)+beta(i,j)*n(j);
        end
        dn(i)=dn(i)-n(i)*sums(3);
   end   
   dn(i)=dn(i)/No(1); 
   
   sums=zeros(3,1);
end

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Jan 2021-6-30

编辑：Jan 2021-6-30

Whenever you mention an error in the forum, post the complete message. It is hard to guess, in which line your code fails, but you have this valuable information directly on the screen. So please share it with the readers.

Is beta a function or an array?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Jan 2021-6-30

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/868633-i-am-getting-error-not-enough-input-arguments-i-don-t-know-how-to-resolve-this#answer_736503

编辑：Jan 2021-6-30

在 MATLAB Online 中打开

How do you call this function? If you use the green triangle in the editor, the function is called without inputs and you will see the mentioned error.

You have to call it with inputs from the command window or from another script or function:

t = ???;
n = ???;
dn = numparticles(t, n)

3 个评论
显示 1更早的评论隐藏 1更早的评论

KAJAL KUMARI 2021-6-30

在 MATLAB Online 中打开

I am using this numparticles in another code which is dn.m

global SECT funcall m No 
relerr=1e-5;
funcall=0;
massinput = 8; %total mass of solid (kg)
d0= 0.2e-3; %particle diameter (m)
densp = 2500; %particle density (kg/m3)
Vo = pi/6*d0^3; %particle volume
massofone = densp*Vo; %mass of one particle (density*volume)
startnum= massinput/massofone; %Initial N naught
SECT = 7; % Number of sections
m=1;
No=zeros(SECT,1);
No(1)=startnum;
[t, N]=ode15s('numparticles',[0 120],No);
%%
% Plots the
figure
for ii=1:size(N)
 plot(1:SECT,N(ii,:),'r');
 xlabel('No. of Sections');
 ylabel('No. of Particles');
 pause(3/100);
end
function  [dn]=numparticles (t, n)
global SECT funcall No m
funcall=funcall+1;
dn=zeros(SECT,1);
sums=zeros(3,1);
 disp(['Function call ', num2str(funcall),' t value ', num2str(t) ,' of',num2str(120),'.'])
for i=1:SECT
     if (i-2 >= 1)
         for j=1:i-2
             sums(1)=sums(1)+beta(i-1,j)*2^(j+1-i)*n(j);
         end
     end
 if (i-1 >= 1)
      for j=1:i-1
           sums(2)=sums(2)+beta(i,j)*2^(j-i)*n(j);
      end
 dn(i)=dn(i)+ n(i-1)*sum(1)+ 0.5*beta(i-1,i-1)*n(i-1)^2-n(i)*sums(2);
 end
 if (i+1 <=SECT)
% dn(i)=dn(i)+2*gamma(i+1)*n(i+1);
      for j=1:SECT
           sums(3)=sums(3)+beta(i,j)*n(j);
      end
 dn(i)=dn(i)-n(i)*sums(3);
 end
 dn(i)=dn(i)/No(1);
 sums=zeros(3,1);
end
function B=beta(i,j)
beta0 = 1; %initial guess
 B = beta0.*1/4.*(2+((i/j)^(1/3))+((j/i)^(1/3))); %brownian motion kernel 
end							
end

Jan 2021-6-30

编辑：Jan 2021-6-30

在 MATLAB Online 中打开

If I change the line:

[t, N]=ode15s('numparticles',[0 120],No);

to

[t, N]=ode15s(@numparticles,[0 120],No);

the code runs in Matlab r2021a.

A small acceleration of the code:

function B = beta(i,j)
k = (i/j) ^ (1/3); 
B = 0.25 * (2 + k + 1/k)
end							

KAJAL KUMARI 2021-6-30

thank you so much, now it is working

请先登录，再进行评论。

I am getting error "not enough input arguments" , I don't know how to resolve this.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

回答（1 个）

3 个评论
显示 1更早的评论隐藏 1更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

I am getting error "not enough input arguments" , I don't know how to resolve this.

1 个评论 显示 -1更早的评论隐藏 -1更早的评论

回答（1 个）

3 个评论 显示 1更早的评论隐藏 1更早的评论

另请参阅

类别

标签

产品

版本

Community Treasure Hunt

WeChat

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论