```c++
clear;
clc;x=10:0.1:63.45842269
a=[]
for x1=x
if (10<=x1<=53.73870947)
y11=-1.824e-5*x1.^3+0.05575*x1.^2-0.1125*x1-4.432;
else
y11=0;end
a=[a,y11]
end
plot(x,a)
clear;
clc;
x =10:0.1:63.45842269a =[]
for x1=x
if (10<=x1)&&(x1<=53.73870947)
y11=-1.824e-5*x1.^3+0.05575*x1.^2-0.1125*x1-4.432;
elsey11=0;
enda = [a,y11]endplot(x,a)