updated the code

1 files changed, 28 insertions, 38 deletions

diff --git a/605/CH9/EX9.7/9_7.sce b/605/CH9/EX9.7/9_7.sce
index fc55c51bc..7b455df66 100755
--- a/605/CH9/EX9.7/9_7.sce
+++ b/605/CH9/EX9.7/9_7.sce
@@ -1,38 +1,28 @@
-

-Gr=0.01

-m=acosh(sqrt(10^(0.1*L)-1)/sqrt(10^(0.1*Gr-1)))/acosh(W/Wc)

-printf("\nm=%.0f\n",m)

-m=3

-E=log(coth(Gr/17.37))

-X=sinh(E/2/m)

-n=3

-gp=1

-

-for p=1:1:n

-    ap=sin((2*p-1)*%pi/2/m)

-    bp=X^2+sin(p*%pi/m)^2

-    printf("\nap=%.4f\nbp=%.4f\n",ap,bp)

-end

-gp=0.62425

-printf("\ng0=g4=1")

-printf("\np=1\tgp=0.62425")

-for p=2:1:n

-    gp=4*sin((2*(p-1)-1)*%pi/2/m)*sin((2*p-1)*%pi/2/m)/(X^2+sin((p-1)*%pi/m)^2)/gp

-     printf("\np=%.0f\tgp=%.5f",p,gp)

-end

-

-printf("\nL1=L3=%.4e H\nC1=%.4e F",75*0.62425/(2*%pi*10^8),0.9662/(75*2*%pi*10^8))

-

-

-

-

-

-

-

-

-

-

-

-

-

-

+clc
+clear
+L = 5;
+W  =400;
+ Wc = 100;
+Gr=0.01
+m=acosh(sqrt(10^(0.1*L)-1)/sqrt(10^(0.1*Gr-1)))/acosh(W/Wc)
+printf("\nm=%.0f\n",m)
+m=3
+E=log(coth(Gr/17.37))
+X=sinh(E/2/m)
+n=3
+gp=1
+
+for p=1:1:n
+    ap=sin((2*p-1)*%pi/2/m)
+    bp=X^2+sin(p*%pi/m)^2
+    printf("\nap=%.4f\nbp=%.4f\n",ap,bp)
+end
+gp=0.62425
+printf("\ng0=g4=1")
+printf("\np=1\tgp=0.62425")
+for p=2:1:n
+    gp=4*sin((2*(p-1)-1)*%pi/2/m)*sin((2*p-1)*%pi/2/m)/(X^2+sin((p-1)*%pi/m)^2)/gp
+     printf("\np=%.0f\tgp=%.5f",p,gp)
+end
+
+printf("\nL1=L3=%.4e H\nC1=%.4e F",75*0.62425/(2*%pi*10^8),0.9662/(75*2*%pi*10^8))
\ No newline at end of file