What's wrong with my code?

1 view (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
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
Ani Asoyan
Ani Asoyan on 6 Jun 2020
Thank you! it works

Sign in to comment.

Community Treasure Hunt

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

Start Hunting!