What's wrong with my code?

2 views (last 30 days)
Ani Asoyan
Ani Asoyan on 6 Jun 2020
Commented: Star Strider on 6 Jun 2020
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?

Accepted Answer

Star Strider
Star Strider on 6 Jun 2020
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 Comments
Ani Asoyan
Ani Asoyan on 6 Jun 2020
Thank you! It worked
Star Strider
Star Strider on 6 Jun 2020
As always, my pleasure!

Sign in to comment.

More Answers (1)

madhan ravi
madhan ravi on 6 Jun 2020
Naming a variable pi is a bad idea
doc pi
  3 Comments
madhan ravi
madhan ravi on 6 Jun 2020
a=1:5 should be after N(a)... and implying ./
Ani Asoyan
Ani Asoyan on 6 Jun 2020
Thank you! it works

Sign in to comment.

Categories

Find more on Line Plots in Help Center and File Exchange

Community Treasure Hunt

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

Start Hunting!