phase portrait for system of ode

dx/dt=y

dy/dt=-(1/r+g+a+b*|x|)y/c-x/(lc)

for this system i am finding phase portraits ,my code is

clc
clear all

%data:
 r = 1430; 
    a = -0.0683; 
    b = 0.0676;
    c = 36*10^-9;
    g = -0.0676;
    l=27
[x,y]=meshgrid(-2:0.01:2.05, -2:0.01:2.05);
dx = y;
dy = -(1/c)*((1/r)+g+a+b*abs(x)).*y-x/(l*c)
streamslice(x, y, dx, dy,'filled');
% hold on
% figure
quiver(x, y, dx, dy);
% plot(dx,dy)
% figure

title('Phase portrait')
% axis tight

my output is coming fizzy but i want one single elliptical nature phase portrait

[darova]

What is the question?

Answer 1

[SaiDileep Kola]

Refer the answer to similar topic here