Analytical Model Formulation. In this project, you are expected to undertake computational modeling of a predator-prey ecosystem in a pond consisting of sharks (predators) and tuna (prey). Your task is to formulate a mathematical model and solve an optimization problem for a harvesting strategy (i.e., policy making).Let N(t) be the tuna population and P(t) be the shark population at time t. The system of equations governing their interaction is:dNdt=N[r(1-NK)-kPND]dPdt=P[s(1-hPN)]Where:r=0.6(tuna intrinsic growth rate per year) K=5000(carrying capacity of the pond for tuna) k=0.03(predation rate per interaction) D=100(half-saturation constant for predator-prey interaction) s=0.4(shark growth rate per year) h=0.02(harvesting factor of sharks per tuna unit) Consider an initial condition of N(0)=2000 and P(0)=150.(a) Find equilibrium points (N**,P**) by solving dNdt=0 and dPdt=0, analytically either using MATLAB symbolic or by hand.I need help with this I get different value everytime