function solve_odes
clear all;clc
format long
tspan=0:10:10000;
length(tspan);
a=input('ÇëÊäÈët=0ʱxµÄÊýÖµ£¬¼´x(0)=');
b=input('ÇëÊäÈët=0ʱyµÄÊýÖµ£¬¼´y(0)=');
u0=[a b];
[t u]=ode45(@odes,tspan,u0);
figure(1)
plot(t,u(:,1),'b-',t,u(:,2),'r');
ans=[600 u(600/10+1,1) u(600/10+1,2);
    1800 u(1800/10+1,1) u(1800/10+1,2);
    3600 u(3600/10+1,1) u(3600/10+1,2);
    5400 u(5400/10+1,1) u(5400/10+1,2);
    7200 u(7200/10+1,1) u(7200/10+1,2);]



function f=odes(t,u)
x=u(1);
y=u(2);
f(1)=2.44e-14/x-1.8e-15/y-1.94e-7*exp(-0.11*t);
f(2)=2.73e-15/y-1.6e-13/x;
f=f';