So I have been playing with code and I have came up with the following:
%%Given Definitions
dV=1.4*10^4; %units are in m/s
g=9.81;%units are in m/s^2
T=0.3;% units are in N
Mo=50000;% units are in kg
alpha= 10/1000; %units are in kg/W
%Isp=.1;
%%Given Equations
Eff=@(Isp) 1-(((Isp-5000)^2)/(5000^2));
Mp=@(Isp, Eff) (alpha*g*T*Isp)/(2*Eff);
dM=@(Isp) Mo*(1-exp(-dV/(Isp*g)));
%%Analytical Soultion
syms Isp
diff((alpha*g*T*Isp)/(2*(1-(((Isp-5000)^2)/(5000^2)))) + Mo*(1-exp(-dV/(Isp*g))))
Am I on the right track?