function solve_ode
clear all;clc
format long
global VG u
VG=3.42;
u=0.0617;
tspan=1:60;

y0=[1 1];
[t y]=ode45(@odes,tspan,y0);
figure(1),plot(t,y(:,1),'b'),legend('G')
figure(2),plot(t,y(:,2),'r'),legend('X')



function f=odes(t,y)
global VG u
G=y(1);
X=y(2);

f(1)=VG*180*X;
f(2)=u*X;
f=f';