function dy=lxn(x,y)
dy=y^2+x^2+x*y; % 保存为 lxn.m 文件后
-----------------------------------------------------------------
运行以下程序：
clear;clc
[x,y]=ode45('lxn',[0 1],0); % y(0)=0
plot(x,y,'--')
先要求解微分方程，然后得到y关于x的函数，然后在画图
