Integrate a System of first order ODEs Help

Hello!

So basically I have this system of ODes

where k is a constant

which should be solved over the range [-180,0]

I am to use a function shown below that i have created in the past which outputs and V that i guess will be passed through when solving the two ODEs above.

function [V,A,Up,dV,dA] = Assignment2function(Theta,B,L,r,R,Ach,Ap,N)

Vc=269.4
V = (1+(0.5*(r-0.5))*(R + 1 - cosd(Theta) - (R^2 - (sind(Theta))^2)^0.5))*Vc
A= Ach + Ap + ((pi.*B.*L)/2).*(R+1 - cosd(Theta) + sqrt(R.^2 - (sind(Theta).^2)))
Up= pi.*L.*N.*sind(Theta) .* ( (cosd(Theta)/(R.^2 - sind(Theta).^2) + 1))
dVdtheta= (Vc.*(r-1))/2 .*(0.5.*(R.^2 - (sind(Theta)).^2).^-0.5 .* (2.*R-2.*sind(Theta).*cosd(Theta)))
dAdtheta= (pi.*B*L)/2 .* (0.5.*(R.^2 - (sind(Theta)).^2).^-0.5 .* (2.*R-2.*sind(Theta).*cosd(Theta)))

end

Having trouble coming up with the script and new function to solve the system above. Can anyone help?

thank you