disp('Solution Domain Range Input')
prompt='lims='; lims=input(prompt);
disp(' ')

lims=[-3 3 ] works fine, x axis range is from -3 to 3

lims=[-3 3 -3 3 ] gives error message

odes = diff(Y)== A*Y + B
[xSol(t), ySol(t)] = dsolve(odes);
xSolvpa = vpa(xSol)
ySolvpa = vpa(ySol)
fplot(xSolvpa, lims)
hold on
fplot(ySolvpa,lims)

Error using fplot (line 145)
Invalid parameter '-3  3 -...'.

lims=[-1.1 1.5 -1 3]
fplot(xSolvpa)
axis(lims)