function dy = f1(B,n,b,b1,om,c,ff,K,meta,y,ee,a,r,fspPosition,Mnti)
for i=1:n
G(i)=1-B(i)/K;
sumPrey(i)=0;
for k=1:n
sumPrey(i)=sumPrey(i)+a(i,k)*om(i)*B(k)^b;
end
end

for i=1:n
for j=1:n
F(i,j)=(a(i,j)*om(i)*B(j)^b)/(b1^b+c(i)*B(i)*b1^b+sumPrey(i));
end
end

for i=1:n
prey(i)=0;
pred(i)=0;
F2=F';
for j=1:n
prey(i)=prey(i)+(meta(i)*y(i)*B(i)*F(i,j));
pred(i)=pred(i)+(meta(j)*y(j)*B(j)*(F2(i,j)))/(ee(i,j)*ff(i,j));
end
end

dy=r'.*G'.*B-meta.*B+prey'-pred';