diff --git a/src/sphere/edfb.cpp b/src/sphere/edfb.cpp
index 04ded87d1a91de67f47373058ed33598388b4569..244c11b2a4432052e0f5c220cf4ef4b871b48b75 100644
--- a/src/sphere/edfb.cpp
+++ b/src/sphere/edfb.cpp
@@ -345,7 +345,7 @@ int main(int argc, char **argv) {
       }
       vss = xi_vector.to_array();
       xiv = xi_vector.to_array();
-      pu = xip + wp;
+      pu = xip * wp;
       wn = pu / 3.0e8;
       fprintf(output, "          XIP          WN           WL           PU           EV\n");
       fprintf(output, "     %13.4lE", xip);