MATLAB: How to fix odearguments error? MATLAB

ode45

I have a homework. I write a code but it gives some error.
Error & output:
dy1 =
-0.1300
0.1300
0.0083
Output argument "y1" (and maybe others) not assigned during call to "problem8>Bqfun1".
Error in odearguments (line 90)
f0 = feval(ode,t0,y0,args{:}); % ODE15I sets args{1} to yp0.
Error in ode45 (line 115)
odearguments(FcnHandlesUsed, solver_name, ode, tspan, y0, options, varargin);
Error in problem8 (line 2)
[T,Y] = ode45(@Bqfun1,[0 720], [50 1 0]);
And this is my code:
function problem8
[T,Y] = ode45(@Bqfun1,[0 720], [50 1 0]);
plot(T,Y(:,1), 'b-', 'linewidth', 1.5)
hold on
plot(T,Y(:,2), 'g-' , 'linewidth', 1.5)
xlim([0 720])
xlabel('t', 'fontsize', 12)
ylabel('S,I,R', 'fontsize', 12)
legend('S,I,R')
function y1 = Bqfun1(t,y)
dy1 = zeros(5,1);
beta = 0.0026; %try for 0.0006 and 0.0013 too
gamma = 0.0083;
dy1 = [-beta*y(1)*y(2)
beta*y(1)*y(2)
gamma*y(2)]
Where is my mistake and how can I fix it?

Best Answer

You are missing 'd' from the function definition
function dy1 = Bqfun1(t,y)
%^ this is missing