zxy2_5.m

最近一段数学建模源程序汇编,适用于建模初学者练习

clf,clear
a=0;b=4;c=0;d=4;n=15;
[X,Y]=meshgrid(linspace(a,b,n),linspace(c,d,n));
z=X.*Y;
Fx=cos(atan(X.*Y));Fy=sqrt(1-Fx.^2);
quiver(X,Y,Fx,Fy,0.5),hold on,axis([a,b,c,d])
%在每一网格点绘制相应的二维斜率矢量;0.5是控制矢量大小的控制参数,可以调整
[x,y]=ade45('zxy2_5f',[0,4],0.4);
plot(x,y,'r.-')