MOM
Never
%---------FUNZIONE------------- %Copiare la funzione in un file a parte chiamato integrando.m function f_int=integrando(delta_z,k,a,z_mp,zm,x) f_int=(exp(-1i*k*sqrt(a.^2+(z_mp-x*delta_z-zm)^2))/(sqrt(a.^2+(z_mp-x*delta_z-zm)^2))) end %---------------------- clear close all clc %Definizioni parametri elettromagnetici freq=300E6; c=3E8; lambda=c/freq; mi=4*pi*1E-7; k=2*pi/lambda; l=lambda/4; a=lambda/100; epsilon=8.854E-12; omega=2*pi*freq; Vg=1; %Discretizzazione incognite delta_z=lambda/30; N=2*l/delta_z; zm=linspace(-l,l-delta_z,N) z_mp=zm+delta_z/2; w1=0.5; w2=w1; x1=(-1/sqrt(3)+1)*0.5; x2=(1/sqrt(3)+1)*0.5; %Riempo la matrice A=zeros(N); for zi=1:N for zj=1:N if zj==zi A(zi,zj)= (mi/(4*pi))*log((sqrt(1+((4*a^2)/delta_z^2))+1)/(sqrt(1+(4*a^2)/delta_z^2)-1))-((1i*k*delta_z)/(4*pi)); else A(zi,zj) =mi/(4*pi)*((w1*delta_z*integrando(delta_z,k,a,z_mp(zj),zm(zi),x1))+ w2*delta_z*(integrando(delta_z,k,a,z_mp(zj),zm(zi),x2))); end end end %Calcolo il termine noto u = zeros(1, N); u(1) = 1; u(N) = 1; c=(cos(k*z_mp))' s=(-1i*((2*pi)*freq)*epsilon*Vg)/(2*k)*sin(k*abs(zm))'; C=-(u*inv(A)*s)/(u*inv(A)*c) I=inv(A)*C*c+inv(A)*s plot(abs(I));
Raw Text
-
Untitled
14 min ago
-
~@~[[Official!]] Men's Basketball Championship Live Streams@ReddiT at Home?
14 min ago
-
Phoenix SEO
15 min ago
-
Plumber fucks someone else's wife while her husband is at work
16 min ago
-
Dr. Oz's Guava Glow CBD Gummies
38 min ago
-
Fucked and cum on the face of an adorable anime girl
52 min ago
-
Adult Telegram Channels
55 min ago
-
No Deposit Bonus
58 min ago
-
Movstream
1 hour ago
-
Makers CBD Gummies Reviews
1 hour ago