frcts = ...
@(a,b,c)rats(union([a:c:b],b))

frcts(-2,3,1/2)
frcts(1,2,2/3)