clc
a=2500
//doing for the first values only
Bf=4
Bp=0.305
q=a/Bf^2
Sep=4
Sef=Sep*(2*Bf/(Bf+Bp))^2
printf('Sef = %f mm',Sef)