%Программа для построения огибающей семейства нормалей к параболе
t=-4:0.01:5;
yp=t.^2-2*t+1;
hold on;
plot(t,yp,'r-'); % График y=t^2
dt=0.1;
% Построение нормалей
xn=-10:0.01:10;
for T=-1.5:dt:3.5
    yn=T^2-2*T+1-(xn-T)/(2*(T-1));
    plot(xn,yn); % Нормали
    %pause(0.2*(1-i/(n+1)));
end

xo=t-6*(t-1).^3-(t-1)+2*(t-1).*(t.^2-2*t+1);
yo=3*(t-1).^2+0.5;
%plot(xo,yo,'k');  % График y = ax(1-x)
% a2=0.9; % Значение коэффициента a
% y2=a2*x.*(1-x);
% plot(x,y2,'k');  % График y = ax(1-x)
axis([-4 4 -1. 7])

