function plotpert(x,i,j,a,b,m)
f=[]; for t=a:.01:b
y=restrict(x+t.*omatrix(x,i,j),m,0);
f=[f y(1)];
end
t=a:.01:b; plot(t,f)