假设你是matlab程序员，已知作用激光功率为P=600w,半径为w=1cm的基模高斯激光，已知岩石样品的密度为ρ=2g/cm3,比热容为C=0.75J/(g.K),热传导系数为K=4.4W/(m.K),假设岩石对光吸收率为η=0.6，岩石为长10cm，宽10cm，高15cm的长方体体,初始温度T0＝300K，边界条件为第二类边界条件，根据matlab软件利用有限差分法计算激光以v=0.13cm/s的速度，照射3s后的岩石表面沿移动方向温度场和应力场

这是一个比较复杂的问题，需要对热传导方程进行数值求解。以下是一个可能的解决方案：

1. 建立模型和边界条件

首先，我们需要建立一个模型来描述热传导方程。根据热传导方程，我们可以得到如下的偏微分方程：

ρC∂T/∂t - K(∂^2T/∂x^2 + ∂^2T/∂y^2 + ∂^2T/∂z^2) = Pηδ(x-vt)δ(y)δ(z)

其中，ρ是岩石样品的密度，C是比热容，K是热传导系数，T是温度场，P是激光功率，η是岩石对光的吸收率，v是激光照射的速度，t是时间，δ是狄拉克函数。

边界条件为第二类边界条件，即表面温度梯度为0。这可以表示为：

∂T/∂n = 0

其中，n是法向量。

2. 离散化

为了将偏微分方程转换为差分方程，我们需要对空间和时间进行离散化。我们将区域离散化为nx×ny×nz个网格，时间离散化为nt个时间步长。假设空间步长为Δx，Δy，Δz，时间步长为Δt。

我们可以使用中心差分法来近似偏微分方程中的二阶导数，得到如下的差分方程：

ρC(T(i,j,k,n+1) - T(i,j,k,n))/Δt - K((T(i+1,j,k,n) - 2T(i,j,k,n) + T(i-1,j,k,n))/Δx^2 + (T(i,j+1,k,n) - 2T(i,j,k,n) + T(i,j-1,k,n))/Δy^2 + (T(i,j,k+1,n) - 2T(i,j,k,n) + T(i,j,k-1,n))/Δz^2) = Pηδ(i,j,k,n)δ(y)δ(z)

其中，i，j，k是网格索引，n是时间步长。

3. 边界处理

对于边界，我们需要使用特殊的处理方法。由于表面温度梯度为0，我们可以使用向后差分法来近似表面的导数，得到：

(T(i,j,k,n) - T(i,j,k,n-1))/Δt = 0

这意味着我们可以将表面温度设置为前一个时间步长的温度值。

4. 初始条件

我们需要设置初始条件，即将所有网格的温度设置为初始温度T0。

5. 数值求解

使用matlab软件，我们可以编写一个程序来求解差分方程。在每个时间步长内，我们可以使用矩阵运算来更新温度场。最终，我们可以得到岩石表面沿移动方向的温度场和应力场。

需要注意的是，这只是一个简单的解决方案，实际上可能会有更多的细节和问题需要考虑。