How do I solve this differential Equation?

I have defined my variables a, b and c for a differential equation using the code below, where a is a 31x1 matrix, how do I now solve the differential equation: a*x + b*x^1/2 + c = 0 and then plot the solutions of the differential equation with respect to the thrustloading variable?

Code:

thrustloading = transpose(0:0.1:3);

delta = 1;

theta =1.0849;

CLmax=2.4;

kto=1.2;

friction=0.07;

Cdr=0.1;

Cd0=0.014;

K1=0.0375;

Mto=0.1;

tr=3;

beta=1;

gamma=1.4;

TR=1.087;

g=9.81;

sigma = delta/theta;

theta0=theta*(1+((gamma-1)/2)*Mto^2);

delta0=delta*(1+((gamma-1)/2)*Mto^2)^(gamma/(gamma-1));

density=1.225*sigma;

alpha=delta0*(1-0.49*(sqrt(Mto)));

CL=CLmax/kto^2;

Cd=(K1*CL^2)+Cd0;

XI=Cd-friction*CL;

a=-(beta/(sigma*density*g*XI))*log(1-XI./((((alpha/beta)*thrustloading)-friction)*CL))

b=tr*kto*sqrt(2*beta/(sigma*density*CLmax))

c = -3048