end T（ⅰ，j）=T1 end 其三，自定义被积函数，即 [email protected](fai)fai_s./(t-fai).^(3/2).*exp(-x^2./exp(a2.*(t-fai)））； 其四，使用integral数值积分函数求T(x,t)，即 T(x,t)=x/(8.91e-4*sqrt(pi))*integral(func,0,t)； 其五，使用surf绘制其【x,t,T(x,t)】三维曲面图，即surf(x, t, T)； 按上述思路编程，运行可以得到如下结果 。。
【matlab怎么写代码】
文章插图