写一段matlab程序求解常微分方程组 f'''+(f')^2-1=0 ， f(0)=0,f'(0)=0,f'(inf)=1，并分别画出f，f',f'',f'''在0到1000上的图像

以下是解常微分方程组并画图的 MATLAB 代码：

% 定义常微分方程组
fun = @(x, y) [y(2); y(3); 1 - y(2)^2 - y(3)*y(1)];

% 设置初始条件
x0 = 0;
y0 = [0; 0; 0];

% 求解常微分方程组
[x, y] = ode45(fun, [x0, 1000], y0);

% 绘制图像
subplot(2, 2, 1);
plot(x, y(:, 1));
title('f(x)');

subplot(2, 2, 2);
plot(x, y(:, 2));
title('f''(x)');

subplot(2, 2, 3);
plot(x, y(:, 3));
title('f''''(x)');

subplot(2, 2, 4);
f3 = gradient(y(:, 2), x);
plot(x, f3);
title('f''''''(x)');

该代码使用了 MATLAB 中的 `ode45` 函数来求解常微分方程组，然后使用 `subplot` 函数将四个图像绘制在一张图中。其中，最后一个图像使用了 `gradient` 函数来求解 $f''''(x)$。

注意，由于 $f'(inf)=1$，所以我们将求解区间设置为 $[0,1000]$，这个值可以根据具体情况进行调整。