function inn_og_trans

%pi = 3.141592654;
% Strekk-kraft (N = kg*m/s*s):
S = 4.0;
% Massetettheter (kg/m);
mu1 = 0.010;
mu2 = 0.090;
% Bolgehastigheter (m/s):
v1 = sqrt(S/mu1);
v2 = sqrt(S/mu2);
% Vinkelfrekvens (1/s):
w = 10*pi;
% Dette tilsvarer en periode T = 2*pi/w = 0.2 s.
% Bolgetall (1/m):
k1 = w/v1;
k2 = w/v2;
% Amplituden til innkommende bolge (mm):
A = 5.0;
% Amplituden til reflektert bolge (mm):
B = A * (sqrt(mu2)-sqrt(mu1))/(sqrt(mu2)+sqrt(mu1));
% Amplituden til transmittert bolge (mm):
C = A * 2.0 * sqrt(mu1)/(sqrt(mu2)+sqrt(mu1));
% Det halve av antall punkter langs x-aksen:
N = 400;
% Antall tidssteg:
M = 400;
% Tabeller for x- og y-verdier, begge med 2*N elementer.
% Velger x-intervallet slik at vi faar minst 2 bolgelengder paa
% begge sider av x=0
kmin = min(k1,k2);
xmax = 2*2*pi/kmin;
ymax = 1.25*(abs(A) + abs(B));
x = linspace(-xmax,xmax,2*N);
y = linspace(-ymax,ymax,2*N);
y1 = linspace(-ymax,ymax,2*N);
y2 = linspace(-ymax,ymax,2*N);
y3 = linspace(-ymax,ymax,2*N);

axis([-xmax xmax -ymax ymax])
%hold on
xlabel('posisjon (m)');
ylabel('utsving (mm)');
title('Transmisjon og refleksjon i x = 0 av harmonisk b�lge p� streng');
% Velger tidspunkt...
t = i*10*2*pi/(w*M);
% For j = 1:N er x<0, og der er y = y_inn + y_refl:
for j = 1:N
    y(j) = A*sin(k1*x(j)-w*t);
    y1(j) = B*sin(k1*x(j)+w*t);
    y2(j) = A*sin(k1*x(j)-w*t) + B*sin(k1*x(j)+w*t);
end;
% For j = N+1:2*N er x>0, og der er y = y_trans:
for j = N+1:2*N
    y(j) = C*sin(k2*x(j)-w*t);
    y1(j) = C*sin(k2*x(j)-w*t);
    y2(j) = C*sin(k2*x(j)-w*t);
end;
plot(x,y2,x,y1,x,y);
end
