编写代码:解边值问题
,并写出结果,保留5个有效数字。[文件题]
参考答案
程序:
clear all;clc
clear all;
syms t x(t);
s=dsolve('D2x==-2/t*Dx+2/t^2*x+10*cos(log(t))/t^2','x(1)==1','x(3)==3');
t=vpa(s,5)
结果
t=
2.9513/t^2 - (3.1623*t^2*cos(log(t) + 0.32175) - 1.0487*t^3)/t^2
编写代码:解边值问题
,并写出结果,保留5个有效数字。[文件题]
参考答案
程序:
clear all;clc
clear all;
syms t x(t);
s=dsolve('D2x==-2/t*Dx+2/t^2*x+10*cos(log(t))/t^2','x(1)==1','x(3)==3');
t=vpa(s,5)
结果
t=
2.9513/t^2 - (3.1623*t^2*cos(log(t) + 0.32175) - 1.0487*t^3)/t^2