Thanks for the answer, also If I have a quadratic expression on the rhs will this work. I tried to write like this. will this work too?
X=zeros (1,1000);
time =0:0.05:50;
k=0.095;
for i=2:1000
X(i) = (k*(0.305*X(i-1)^2)-(0.554X(i-1))+0.249)*0.05)+X(i-1);
end
plot (time,X)
hold on
plot(time25, alpha25)
end