Plotting a function with 3 differing initial values

Question

Nicholas Davis 2021-1-29

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/730753-plotting-a-function-with-3-differing-initial-values

回答： Shubham Khatri 2021-2-3

Hi all. I am trying to plot a function with 3 different scale values, so I'd like to produce a line for each scale value (included in code as a vector, notated as lambda on plot legend), though I am unsure of how to go about it. Thanks!

% Initial Values
m = linspace(0,1,10^6);
scale = [ 0.1 0.04 0.02 ];
% The Workzone
[K,E] = ellipke(m); % Elliptical Integral of First and Second Kind
% Axes for Plotting
xaxis = -1*log10(1-m);
yrecip = sqrt(8*K.*scale^2.*(K-E));
yaxis = 1./yrecip;
% Plotting
clf;   figure(1);
plot(xaxis,yaxis,'-k');
% Axes
colororder({'k','k'});
xlim([-.25,15.5]);   ylim([0,8]);
xlabel('-log10(1-m)');   ylabel('(8K(m)(\lambda)^2(K(m)-E(m))^-1/2');
yyaxis right 
ylim([0,8]);
ylabel('j');
grid on
legend('\lambda = 1/25')
% title('Graphical Solution of Eq.15');

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Shubham Khatri 2021-2-3

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/730753-plotting-a-function-with-3-differing-initial-values#answer_613427

在 MATLAB Online 中打开

Hello,

To my understanding of the question and code, I find that the dimension of scale is not correct thats why it is creating error. Also, while defining the variable yrecip the 'scale' should be dot multiplied/squared. I have pasted the corrected code below.

% Initial Values
m = linspace(0,1,10^6);
scale = 0.0000001:0.000001:1;
% The Workzone
[K,E] = ellipke(m); % Elliptical Integral of First and Second Kind
% Axes for Plotting
xaxis = -1*log10(1-m);
yrecip = sqrt(8*K.*(scale.^2).*(K-E));
yaxis = 1./yrecip;
% Plotting
clf;   figure(1);
plot(xaxis,yaxis,'-k');
% Axes
colororder({'k','k'});
xlim([-.25,15.5]);   ylim([0,8]);
xlabel('-log10(1-m)');   ylabel('(8K(m)(\lambda)^2(K(m)-E(m))^-1/2');
yyaxis right 
ylim([0,8]);
ylabel('j');
grid on
legend('\lambda = 1/25')
% title('Graphical Solution of Eq.15');