x = linspace(-25, 25, 250);
part1 = @(x) (x >= 5).*(10*sqrt(x) + 25);
part2 = @(x) (x >= 0 & x <= 4).*(10*x + 25);
part3 = @(x) (x < 2) .*25;
figure
plot(x, part1(x))
hold on
plot(x, part2(x))
plot(x, part3(x))
hold off
grid
legend('part1','part2','part3', 'Location','best')
.

