What's wrong with my code?

here's my code..
syms a
b=2; e=0.9; q=0.4; g=0.6; k=0; d=0.5;
a=1:1:5;
N(a) = (q*b*e/a-d*e)/((q^2)*(e^2)/a+2*e*g)
pi(a)=(b-q*N*e)/a
m=N(a)
n=pi(a)
plot (a,N(a),'b')
hold on
plot (a, pi(a),'red')
When I run it, it gives me these errors.
Error using /
Matrix dimensions must agree.
Error in sev (line 31)
N(a) = (q*b*e/a-d*e)/((q^2)*(e^2)/a+2*e*g)
When I use "./", it plots only pi, but N is on x axis fixed at 0, what can I do to make it work?

 采纳的回答

Try this:
syms a
b=2; e=0.9; q=0.4; g=0.6; k=0; d=0.5;
% a=1:1:5;
N(a) = (q*b*e/a-d*e)/((q^2)*(e^2)/a+2*e*g)
fpi(a)=(b-q*N*e)/a
m=N(a)
n=fpi(a)
figure
fplot (N(a),[1 5],'b')
hold on
fplot (fpi(a), [1 5],'red')
hold off
.

5 个评论

it gives the same error :/
Error using /
Matrix dimensions must agree.
Error in sev (line 30)
N(a) = (q*b*e/a-d*e)/((q^2)*(e^2)/a+2*e*g)
Note : I commented-out this line so that it does not execute:
% a=1:1:5;
So the code I posted with that change runs without error (in R2020a) and produces:
Remove (or comment-out) the ‘a’ assighment as I did, and my code works.
Thank you! It worked
As always, my pleasure!

请先登录,再进行评论。

更多回答(1 个)

Naming a variable pi is a bad idea
doc pi

3 个评论

ok but the problem remains
a=1:5 should be after N(a)... and implying ./
Thank you! it works

请先登录,再进行评论。

类别

帮助中心File Exchange 中查找有关 Axis Labels 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by