matlab的微积分编程问题
求解该函数在两种情形下对应点x处的二阶导数值，请给出思路与程序，谢谢
syms a x
y = a*exp(x)/sqrt(a+x^2);
y_1 = diff(y, x);
y_2 = diff(y_1, x);
a1 = 2;
x1 = 3*a1;
y2_1 = double(subs(y_2, [a, x], [a1, x1]));
a2 = 3;
x2 = 2*a2;
y2_2 = double(subs(y_2, [a, x], [a2, x2]));
disp(['a=2, x=3a, y2 = ' num2str(y2_1)]);
disp(['a=3, x=2a, y2 = ' num2str(y2_2)]);
以下内容部分参考ChatGPT模型：
思路：
根据二阶导数的定义，可以使用以下公式求解：
f''(x) = [f(x+h) - 2f(x) + f(x-h)] / h^2
其中，h为步长，一般取很小的值，如0.001。如果需要求解在x处的二阶导数，则将公式中的x替换为该值即可。
程序代码：
假设需要求解函数y = x^3在x=2处的二阶导数值，则可以编写以下代码：
% 定义函数
syms x
y = x^3;
% 求解二阶导数
h = 0.001;
y2_1 = subs(y, x, 2+h);
y2_2 = subs(y, x, 2-h);
y2 = (y2_1 - 2*y + y2_2) / h^2;
y2_value = double(y2);
其中，syms x表示定义一个符号变量x，y = x^3表示定义一个函数y = x^3。subs函数可以将x替换为指定的值，求解出对应的函数值。
最后，将求解出的二阶导数值y2_value输出即可。