展开全部
^
看标题以为你要求微分方程呐,结果是32313133353236313431303231363533e58685e5aeb931333363396464画dR/dr vs. r% 画出图中的公式
% 定义微分方程函数
dRdr = @(r) 0.89 ./ r .* exp(-(log(r) + 0.84).^2 / 0.086);
% 在(0, 10]上画图
r = linspace(0.01, 3, 500);
dR = dRdr(r);
figure(10);
plot(r, dR, '^b', 'MarkerFaceColor', 'blue', 'DisplayName', 'dR/dr vs r');
xlabel('r', 'FontSize', 16);
ylabel('dR/dr', 'FontSize', 16);
legend('show');
作图的结果是
-----------------如果要求解微分方程的话-------------------------
这是一个一阶线性微分方程, 可以用龙格库塔求解. 首先用一个.m文件
定义图中的微分方程. 然后再用matlab的ode45函数求解这个方程.
举个例子, 建立一个solveFcn.m的文件如下% 调用MATLAB的`ode45`函数实现求解. 具体说明在
% MATLAB命令行中输入`doc ode45`查询.
% 简单来说, 需要两