1 files changed, 38 insertions, 0 deletions
diff --git a/605/CH5/EX5.6/5_6.sce b/605/CH5/EX5.6/5_6.sce
new file mode 100755
index 000000000..35e100830
--- /dev/null
+++ b/605/CH5/EX5.6/5_6.sce
@@ -0,0 +1,38 @@
+
+Zo=50
+Er=2.08
+B=60*%pi^2/Zo/sqrt(Er)
+c=3*10^8
+f=2*10^9
+Ee=1.7875
+t=0.159*10^-6
+a=5.813*10^7
+fGHz=2
+A=Zo/60*((Er+1)/2)^(1/2)+(Er-1)/(Er+1)*(0.23+0.11/Er)
+printf("\nA=%.4f\nB=%.4f",A,B)
+h=1.59*10^-3
+if A>1.52 then
+    W=h*(8*exp(A)/(exp(2*A)-2))
+else
+    W=h*(2/%pi*(B-1-log(2*B-1)+(Er-1)/2/Er*(log(B-1)+0.39-0.61/Er)))
+end
+printf("\nW=%.6f cm",W)
+
+Ere=0.5*(2.08+1+(2.08-1)/sqrt(1+12/3.192094))-(2.08-1)*0.0001/4.6/sqrt(3.1921)
+printf("\nEre=%.4f",Ere)
+l=c/(2*f*sqrt(Ee))
+printf("\nl=%.6f m",l)
+We=W/h+0.3979*t/h*(1+log(2*h/t))
+printf("\nWe/h=%.4f",We)
+F=1+1/We*(1-1.25*t/(%pi*h)+1.25/%pi*log(2*h/t))
+ac=44.1255*10^-5*F*Zo*Ere/h*sqrt(fGHz/a)*(We+0.667*We/(We+1.444))
+printf("\nac=%.4f Np/m",ac)
+del=0.00040
+ad=10.4766*Er/(Er-1)*(Ere-1)/sqrt(Ere)*fGHz*tan(del)
+printf("\nad=%.4f Np/m",ad)
+B=2*%pi/(2*l)
+printf("\nB=%.4f rad/m",B)
+Q=B/(2*(ac+ad))
+printf("\nQ=%.1f",Q)
+
+