这样一个函数该怎么用matlab编程画图

其中1+coswt可以直接看成1，主要需要画出T和r的关系图像，K也是1，pc是常数，ae是可以改变的参量

这个这么写

kappa= 1;
rhoc = 1;
ae = 1;
alpha = 1;
P = 1;
f = @(r, t, tau) (1+0)*exp(-2*r^2./(8*kappa*(t-tau)/rhoc + ae^2))./(8*kappa*(t-tau)+ae^2*rhoc);
g = @(r, t) alpha*P/pi*integral(@(tau)f(r, t, tau), 0, t);
[r, t] = meshgrid(0:0.1:1);
T =  zeros(size(r));
for i = 1:numel(r)
   T(i) =  g(r(i), t(i));
end
surf(r, t, T)

效果

答题不易，如有帮助还望题主给个宝贵的采纳支持一下呢