> restart;
> with(plots):
# Das Programm untersucht Formeln fuer resonaante Wirkungsquerschnitte.
# Es ladet dazu ein, die Parameter zu variieren. Der Befehl
# restart hilft, wenn man sich verheddert hat. Zunaechst also die
# nichtrelativistische Behandlung (Formel 2.258 im Buch). Alle Energien
# sind in GeV einzusetzen. Der Faktor 389.4 rechnet 1/GeV^2 in mikrobarn
# um.
> fBW:=1/(4*(W-M_R)^2+Gamma^2);

                                      1
                     fBW := ---------------------
                                       2        2
                            4 (W - M_R)  + Gamma

> fakBW:=12*Pi/W^2*Gamma*Gamma_i*389.4;

                                   Pi Gamma Gamma_i
                   fakBW := 4672.8 ----------------
                                           2
                                          W

> sigBW:=fakBW*fBW;

                                   Pi Gamma Gamma_i
              sigBW := 4672.8 --------------------------
                               2             2        2
                              W  (4 (W - M_R)  + Gamma )

> p1:=plot(subs(M_R=2,Gamma=0.1,Gamma_i=0.00006,sigBW),W=1.5..2.5):
> display(p1);

# Nun die relativistischen Formeln (Formel 2.261 im Buch).
> frBW:=1/((W^2-M_R^2)^2+M_R^2*Gamma^2);

                                      1
                  frBW := --------------------------
                            2      2 2      2      2
                          (W  - M_R )  + M_R  Gamma

> fakrBW:=12*Pi*W^2*Gamma_i*Gamma/M_R^2*389.4;

                                      2
                                  Pi W  Gamma_i Gamma
                 fakrBW := 4672.8 -------------------
                                            2
                                         M_R

> sigrBW:=fakrBW*frBW;

                                      2
                                  Pi W  Gamma_i Gamma
          sigrBW := 4672.8 ---------------------------------
                              2    2      2 2      2      2
                           M_R  ((W  - M_R )  + M_R  Gamma )

> p2:=plot(subs(M_R=2,Gamma=0.1,Gamma_i=0.00006,sigrBW),W=1.5..2.5,color
> =blue):
> display(p2);

> display(p1,p2);

# Der Effekt ist also nicht sehr gross.