%FY1001/TFY4145 Mekanisk fysikk høsten 2012.
%Oving 5, oppgave 3c: Plotting av funksjonen
%u(x) = 6x^2 - x^3 sammen med noen utvalgte
%verdier at totalenergien E. Vi antar at
%partikkelen slippes med null starthastighet.
clear all;
%N = antall punkter langs x-aksen
N = 100;
%tabell med N x-verdier mellom x_min og x_max
x_min = -3;
x_max = 7;
%har valgt x_min og x_max slik at 
%u_min = u(x_max) og u_max = u(x_min)
u_max = 6*x_min^2 - x_min^3;
u_min = 6*x_max^2 - x_max^3;
x = linspace(x_min,x_max,N);
%enkleste maate aa lage tabell med u-verdier;
%legg merke til at det skrives x.^2 og ikke x^2;
%punktumet foran ^ resulterer i at
%u(1) = 6x(1)^2-x(1)^3,
%u(2) = 6x(2)^2-x(2)^3 osv
u = 6*x.^2 - x.^3;
x1_min = -1.00;
x1_max = 1.25;
x1 = linspace(x1_min,x1_max,N);
for i = 1:N
    E1(i) = 6*x1_min^2 - x1_min^3;
end;
x2_min = -2.50;
x2_max = 7;
x2 = linspace(x2_min,x2_max,N);
for i = 1:N
    E2(i) = 6*x2_min^2 - x2_min^3;
end;
x3_min = 6;
x3_max = 7;
x3 = linspace(x3_min,x3_max,N);
for i = 1:N
    E3(i) = 6*x3_min^2 - x3_min^3;
end;
%Plotter u(x) sammen med E1, E2 og E3:
plot(x,u,x1,E1,x2,E2,x3,E3,'linewidth',2);
axis([x_min x_max u_min u_max]);
% Aksetekster:
xlabel('x (dim.los lengde)');
ylabel('U(x)/U_0 (dim.lost potensial)');
% Tittel:
title('Oving 5, oppgave 3c');