function dxdt = derivfunc(a,b,c,x)
% Function to calculate the output of a first order differential
% equation

dxdt = b-(c/a)*x^2;