function[u,I]=s174(m,N,r,a,c)
%m為第一年所所使用的設備在年初時的役齡
%N為所考慮的年數
%r為具有不同役齡的設備使用一年可得到的經濟收入向量
%a為具有不同役齡的設備使用一年所需的維修費用向量
%c為賣掉具有不同役齡的設備同時買進一臺新設備的純支出費用向量
for i=1:m+N
    V(i,N+1)=0;
end
for k=N:-1:2
    for i=1:m+k-1
       if r(i+1)-a(i+1)+V(i+1,k+1)<r(1)-a(1)-c(i+1)+V(1,k+1)
           V(i,k)=r(1)-a(1)-c(i+1)+V(1,k+1);
           w(i,k)=0;
       else
           V(i,k)=r(i+1)-a(i+1)+V(i+1,k+1);
           w(i,k)=1;
       end
    end
end
if r(m+1)-a(m+1)+V(m+1,2)<r(1)-a(1)-c(m+1)+V(1,2)
    I =r(1)-a(1)-c(m+1)+V(1,2);
    u(1)=0;
else
     I = r(m+1)-a(m+1)+V(m+1,2);
     u(1)=1;
end
x(1)=m;
for k=2:N
    x(k)=x(k-1)*u(k-1)+1;
    u(k)=w(x(k),k);
end