close all
clear all
clc

%parametri
f=400*10^6;
epsr1=1;
epsr2=9;
ur1=1;
ur2=1;
u0=4*pi*10^(-7);
eps0=8.854187871*10^(-12);
amp1= input("Inserire un ampiezza per l'onda incidente: ");
z= input("Inserire una fase per l'onda incidente: ");
p= (sqrt(epsr1) - sqrt(epsr2))/(sqrt(epsr1) + sqrt(epsr2));
tao=1+p;
w=2*pi*f;
B1=w*sqrt(u0*eps0*epsr1);
B2=w*sqrt(u0*eps0*epsr2);

%plot
t= linspace(-10,0,1000);
x=linspace(0,10,1000);
E1= amp1*cos(w*t - B1*z);
E2= amp1*p*cos(w*t + B1*z);
E3= amp1*tao*cos(w*x - B2*z);

axis([-10 10 -10 10]);
plot(t,E1,t,E2,x,E3,'linewidth',2);
legend('onda incidente','onda riflessa','onda trasmessa');
grid on