function [vec] = fiboe( n,va,vb,e )
if n==0
    vec = [vb-e];
    return
elseif n==1
    vec = [va-e];
    return
else
    vec = [fiboe(n-1,va,vb,e),fiboe(n-2,va,vb,e)];
end
    
end
