function xODE
% xODE.m
% Demonstrate how to use ode45() and ode23()
%
%   Author: HUANG Huajiang
%   Copyright 2003 UNILAB Research Center, 
%   East China University of Science and Technology, Shanghai, PRC
%   $Revision: 1.0 $  $Date: 2003/02/24 $

clear all
clc

format long

y0 = 1;
[x1,y1] = ode45(@f,[0,1],y0)
[x2,y2] = ode23(@f,[0,1],y0)
plot(x1,y1,'k-',x2,y2,'b--')
xlabel('x')
ylabel('y')

% ------------------------------------------------------------------
function dydx = f(x,y)
dydx = y - 2*x/y;