Write a MATLAB code using the Step-Halving method along with Eulers method to solve the system of ODEs above with an adaptive step size and a desired approximate percent relative error of 0.001%.

The initial conditions are Ca(t=0)=100 and Cb(t=0)=1000. It must begin at t=0 and complete at t=10. It must end exactly at t=10. The output will be one figure with the two concentrations as a function of time. The y-axis should be labeled Concentration, and the x-axis should be labeled Time.

dca dt -Ca + CD dt Ca - 363