MATLAB: Phase portrait for system of ode

MATLABphase portraits

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

Best Answer

Refer the answer to similar topic here
Related Question