lambda[1] = -1; lambda[2] = -3; lambda[3] = 14;
d[1] = 1/(exp[
lambda[1]]*(lambda[1] - lambda[2])*(lambda[1] - lambda[3])) +
1/(exp[lambda[2]]*(lambda[2] - lambda[1])*(lambda[2] -
lambda[3])) +
1/(exp[lambda[3]]*(lambda[3] - lambda[1])*(lambda[3] - lambda[2]));
d[2] = exp[
lambda[1]]*(1 -
exp[-lambda[1]/2])/((lambda[1] - lambda[2])*(lambda[1] -
lambda[3])) +
exp[lambda[
2]]*(1 -
exp[-lambda[2]/2])/((lambda[2] - lambda[1])*(lambda[2] -
lambda[3])) +
exp[lambda[
3]]*(1 -
exp[-lambda[3]/2])/((lambda[3] - lambda[1])*(lambda[3] -
lambda[2]));
d[3] = exp[lambda[1]]*lambda[1]*lambda[2] -
exp[lambda[2]]*lambda[1]*lambda[2] +
exp[lambda[2]]*lambda[3]*lambda[2] -
exp[lambda[3]]*lambda[3]*lambda[2] +
exp[lambda[3]]*lambda[1]*lambda[3] -
exp[lambda[1]]*lambda[1]*lambda[3];
Plot[(exp[lambda[2]]*lambda[2]*d[1] - exp[lambda[3]]*lambda[3]*d[1] +
lambda[3]*d[2] - lambda[2]*d[2])*
exp[lambda[1]*x]/
d[3] + (exp[lambda[3]]*lambda[3]*d[1] -
exp[lambda[1]]*lambda[1]*d[1] + lambda[1]*d[2] - lambda[3]*d[2])*
exp[lambda[2]*x]/
d[3] + (exp[lambda[1]]*lambda[1]*d[1] -
exp[lambda[2]]*lambda[2]*d[1] + lambda[2]*d[2] - lambda[1]*d[2])*
exp[lambda[3]*x]/d[3] +
exp[lambda[1]*x]*
NIntegrate[
exp[-lambda[1]*s], {s, 0, x}]/(lambda[1] - lambda[2])/(lambda[
1] - lambda[3]) +
exp[lambda[2]*x]*
NIntegrate[
exp[-lambda[1]*s], {s, 0, x}]/(lambda[1] - lambda[2])/(lambda[
3] - lambda[2]) +
exp[lambda[3]*x]*
NIntegrate[
exp[-lambda[1]*s], {s, 0, x}]/(lambda[3] - lambda[1])/(lambda[
3] - lambda[2]), {x, 0, 0.5}, PlotStyle -> {Red, Thick}]
此时NIntegrate报错,且没有显示图形,
d[1] = 1/(exp[
lambda[1]]*(lambda[1] - lambda[2])*(lambda[1] - lambda[3])) +
1/(exp[lambda[2]]*(lambda[2] - lambda[1])*(lambda[2] -
lambda[3])) +
1/(exp[lambda[3]]*(lambda[3] - lambda[1])*(lambda[3] - lambda[2]));
d[2] = exp[
lambda[1]]*(1 -
exp[-lambda[1]/2])/((lambda[1] - lambda[2])*(lambda[1] -
lambda[3])) +
exp[lambda[
2]]*(1 -
exp[-lambda[2]/2])/((lambda[2] - lambda[1])*(lambda[2] -
lambda[3])) +
exp[lambda[
3]]*(1 -
exp[-lambda[3]/2])/((lambda[3] - lambda[1])*(lambda[3] -
lambda[2]));
d[3] = exp[lambda[1]]*lambda[1]*lambda[2] -
exp[lambda[2]]*lambda[1]*lambda[2] +
exp[lambda[2]]*lambda[3]*lambda[2] -
exp[lambda[3]]*lambda[3]*lambda[2] +
exp[lambda[3]]*lambda[1]*lambda[3] -
exp[lambda[1]]*lambda[1]*lambda[3];
Plot[(exp[lambda[2]]*lambda[2]*d[1] - exp[lambda[3]]*lambda[3]*d[1] +
lambda[3]*d[2] - lambda[2]*d[2])*
exp[lambda[1]*x]/
d[3] + (exp[lambda[3]]*lambda[3]*d[1] -
exp[lambda[1]]*lambda[1]*d[1] + lambda[1]*d[2] - lambda[3]*d[2])*
exp[lambda[2]*x]/
d[3] + (exp[lambda[1]]*lambda[1]*d[1] -
exp[lambda[2]]*lambda[2]*d[1] + lambda[2]*d[2] - lambda[1]*d[2])*
exp[lambda[3]*x]/d[3] +
exp[lambda[1]*x]*
NIntegrate[
exp[-lambda[1]*s], {s, 0, x}]/(lambda[1] - lambda[2])/(lambda[
1] - lambda[3]) +
exp[lambda[2]*x]*
NIntegrate[
exp[-lambda[1]*s], {s, 0, x}]/(lambda[1] - lambda[2])/(lambda[
3] - lambda[2]) +
exp[lambda[3]*x]*
NIntegrate[
exp[-lambda[1]*s], {s, 0, x}]/(lambda[3] - lambda[1])/(lambda[
3] - lambda[2]), {x, 0, 0.5}, PlotStyle -> {Red, Thick}]
此时NIntegrate报错,且没有显示图形,