function p=p(xa,xb,xc,xd,xf)
p=(sqrt(-((2.5992e-26)*(1/(0.081*xc + 0.133361783847*xb) +0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*(1/(0.0666808919236*xa +1/xf) + 0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*xa*xb + 24*(1/(0.081*xc + 0.133361783847*xb)+ 0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*((2.736e-27)*xc +(9.5304e-28)*xa)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa))*xb - 24*(1/(0.081*xc + 0.133361783847*xb) +0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*((5.472e-27)*xc +(1.90608e-27)*xa)*((-0.04*xc/(0.081*xc + 0.133361783847*xb) -0.000355707307818*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 1/(0.081*xc + 0.133361783847*xb) +0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*xb +24*(1/(0.0666808919236*xa + 1/xf) +0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*((1.368e-27)*xc +(4.7652e-28)*xa)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa))*xa + 24*((2.4e-12)*xc +(8.36e-13)*xa)*((1.5625e-17)*xd^2 +(2.85e-14)*xb)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)) + 24*((2.96875e-29)*xd^2 +(5.415e-26)*xb)*(1/(0.081*xc + 0.133361783847*xb) +0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*xb +(2.28e-11)*((1.5625e-17)*xd^2 + (2.85e-14)*xb)*(1/(0.0666808919236*xa+ 1/xf) + 0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*xa - 24*(-0.00266723567694*xc*xb/((0.081*xc+ 0.133361783847*xb)*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)*(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa)) + 0.0666808919236*xb/((0.081*xc +0.133361783847*xb)*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*((1.9204096874e-14)*xc +(6.68942707778e-15)*xa)*((2.4e-14)*xd +(1.672e-14)*xb)*xb/((0.081*xc +0.133361783847*xb)*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) + 24*(1/(0.081*xc +0.133361783847*xb) + 0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*((1.368e-27)*xd +(9.5304e-28)*xb)*xb/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb) +(6.84e-13)*(1/(0.0666808919236*xa + 1/xf) +0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*((2.4e-14)*xd +(1.672e-14)*xb)*xa/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb) -24*(5.33447135389*(-0.001*xc/(0.081*xc + 0.133361783847*xb) -(8.89268269546e-06)*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*xa/((0.0666808919236*xa +1/xf)*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)) + 0.0666808919236*xa/((0.0666808919236*xa +1/xf)*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*((9.12194601515e-29)*xc +(3.17747786194e-29)*xa)*xa^2/((0.0666808919236*xa +1/xf)*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)) + 24*((7.2e-14)*xc +(2.508e-14)*xa)*((2.4e-14)*xd +(1.672e-14)*xb)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa))/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb) + 24*((1.2e-12)*xd +(8.36e-13)*xb)*((1.5625e-17)*xd^2 +(2.85e-14)*xb)/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb) - 24*((1.21651899274e-29)*xd+ (8.47508231608e-30)*xb)*xb^3/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)^2))/((1.14e-13)*(1/(0.081*xc +0.133361783847*xb) + 0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*xb +(5.7e-14)*(1/(0.0666808919236*xa + 1/xf) +0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*xa + ((1.44e-13)*xc +(5.016e-14)*xa)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)) + (1.5625e-15)*xd^2 + (2.85e-12)*xb +((7.2e-14)*xd + (5.016e-14)*xb)/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb))^2 + 1) -1)*((1.14e-13)*(1/(0.081*xc + 0.133361783847*xb) +0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*xb +(5.7e-14)*(1/(0.0666808919236*xa + 1/xf) +0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*xa + ((1.44e-13)*xc +(5.016e-14)*xa)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)) + (1.5625e-15)*xd^2 + (2.85e-12)*xb +((7.2e-14)*xd + (5.016e-14)*xb)/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb))/((2.166e-27)*(1/(0.081*xc+ 0.133361783847*xb) + 0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*(1/(0.0666808919236*xa +1/xf) + 0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*xa*xb + 2*(1/(0.081*xc + 0.133361783847*xb)+ 0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*((2.736e-27)*xc +(9.5304e-28)*xa)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa))*xb - 2*(1/(0.081*xc + 0.133361783847*xb) +0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*((5.472e-27)*xc +(1.90608e-27)*xa)*((-0.04*xc/(0.081*xc + 0.133361783847*xb) -0.000355707307818*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 1/(0.081*xc + 0.133361783847*xb) +0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*xb +2*(1/(0.0666808919236*xa + 1/xf) +0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*((1.368e-27)*xc +(4.7652e-28)*xa)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa))*xa + 2*((2.4e-12)*xc +(8.36e-13)*xa)*((1.5625e-17)*xd^2 +(2.85e-14)*xb)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)) + 2*((2.96875e-29)*xd^2 +(5.415e-26)*xb)*(1/(0.081*xc + 0.133361783847*xb) +0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*xb +(1.9e-12)*((1.5625e-17)*xd^2 + (2.85e-14)*xb)*(1/(0.0666808919236*xa+ 1/xf) + 0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*xa - 2*(-0.00266723567694*xc*xb/((0.081*xc+ 0.133361783847*xb)*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)*(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa)) + 0.0666808919236*xb/((0.081*xc +0.133361783847*xb)*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*((1.9204096874e-14)*xc +(6.68942707778e-15)*xa)*((2.4e-14)*xd +(1.672e-14)*xb)*xb/((0.081*xc +0.133361783847*xb)*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) + 2*(1/(0.081*xc +0.133361783847*xb) + 0.00889268269546*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*((1.368e-27)*xd +(9.5304e-28)*xb)*xb/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb) +(5.7e-14)*(1/(0.0666808919236*xa + 1/xf) +0.00444634134773*xa^2/((0.0666808919236*xa +1/xf)^2*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*((2.4e-14)*xd +(1.672e-14)*xb)*xa/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb) -2*(5.33447135389*(-0.001*xc/(0.081*xc + 0.133361783847*xb) -(8.89268269546e-06)*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)))*xa/((0.0666808919236*xa +1/xf)*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)) + 0.0666808919236*xa/((0.0666808919236*xa +1/xf)*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)))*((9.12194601515e-29)*xc +(3.17747786194e-29)*xa)*xa^2/((0.0666808919236*xa +1/xf)*(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa)) + 2*((7.2e-14)*xc +(2.508e-14)*xa)*((2.4e-14)*xd +(1.672e-14)*xb)*((-0.08*xc/(0.081*xc + 0.133361783847*xb) -0.000711414615637*xc*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) +0.0666808919236*xb)))/(-0.00444634134773*xa^2/(0.0666808919236*xa +1/xf) + 0.0666808919236*xa) + 2/(0.081*xc + 0.133361783847*xb) +0.0177853653909*xb^2/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)) +1/(-0.00444634134773*xa^2/(0.0666808919236*xa + 1/xf) +0.0666808919236*xa))/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb) + 2*((1.2e-12)*xd +(8.36e-13)*xb)*((1.5625e-17)*xd^2 +(2.85e-14)*xb)/(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb) - 2*((1.21651899274e-29)*xd+ (8.47508231608e-30)*xb)*xb^3/((0.081*xc +0.133361783847*xb)^2*(-0.00889268269546*xb^2/(0.081*xc +0.133361783847*xb) + 0.0666808919236*xb)^2))
endfunction