t = 0:0.1:10;
y = sin(t);
y2= 0.8*sin(t);
plot(t,y,t,y2,'r')
hold on
f=zeros(length(t),2);
i=1;
for t = 0:0.1:10;
f(i,1)=max([y(i),y2(i)]);
f(i,2)=min([y(i),y2(i)]);
i=i+1;
end
t = 0:0.1:10;
fy=cat(2,f(:,1)',flipdim(f(:,2),1)')';
fx=cat(2,t,flipdim(t',1)')';
fill(fx,fy,[0.2 0.5 0.9])
|