The main problem from what i can tell is that the number of '(' and ')' you use in declarations (and maybe even functions) are not correct.<br><br>Take for instance :<br><br>u0prime = beta*(sqrt(d**2 +(h +length1)**2) - h +length1))
<br><br>You open 3 '(' and close 4 ')' .<br><br>The problem is not the little test code at the end (as you illustrated yourself by moving it up and getting a different error).<br>The "Token Error: EOF in multi-line statement" usually means you made an error using too many or too little ()'s.
<br><br>I suggest carefully re-examining your code and check if everything is entered correctly using the right amount of ()'s ;-)<br><br>Hope this helps some.<br><br>Ciao - Geofram<br><br><div><span class="gmail_quote">On 10/7/06,
<b class="gmail_sendername">Chris Smith</b> <<a href="mailto:mc_anjo@tamu.edu">mc_anjo@tamu.edu</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
I'm writing a numerical program for an assignment at school. So the code<br>of my program isn't too long I've coded the formulas, which are rather<br>long, as funcions. However when I try to run my program I keep getting<br>
one of two errors. The first happens when I do a test run of my code<br>with the test portion of the code at the bottom. It keeps popping up an<br>error message that says that my import statement is spaced incorrectly.<br>
It's not supposed to be indented at all and I can't figure out why it's<br>popping up at all. If I try moving the test portion of the code up to<br>the top it gives me "Token Error: EOF in multi-line statement". I don't
<br>understand this one because I try to have the last line be the one with<br>the return statement of my last function and when the error happens it<br>adds a line to my code and the error pops up.<br><br>Can anyone tell me why I'm having these error or what I can do to get
<br>around them?<br><br>Chris Smith<br><br><br>#Functions for Numerical Program<br>#----------------------------------<br>### The sine and cosine integrals are taken from Abramowitz and Stegun.<br>### Only use the first 6 terms of the summation in the sine and cosine
<br>### integrals.<br><br><br>def Si(x):<br> sine_integral = x - x**3/18. + x**5/600. - x**7/35280. \<br> + x**9/3265920. + x**11/439084800.<br> return sine_integral<br><br>def Ci(x):<br> # Euler's constant
<br> Euler_const = 0.5772156649<br><br> cosine_integral = Euler_const + log(x) - x**2/4. + x**4/96. \<br> - x**6/4320. + x**8/322560. + x**10/36288000<br> return cosine_integral<br><br><br>def Mutual_impedance(length1, length2, stagger, d):
<br> """<br> Mutual impedance formulas for Parallel in Echelon Configuration<br> The formulas are taken from a paper by Howard King, "Mutual Impedance<br> of Unequal Length Antennas in Echelon"
<br><br> NOTE: all measurements should be entered in wavelengths<br> """<br><br> # stagger (this is the vertical separation between antenna centers)<br> # d (this is the horizontal separation between the antennas)
<br> # length1 and length2 (this is the half length of the antennas)<br><br> # vertical separation between center of antenna 1 and bottom of antenna 2<br> h = stagger - length2<br><br> # wave propagation constant and eta
<br> beta = 2*pi<br><br> # formulas to put into mutual impedance equation<br> u0 = beta*(sqrt(d**2 +(h -length1)**2) +(h -length1))<br> v0 = beta*(sqrt(d**2 +(h -length1)**2) -(h -length1))<br> u0prime = beta*(sqrt(d**2 +(h +length1)**2) - h +length1))
<br> v0prime = beta*(sqrt(d**2 +(h +length1)**2) +(h +length1))<br> u1 = beta*(sqrt(d**2 +(h -length1 +length2)**2) +(h -length1 +length2))<br> v1 = beta*(sqrt(d**2 +(h -length1 +length2)**2) - h -length1 +length2))
<br> u2 = beta*(sqrt(d**2 +(h +length1 +length2)**2) -(h +length1 +length2))<br> v2 = beta*(sqrt(d**2 +(h +length1 +length2)**2) +(h +length1 +length2))<br> u3 = beta*(sqrt(d**2 +(h -length1 +2*length2)**2) +(h -length1 +2*length2))
<br> v3 = beta*(sqrt(d**2 +(h -length1 +2*length2)**2) -(h -length1 +2*length2))<br> u4 = beta*(sqrt(d**2 +(h +length1 +2*length2)**2) -(h +length1 +2*length2))<br> v4 = beta*(sqrt(d**2 +(h +length1 +2*length2)**2) +(h +length1 +2*length2))
<br> w1 = beta*(sqrt(d**2 +h**2) -h)<br> y1 = beta*(sqrt(d**2 +h**2) +h)<br> w2 = beta*(sqrt(d**2 +(h +length2)**2) -(h +length2))<br> y2 = beta*(sqrt(d**2 +(h +length2)**2) +(h +length2))
<br> w3 = beta*(sqrt(d**2 +(h +2*length2)**2) -(h +2*length2))<br> y3 = beta*(sqrt(d**2 +(h +2*length2)**2) +(h +2*length2))<br><br> R12 = 15*(cos(beta*(length1 - h))*(Ci(u0) +Ci(v0) -Ci(u1) -Ci(v1)) \
<br> +sin(beta*(length1 - h))*(-Si(u0) +Si(v0) +Si(u1) -Si(v1)) \<br> +cos(beta*(length1 + h))*(Ci(u0prime) +Ci(v0prime) -Ci(u2) -Ci(v2)) \<br> +sin(beta*(length1 +h))*(-Si(u0prime) +Si(v0prime) +Si(u2) -Si(v2)) \
<br> +cos(beta*(length1 -2*length2 -h))*(-Ci(u1) -Ci(v1) +Ci(u3) +Ci(v3)) \<br> +sin(beta*(length1 -2*length2 -h))*(Si(u1) -Si(v1) -Si(u3) +Si(v3)) \<br> +cos(beta*(length1 +2*length2 +h))*(-Ci(u2) -Ci(v2) +Ci(u4) +Ci(v4)) \
<br> +sin(beta*(length1 +2*length2 +h))*(Si(u2) -Si(v2) -Si(u4) +Si(v4)) \<br> +2*cos(beta*length1)*cos(beta*h)*(-Ci(w1) -Ci(y1) +Ci(w2) +Ci(y2)) \<br> +2*cos(beta*length1)*sin(beta*h)*(Si(w1) -Si(y1) -Si(w2) +Si(y2)) \
<br> +2*cos(beta*length1)*cos(beta*(2*length2 +h))*(Ci(w2) +Ci(y2) -Ci(w3) -Ci(y3)) \<br> +2*cos(beta*length1)*sin(beta*h*(2*length2 +h))*(-Si(w2) +Si(y2) -Si(w3) +Si(y3)))<br><br> X12 = 15*(cos(beta*(length1 - h))*(-Si(u0) -Si(v0) +Si(u1) +Si(v1)) \
<br> +sin(beta*(length1 - h))*(-Ci(u0) +Ci(v0) +Ci(u1) -Ci(v1)) \<br> +cos(beta*(length1 + h))*(-Si(u0prime) -Si(v0prime) +Si(u2) +Si(v2)) \<br> +sin(beta*(length1 +h))*(-Ci(u0prime) +Ci(v0prime) +Ci(u2) -Ci(v2)) \
<br> +cos(beta*(length1 -2*length2 -h))*(Si(u1) +Si(v1) -Si(u3) -Si(v3)) \<br> +sin(beta*(length1 -2*length2 -h))*(Ci(u1) -Ci(v1) -Ci(u3) +Ci(v3)) \<br> +cos(beta*(length1 +2*length2 +h))*(Si(u2) +Si(v2) -Si(u4) -Si(v4)) \
<br> +sin(beta*(length1 +2*length2 +h))*(Ci(u2) -Ci(v2) -Ci(u4) +Ci(v4)) \<br> +2*cos(beta*length1)*cos(beta*h)*(Si(w1) +Si(y1) -Si(w2) -Si(y2)) \<br> +2*cos(beta*length1)*sin(beta*h)*(Ci(w1) -Ci(y1) -Ci(w2) +Ci(y2)) \
<br> +2*cos(beta*length1)*cos(beta*(2*length2 +h))*(-Si(w2) -Si(y2) +Si(w3) +Si(y3)) \<br> +2*cos(beta*length1)*sin(beta*h*(2*length2 +h))*(-Ci(w2) +Ci(y2) -Ci(w3) +Ci(y3)))<br><br> mut_imp = complex(R12, X12)
<br> return mut_imp<br><br>from math import *<br>length1 = 0.45<br>length2 = 0.65<br>stagger = 0.1<br>d = 0.2<br><br>impedance = Mutual_impedance(length1, length2, stagger, d)<br>print impedance<br><br><br>_______________________________________________
<br>Tutor maillist - <a href="mailto:Tutor@python.org">Tutor@python.org</a><br><a href="http://mail.python.org/mailman/listinfo/tutor">http://mail.python.org/mailman/listinfo/tutor</a><br><br><br></blockquote></div><br>