matlab 离散的差分方程

S(0)=300,E(0)=120,I(0)=30;b,d,e,p,v,u,c都为常数（数太繁琐，我就不都打出来了），k=0.5，用matlab求解，当t=7时，S(8),E(8),E(9)各多少S(t+1)=S(t)+(1-p)*bN(t)-dS(t)-uS(t)E(t)/(1+cS(t)^1.2)+kS(t)E(t+1)=E(t)+pbN(t)+uS(t)E(t)/(1+cS(t)^1.2)-dE(t)-eE(t)I(t+1)=I(t)+eE(t)-dI(t)-vI(t)N(t+1)=S(t)+E(t)+I(t)=(b-d)*N(t)-vI(t)+kS(t)请问，能否给出一些主要代码，或者基本的逻辑分析，用的什么方法，实在不行给出代码的大概框架就行，细节的我自己整理如果详细，还可以追加100，谢谢