四阶Runge-Kutta法求解初值问题，dy/dx=x^2-y^2,y(-1)=0,x属于【-1，0】，h=0.1,精确到小数点后四位，请用MATLAB程序实现并给出程序

四阶Runge-Kutta方法是一种常用的数值积分方法，可以应用于微分方程的求解。对于给定的问题 `dy/dx = x^2 - y^2` 和初始条件 `y(-1) = 0`，我们可以在MATLAB中编写程序如下：

syms y(x) % 宣告y为符号变量
eqn = diff(y, x) == x^2 - y^2; % 微分方程

% 初始条件
y0 = sym('0'); % y(-1) = 0
x0 = -1;
xfinal = 0; % 结束区间
h = 0.1; % 步长

% 计算函数值的辅助函数
auxFcn = @(y) feval(eqn, x0, y);
% 四阶Runge-Kutta公式
rk4FourthOrder = @(yi, xi, h) (h/6)*(auxFcn(yi) + 2*auxFcn(xi + h/2) + 2*auxFcn(xi + h) + auxFcn(xi + h));

% 求解过程
[x, ySol] = ode45(@(x,y) rk4FourthOrder(y, x, h), x0:h:xfinal, y0);

% 精确到小数点后四位
ySol = double(num2str(ySol, '%.4f'));

% 显示最终解
fprintf('y(%d) = %.4f\n', int2str(find(x == xfinal)), ySol(end));

运行上述代码会得到 `y(0)` 的近似值。注意，在实际应用中，可能需要多次循环以获得更稳定的解，因为数值方法可能会受到初始猜测和步长的影响。