How can solve this system of differential equations with two variables x, y? I will attach a photo to make it clearer.

-dx/dt = -0.6208*y^0.75*x^3; -dy/dt = -3*0.6208*y^0.75*x^3; y(0) == 0.48; x(0) == 1.44;

[Stephan]

Did you notice the edited answer?

Answer 1

[Stephan]

Hi,

this system can be integrated in a range of t=[0... ca. 0.247]:

x0 = [1.44 0.48];
tspan = [0, 0.247];
[t, y] = ode45(@odefun,tspan, x0);
plot(t, y(:,1), t, y(:,2))
function y = odefun(~, x)    
y = [-(-0.6208.*(x(2).^0.75).*(x(1).^3));...
    -(-3*0.6208.*(x(2).^0.75).*(x(1).^3))];
end

Bigger values for t2 will cause problems.

Best regards

Stephan