スプライン関数の作り方

上図のようなものを作ったのですが、曲線が丸印の線を超えないようにしたいのですがどのようなコードを書き足せばいいかわかりません．

わかる方いればお願いいたします．

コードは以下のものを作りました．

% オリジナルデータの作成

x = 0:119;

y = [zeros(1,45) 90*ones(1,30) zeros(1,45)];

% 多項式近似

p12 = polyfit(x,y,12);

yp12 = polyval(p12,x);

% グラフの作成

figure

plot(x,y,'o-')

xlim([0 500])

ylim([0 500])

hold on

plot(x,yp12)

hold off

legend('元データ','2','4','12')