clear all; close all; clc;

h = 0.25;

tlist = 0 : 0.3: 2*pi;

Clist = -5: 0.1 : 5 ;

h1 = figure;

hold on

for s1 = 1 : length(tlist)

    for s2 = 1 : length(Clist)

        t = tlist(s1);

        C = Clist(s2);

        if abs(C + sin(t)) < 2

            quiver(t, C+ sin(t),  h, h * cos(t), 'b')

        end

    end

end

plot(tlist, sin(tlist),'r')

转载本文请联系原作者获取授权，同时请注明本文来自张江敏科学网博客。

链接地址：https://wap.sciencenet.cn/blog-100379-1209407.html?mobile=1

收藏

分享