weifen912.m

来自「国外计算机科学教材系列 numerial method using matla」· M 代码 · 共 12 行

clc,clear
[t,y] = meshgrid(1:5,4:-1:1);
dt = ones(4,5);
dy = (t-y)/2;
quiver(t,y,dt,dy);
hold on;
x = 0:0.01:5;
z1 = 3 * exp(-x/2) - 2 + x;
z2 = 6 * exp(-x/2) - 2 + x;
plot(x,z1,x,z2);
hold off