Can someone make a MATLAB Code to solve this Problem?

Orbit:

$0=30$deg

i$=20$deg

$=60$deg

a$=20,000$km

e$=0\mathrm{.}3$

$0=40$deg

t$0$based on$$tp$=0$s

Earth:

$=398,600$km$3($s$)$

R$0=6378\mathrm{.}14$km

J$2=0\mathrm{.}00108263$

Compute and Graph one complete nominal orbit using t$$as $$the independent variable. $($J$2=0)$

t$0$t$$t$0+$P

Graph Units: Km

$(10,000$points is $$sufficient$)$

Graph should show

Orbit: Solid Blue

Initial Point : Asterisk Blue

Ascending Node

Point and Radius: AsteriskDash Blue

XY $$Plane $1$it

Quadrant Arc: Solid Green

Origin and Axes: AsteriskSolid Black

Graph should use

l$=50000$;$,7$km

aw$=[-$l$+$l$-$l$+$l$-$l$+$l$]*0\mathrm{.}7$;

av$=[130,20]$; view$($av$)$

grid