Question

我的matlab程序无论t取什么值结果都是一样的，不知道是怎么回事？

clear;clc;
syms r;
v=0.25;b=2e-6;d=0.1;t=200;g=30;E=3.5e+10;c=100;B=1e-4;p=10e+6;D=6e+7;
r=0.1:0.001:0.5
X=0.5.*(r-d)./sqrt(b.*t);
i=erfc(X);
T=i.*d.*(g-c)./r+c;
u='(T-c)*r^2';
U=int('u',0.1,'r');
W=subs(U);
M=p.*d.^3./r.^3+(1-d.^3./r.^3).*D;
N=-0.5.*p.*d.^3./r.^3+D.*(1+0.5.*d.^3./r.^3);
R=-2.*B.*E.*W./((r.^3).*(1-v));
F=E.*B.*W/((r.^3).*(1-v))-(B./(1-v)).*E.*(T-c);
RR=R+M;
FF=F+N;
plot(r,R)
我检查后发现问题出在这：
u='(T-c)*r^2';
U=int('u',0.1,'r');
W=subs(U);
就是说无论t取什么值 得到的W的值（是一系列值）没有变化，请教该怎么解决这个问题？
补充，我这程序里面r，t应该是变量，为了看R随r的变化我给定一个t就应该看到一个r，R曲线。但这个程序是无论t怎么取值，都会出现一个相同的r，R曲线。

Answer 1

T=i.*d.*(g-c)./r+c;
T是r的函数，R不是t的函数，r不随T变化而变化

追问

i里面包含了t，每取一个t值就应该有一种T（r）的分布曲线。

追答

i=erfc(X);
i是x的反误差函数，由于误差函数的i值变化很小，无论t怎么取值，都会出现非常接近的r，R曲线，但仔细看该曲线是不一样的...........

追问

我取了t=20和200这两个时间的R值出来看，是一样的。取t=20和200两个时间的T值出来看，是不同的，也就是说，R值结果应该是不同的。

追答

给你看下erfc（t）的曲线图 ，-33时erfc（x）的值相等都为0.............t=20和200时，由X=0.5.*(r-d)./sqrt(b.*t);
i=erfc(X); 知两者erfc值一样，r，R曲线一样....................................