c(t)的取值是0.05到max(t)的固定值，max(t)=y(t)+(1-d)k(t)，y(t)=3k(t)，k(t+1)=max(t)-c(t)，设定t+1期的V(t+1)=0，计算出V(t)，将V(t)赋值给V(t+1)，反复迭代，直到V(t)收敛，取c(t)的步长是0.01，得出V(0)的最大值的matlab代码

以下是使用Matlab实现的代码：

% 初始化参数
k0 = 4; % 初始资本
d = 0.1; % 折旧率
b = 0.9; % 折现率
tol = 0.0001; % 收敛阈值
max_iter = 1000; % 最大迭代次数
N = 100; % 状态空间离散化的取值个数
c = linspace(0.05, y(k0)+(1-d)*k0, N); % 离散化c(t)的取值范围
V = zeros(1, N); % 初始化V(t)为0

% 迭代计算V(t)
for iter = 1:max_iter
    % 计算当前状态下的收益
    for i = 1:N
        k = y(k0) + (1-d)*k0 - c(i); % 根据k(t+1)=max(t)-c(t)计算k(t+1)
        if k <= 0 % 如果k(t+1)小于等于0，则收益为负无穷
            G(i) = -inf;
        else % 否则，根据V(t+1)=0计算收益
            G(i) = log(c(i)) + b*V(find(c >= max(0.05, y(k)+(1-d)*k), 1));
        end
    end
    
    % 更新V(t+1)为当前的V(t)
    V_new = G;
    
    % 判断是否收敛
    if norm(V_new - V) < tol
        break;
    end
    
    % 更新V(t)
    V = V_new;
end

% 输出V(0)的最大值
fprintf('V(0)的最大值为%.4f\n', max(V));

在上述代码中，我们首先初始化了一些参数，包括初始资本$k_0$、折旧率$d$、折现率$b$、收敛阈值$tol$、最大迭代次数$max\_iter$、状态空间离散化的取值个数$N$，以及c(t)的取值范围。然后，我们使用一个循环来迭代计算V(t)。在每一次迭代中，我们首先计算当前状态下的收益，然后根据V(t+1)=0更新V(t+1)为当前的V(t)，最后判断是否收敛。如果V(t)已经收敛，则输出V(0)的最大值。

需要注意的是，在计算当前状态下的收益时，我们使用了一个if语句来判断k(t+1)是否小于等于0，如果是，则收益为负无穷。这是因为如果k(t+1)小于等于0，则意味着无论我们在该状态下选择多少消费，我们都无法获得正收益，因此该状态的价值为负无穷。