% ParabelPlot.m
%----------------------------------------------------------------------------
% Optegner en parabelbue. Parablen har formlen Kurve = a*u^2 + b*u + c
%----------------------------------------------------------------------------
% Johan Clausen
% Institut for Byggeri og Anlæg
% Ålborg Universitet
% Oktober 2010
%----------------------------------------------------------------------------
clear all ;

%----- Indledende indstillinger -----
a = -1.3 ; % Parameter på det kvadratiske led
b = 3.7 ; % Parameter på det linære led
c = -2.2 ; % Skæring med "Kurve"-aksen (y-aksen)

n = 3000 ; % [-] Antal plotpunkter
u_min = 0 ; % [-] Nedre plotgrænse. Ændr på denne for at "zoome" ind på et nulpunkt
u_max = 3.2 ; % Øvre plotgrænse. Ændr på denne for at "zoome" ind på et nulpunkt
%------------------------------------


%----- Beregning af værdier ------------------------------
U = [linspace(u_min,u_max,n)]' ; % [-] Uafhængig variabel
Led1 = a*U.^2 ; % Første led. Bemærk "." i ".^2". Dette bruges når hvert enkelt led i en vektor skal opløftes i en potens
Led2 = b*U ; % Andet led.
Led3 = c ; % Tredje led
Kurve = Led1 + Led2 + Led3 ; %  Dvs. Kurve = a*U.^2 + b*U + c
%---------------------------------------------------------

%----- Plotning ------------------------
figure(1) ; clf ; hold on ;
plot([u_min u_max],[0 0],'-k') ; % Nul-linje
plot(U,Kurve,'-r','linewidth',1.5) ; % Kurve
title('{\it u}-kurve','fontname','times','fontsize',13) ;
xlabel('{\it u} [-]','fontname','times','fontsize',13) ;
ylabel('{\it u}-kurve','fontname','times','fontsize',13) ;
axis tight
%---------------------------------------

%----- Slut -----------------------------------------------------------------------