How do i get a direction field plot?

25 次查看（过去 30 天）
Tuân Nguyen
Tuân Nguyen 2022-4-20
评论： Sam Chak 2022-4-20
Hello, I need to plot a direction field for this differential equation
dP/dt=0.0015*P*(1-(P/6000))
how would I do it?
Thank you for helping.

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

采纳的回答

Sam Chak
Sam Chak 2022-4-20
编辑：Sam Chak 2022-4-20
Hi @Tuân Nguyen
Because there are two equilibrium points and , I think it should look like this.
[T, P] = meshgrid(0:6e3/14:6e3, 0:6e3/14:6e3);
S = 0.0015*P.*(1 - (P/6000));
L = sqrt(1 + S.^2);
U = 1./L;
V = S./L;
quiver(T, P, U, V, 0.5)
axis square
xlabel('t')
ylabel('P(t)')
  3 个评论
显示 1更早的评论
Tuân Nguyen
Tuân Nguyen 2022-4-20
Hi @Sam Chak
As the code that you replied me with, I wonder how can I extend the interval of P(t) to 8000?
thank you!
Sam Chak
Sam Chak 2022-4-20
Hi @Tuân Nguyen
In this line, you can change the range of P.
[T, P] = meshgrid(0:(6000-0)/14:6000, 0:(8000-0)/14:8000);
Here is the syntax:
x_start:h:x_end
where h is the spacing between the points.
(8000-0)/14
will generate 15 points between 0 to 8000 in P

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Programming 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by