function y=difference_eq_1(ni,nf,a,b,y0,x)
n=ni:nf;
yy=[y0 zeros(1,length(n))];
N=length(a); M=length(b)-1;
for k=N+1:N+length(n)
yy(k)=a*yy(k-1:-1:k-N)'+b*x(k:-1:k-M);
end
y=yy(N+1:N+length(n));
위와 같은 함수에서 ni, nf, a, b, y0, x 라는 입력 인수에 해당하는 값을 넣어주어야 함수를 사용할 수 있습니다. 즉, difference_eq_1 라는 함수를 부를 때에는 예를 들면 아래와 같이 불러야 합니다.
y = difference_eq_1(0, 99, 0, 1, 0, 0:100)
함수에 관한 도움말을 통해 함수의 사용법에 대해 익히는 것이 좋을 것 같습니다.