A NEW THEORY OF MATTER
(Updated 1/8-1994)

By O. Tedenstig
Idungatan 37
19551 Maersta/Sweden


THIS ARTICLE GIVES A NEW INSIGHT OF THE EXISTENCE OF ELEMEN-
TARY PARTICLES, BASED ON A NEW HYDROMECHANICAL MODEL OF
INTERACTION BETWEEN MATTER AND SPACE. THE EXISTENCE OF 
THE MUON, THE K, THE PROTON, THE TAU AND SO ON IS MOTIVATED
BY A VERY SIMPLE QUANTUM MECHANICAL PROCESS.


Today, very great resourses are spent in particle physics.
But current theories have many problems and are not able
to explain why elementary particles exist or why they behave
as they do.

Frequent report states that a final solution soon is to wait,
but that cannot inhibit the fact that existing theories
in many aspects are inconsistent and mainly erroneous.

Ever since the early 1970 the idea of quarks as basic
builing block of matter has been the leading idea. But free
quarks have never been found  in experiments or confirmed
in other ways.

And in addition to that, there exists a lot of very frequent
appearing elementary particles which not, according to theory,
contain no quarks at all.

But problems of particle physics is not only of technical or
theoretical art. The enormeous presistige and status of particle
physical research give no place for alternative ideas, even
if these ideas should offer a better platform for a solution.

But this problem is an inherent problem of the scientific society
and shall not be discussed here. Our aim is to try finding the
truth, unregarding how this truth looks like.

THE NEW THEORY
It will here be discussed an idea which may offer a new sight
on the existence of elementary particles. The idea do not
rest on the existence of quarks, but instead of a quantum
process between matter and the free space.

In this model, the main thought is that particle is a part of
the vacuum itself, but having a more local ordered structure
than the vacuum in the selected point.

This idea produces a set of base particles, being very
elementary and simple in shape. To these particles includes
the electron and the proton as to give some examples. The
other group of "elementary particles" are composed by these
base particles, hence being of complex or compound form.

The neutron is, to give an example, a compound particle.
This grouping of particles in basic, singular particles
and complex or compund particle forms give rise to a very
simple and fundamental description of matter.

MAIN IDEAS:
In our model, pointformed particles, like electrons for
instance, is regarded as a collection of vacuum field matter,
in an analogy with a gas or a raindrop in a cloud of rain.

External inertial forces of that gas collection is in a
steady equilibrium state by the vaccum field with density
1/Eo and velocity C (see reference 1).

All these particles create an electromagnetical field of its
own. The particle mass will oscillate into this force field
and create resonance conditions which will determine the particle
mass.

A BRIEF MATHEMATICAL ANALYSIS :
Let's regard a particle collection having an limiting area, A ,
and volume, V.

If oscillations occur into this gas mass (plasma) the volume
pressure product is an invariant entity. Applying the Boyle's law
on it, gives :

1) ================================

   po.Vo = p.V    ; p = po.(Vo/V)

===================================

The volumes Vo and V at the two states are related by the cube of
their particle radii, hence :

2) ===============================
                 3
    Vo/V = (Ro/R)

==================================

The total force which actuate on the particle surface from inside
or by impact from the outer vacuum field is :

3) ==================================

   Po = Ao.po    ; P = A.p

                   2           2
   where Ao = Ka.Ro  ; A = Ka.R

====================================


Hence, in each moment of the oscillating cycle of the particle
plasma the following force or tension is actuated :

4) ==========================================
   Peff= P-Po = A.p - Ao.po =

             2         3       2
         Ka.R.po.(Ro/R) - Ka.Ro.po =

              2
         Ka.Ro. (Ro/R-1).po = Ao.po.(Ro/R-1)
=============================================

The magnitude of the oscillation in the radius direction
will be :

5) ==============================================

   s= R-Ro = -(Ro/R -1).R  ; s/Ro = -(Ro/R -1).Ro
                            (if s/Ro is small)
=================================================

Inserting that in 4) gives :

6) ====================================
Peff =  - Ao.po.s/Ro
=======================================

This force will interact with the inertial mass of the particle
plasma, given from the second Newtonian law :

7) =======================
             2    2
    Pin = M.d s/dt

==========================

Peff and Pin are two forces which in each moment are in balance
with each other. That gives the balance equation :

8) ==============================
     2    2
    d s/dt  - (Po/M).s/Ro = 0

=================================

For making it able to solve this simple harmonical differential
equation, we must known the value of Po and how the particle
mass is related to the partilce's radius.

We begin by calcualting the value of Po.
In each oscillating periode an amount of mass, dm , is exchanged
between particle and field (with density 1/Eo, see reference 1).

That gives :

9) =====================
   dm = q.dV =  q.(A.s)
========================

By using the well known formula of relation between mass and total
              2
energy, E= m.c , this mass can be converted to energy, giving :

10) ==========================
             2            2
    dE = dm.c  = q.(A.s).c

==============================

In accord with Newton the product between force and distance is
energy, giving :

11) ======================================
                   2                 2
    P.s = q.(A.s).c    or   Po = Ao.c /Eo

==========================================

Hence, we have got an expression for Po expressed in parameters of
space and particle.

Mass density of all point particles, as here calculated is the
same. That give the following relation :

12) =======================
                  3
    M = me.(Ro/re)

===========================

where me, and re are parameters of electron.
From electromagnetic theory (reference 1) we know that Eo is the
same as 1/q, where q is the pseudo material density of the
vacuum field.  q is theregiven by the relation :

13) ==========================
             me
    q = ---------------- =1/Eo
               3   3
           8.Pi .re
===============================

Hence, calculating the relation Po/M gives :

14) ==============================================
                 2          3    3
    Po/M = ( Ao.c /Eo)/(me.Ro /re )

                  2  2      3  3           3  3
         =  (Ka.Ro .c /(8.Pi.re/me))/(me.Ro/re )

                    3          2
         = (Ka/(8.Pi ).(1/Ro).c  =

                    3           2
         = (Ka/(8.Pi )).(1/Ro).c

===================================================

Inserting in  8) gives :

15) ===============================================

      2    2            3           2
     d s/dt  - (Ka/(8.Pi )).(1/Ro).c .s/Ro = 0

===================================================

If s/Ro<<1, which we here can presume, the solution of that
differential equation will be :

16) =======================================
                           3      2  2
     Tr = 2.Pi. SQRT( (8.Pi/Ka).Ro /c )

===========================================

These vibration of the particle create disturbances in the
environment electromagnetic field. Then a resonance effect
is created between these vibrations and the vibrations in the
radius direction, claculated in 16).

The particle's own resonance periode is in the simplest way
calculated by the harmonical differential equation, giving
the pendulum solution.

17) ============================================
                                        2
    Tc = 2.Pi.SQRT(Mo.Ro/Fe)  ; Fe = me.c/re

                           2/3
    Tc = 2.Pi.re/c. (Mo/me)

=================================================

The resonance condition is also determined by a quantum number, n ,
being a whole integer relation value between the two
oscillations modes, giving :

18) =====================
    Tc = Tr.n
=========================

Together that gives :

19) =====================================================
                      2/3                   3      2  2
    2.Pi.re/c . (Mo/me)   = 2.Pi.SQRT( (8.Pi/Ka).Ro /c ).n
=========================================================

Using results from 12) then gives :

20) ==================================================
                 3
     M = me.(K.n)
                                                     2
     where the derived value of K = 2.51 for Ka= 4.Pi (the

     particle has torus shape), K = 4.44 for Ka= 4.Pi (the
     particle has spheric form)
     Because the exact interacting area not is known a
     qualified quess is that K lies somewhere between
     these values.
     ( Corrections done  23/9 -1993)
======================================================

where me is the electron mass, 0.51099906 Mev and n is an
integer value with beginning by n=1.
A table based on this relation is given below:

================================================================
 n      Particle           Name      Value of       Obser-
        mass in GeV                  K              vations
----------------------------------------------------------------
 1      0.01584            ?         Pi              Not observed
 2      0.10565839         u         2.957           data
 3      0.4936646          k         3.296           available
 4      0.93827231         p         3.062           today
 5      1.7841             Tau       3.034
 6      2.9796             nc        3.000
 7      5.2776             B         3.111
 8      8.3                zeta*     3.167

 ...
 17     81.0               W         3.184
 18     92.4               Z         3.142
 ---------------------------------------------------------------
 9      11.6               -         Pi             future
 10     15.8               -         Pi             observ-
 11     21.1               -         Pi             ations
 12     27.4               -         Pi
 13     34.8               -         Pi
 14     43.48              obs*      3.141
 15     53.5               -         Pi
 16     64.9               -         Pi
 --------------------------------------------------------------
 * See text.

The data for established particles were obtained from CERN
Particle Data Booklet, April 1988.

For n=8 see New Scientist 16 Aug. 1989, reporting a particle
named "zeta" at 8300 MeV, which is a better candidate than
T(1S).

For n=14, see New Scientist 25 May 1984, where a new
particle with mass 43.450 was reported. That value agrees
exactly with the given formula.

For n=9 to 13 see New Scientist 11 Sept. 1980, where
obswervation in the area 10 GeV to 35 GeV are indicated.
(See reference 2).


COMPLEX PARTICLE FORMS
If we for instance look at a neutron particle, an uncharged
particle with nearly the same mass as for the proton, being the
building block of the atomic nuclues, we see that this particle
has a slightly larger mass than the proton, its positive charged
counterpart. That slightly larger mass is created by a single
electron which move around the proton nuclues with a furious
velocity, increasing its rest mass by the factor 2.5 approximately.
That mass increase correspond with the well known physical law  :

                                        2  2
             m= mo/L, where L = SQRT(1-v /c ) where
             L is the Lorentz factor.

The mass increase effect was discovered in experiments 1901 by
Kaufmann. The physical reason to it we leave outside this
context but is presumed to be a pure electromagnetic effect.

The neutron system can here be regarded as a hydrogen system, but
where the electron move on a much higher energy level and near the
limit velocity of light in its own reference system. It's easy to
accept the neutron system in this way when we know that the neutron
annihilates in a proton and an electron plus extra energy (neutrinos).

Hence, the first effect of the orbiting particle in such a system
is that the orbiting particle's mass increases with velocity, giving

21) ====================================
                                 2  2
    M = Mo/L    ; L = SQRT( 1 - v /c )

========================================

We solve out the orbiting particle's velocity :

22) ===================================
                      2
    v = c.SQRT ( 1 - A )   ; A = Mo/M

=======================================

If we for simplicity define the orbit as a pure circular movement
(only is partly true), we from Newton's mass inertial laws get :

23) ======================
            2
    Fc = M.v/D1

==========================

where D1 is the orbiting distance related to the center of the
particle system.

There are 3 kind of forces which actuate the orbiting particle,
the strong force, the inertial force and the electromagnetic force.
The "strong force" is created by the vacuum pressure which is in
accord with the electromagnetic theory (see reference 1).

24) ========================================

               2           2         2  2
     Fmax = q.C .A  = (me.c /re) . R /re

============================================

The real strong force also is dependent by the distance to the
central particle, D , hence giving the effective working force :

25) =======================================================
                        2        2          2 3   2  2  2
     Fs = Fmax . 4.Pi.Rc /(4.Pi.D )  = (me.c/re).R .Rc/D

============================================================

Using the hypothesis, that all pointformed particles have invariant
mass density, this equation can be rewritten to :

26) ======================================
             2       2    2/3   2/3  4/3
    Fs = me.c.re.(1/D ).Mc .  M  / me

==========================================

If the central particle, as here for instance in an neutron,
is very heavy in relation to the orbiting particle, the common
rotation radius, D1 , is the same as the distance between the
particles. But in a system where the orbiting particle is more
heavy, there will be a discrepance between these entities :

27) =================================================

    D1 = z. D    where z =  1/(1 + M/Mc)

where D=2*D1 when M=Mc.
=====================================================

There exists also electromagnetical forces between the particles
into the system, but these forces are so small related to Fs and
Fc , so they can be neglected here.

The electromagnetic field instead has another more important
function. The existence of this field will give rise to
a quantum mechanical resonance effect between these oscillating
conditions present into the system. From our base particles we
have the follwoing set of base frequencies (or time periodes) :


28) ===============================================
                                         2/3
1)  To1 = 2.Pi.re/c . K1   ;  K1 = (M/me)
                                          1/3
2)  To2 = 2.Pi.re/c . K2    ;  K2 = (M/me) . 2. Pi
                                           2/3
3)  Tc3 = 2.Pi.re/c . K3    ;  K3 = (Mc/me)
                                           1/3
4)  Tc4 = 2.Pi.re/c . K4    ;  K4 = (Mc/me). 2.P
====================================================

Hence in a system of 2 particles, as in a neutron , there
will be 4 different possiblilities of the orbiting particle
to move in time to the electromagnetic field variations.

For instance, for a neutron system in the lowest resoance
state, the system will be a neutron, but in the highest state
it will be the LAMBDA particle, being a neutral but having a
mass 2170 times the electron mass.

The corresponding factor for the neutron is 1838.12.
These electromagnetic vibration create disturbances on the
orbiting particle and forces it to rotate in time to these
changes (a resoance or quantum effect).

In the same way as in the atomic system ,this relation will
be :

29) ==============================

      2.Pi.D1/v = Tk.N

==================================

where N is the number of orbiting particles of same mass in the
same orbit.

The solution of the total mass of the particle system will be
when the forces Fc and Fs are in balance.

30) ============================================================

a)   A = Mo/M                 Mo orbit particle mass in rest
                              M orbit particle mass when moving
                      2
b)   v = c.SQRT (1 - A )      v is the orbit particle velocity
                              in its own reference system
             2
c)   Fc = M.v /D1             The centrifugal force of the orbiting
                              particle

              2    2    2*k/3  2*k/3 4/3
d)   Fs = me.c.re/D . Mc   .  M /  me     The strong force acting
                                          on the orbiting particle
     (the k-factor, see definition below)


e)   D1 = z.D ; z = 1/(1+M/Mc)             Discrepancy between the
                                           particle distance and
                                           the common orbiting
                                           distance

f)   Fc - Fs = 0              The balance condistion for a stabele
                              system.


g)   D1 = Tk.N.v/(2.Pi)


h)   Tk = Toc,Tor,Tcc or Tcr   see (28)

g)   Mtot = Mc + N*M

(The k-factor has been introduced for accomodating the copling
factor between the orbit and center particle, which can vary
slightly in different systems, see table below)

 =============================================================

This set of equations cannot easily be sovled in a convinient way by
analogous calculations. However, by a recursive computer process the
balance point is easily calculated.

In experiments, mass of particles are determined by a distribution
function where the top medium value is determined. Hence there
are some distributions in the measuring values which at the end
will determine the official particle mass. The Gaussian mean
mass value is sanctioned as the official mass value.

To make it possible identifying and comparing values from this
model and experimental founds, we have introduced a correction
factor "k" in the listing below. The k-factor has been
introduced into the model in formula 30d) above in the
exponential factors with the idial value equal to k=1.

Some of the most easy identificable compound particles areas
calculated in the table below.

Mc  Mo  n  State  Mass Mev   Symbole,name  Charge  k
---------------------------------------------------------------
u   e   1  1      0.13497    Pi-zero       0      1.15
           2      "          Pi-zero       0      1.22
           3                 not observed
           4                 nor observed
---------------------------------------------------------------
u   e   2  1      0.139568   Pi +-         +-     1.06
           2         "       Pi +-         +-     1.18
           3                 not observed
           4                 not observed
---------------------------------------------------------------
K   e   1  1      0.5488     eta zero      0      1.05
           2         "       eta zero      0      1.07
           3      0.49767    K zero        0      1.06
           4         "       K zero        0      1.02
--------------------------------------------------------------
K    e  2  1-4                         *)  see comment below
--------------------------------------------------------------
K    u  1  1-4    0.6        epsilon       ?0     1.0   rare
--------------------------------------------------------------
K    u  2  1-4    0.770-0.783  rho,ohmega  0+-    1.05
--------------------------------------------------------------
K    K  1  1-4    1.020       Phi          0      1.00-1.05
--------------------------------------------------------------
p    e  1  1      1.11563     Lambda zero  0      1.135
           2         "            "        0      1.12
           3      0.93956563  neutron      0      1.00
           4         "            "        0      0.95
-------------------------------------------------------------
p    e  2  1-4                   *) see comment below
-------------------------------------------------------------
p    u  1  1-4    1.11563     Lambda zero  0      1.05
-------------------------------------------------------------
p    u  2  1-4    1.18937     Sigma        +-0    0.98
-------------------------------------------------------------
p    u  3  1-4    1.3149      Xsi          +-0    0.98
-------------------------------------------------------------
p  u+K  1  1-4    1.67243     Ohmega       +-   medium p+u.p+K
-------------------------------------------------------------
T    e  1  1      1.8645      D zero       0      1.06
           2        "           "          0      1.06
           3        "           "          0      1.05
           4        "           "          0      1.05
-------------------------------------------------------------
T    e  2  1      1.8693      D +-        +-      1.04
           2        "           "         +-      1.04
           3        "           "         +-      1.04
           4        "           "         +-      1.04
-------------------------------------------------------------
T+(e)+u 1  1-4    1.971       F +-0       +-0     1.05
-------------------------------------------------------------
T    u  2  1-4    2.1127      Ds +-       +-      1.04
-------------------------------------------------------------

*)   It seems that systems containing 2 orbit electrons are
     not so probable

**)  A compound particle where the orbit mass is large related
     to the center particle will be short lived, and rare.


References : 1) A new model of interaction between matter
                and vaccum, Galilean ElectroDynamics, Vol 4,
                no 4 by the author.

             2) Elementary-Particles masses predicted, Galilean
                ElectroDynamics Vol 2, no 3, by the author

             3) Particle Properties Data Booklet, Technical
                information Depratmetn, MS 90-2125, Lawrence
                Berley, CA 94720 USA

             4) A new way to physics 1990, Ove Tedenstig
                (paperback 500 pages).

From:	ERA6::OTG          23-SEP-1993 13:29:26.33
To:	ERA5::ERAOTG
CC:	
Subj:	

********* Program for calculating singular particles ***********
****************************************************************

Program p(input,output);
uses    CRT;
var     n:integer;
        k,PI,m,x:real;
        q:char;

procedure part_mass;
begin

if ((n>1.8) and (n<2.2)) then
writeln('Reference mass u-on 105.65839 Mev = 206.77 emu ');

if ((n>2.8) and (n<3.2)) then
writeln('K-on 493.646 Mev= 966.06 emu ');

if ((n>3.8) and (n<4.2)) then
writeln('proton 938.27231 Mev = 1836.12 emu ');

if ((n>4.8) and (n<5.2)) then
writeln('Tau-on 1784.1 Mev = 3491.5 emu ');

if ((n>5.8) and (n<6.2)) then
writeln(' nc 2979.6 Mev = 5831 emu ');

if ((n>6.2) and (n<7.2)) then
writeln(' B 5277.6 = 10328.2 emu ');

if ((n>7.2) and (n<8.2)) then
writeln(' zeta 9400 Mev= 16243 emu ');

end;


Procedure present;
begin
for n:=1 to 25 do begin
writeln;
end;
writeln('Written by O.Tedenstig/M„rsta Sweden 1991 ');
writeln('This program calculates the mass of pointformed and always ');
writeln('charged particles like the u-on, the K-on, the proton, the ');
writeln('Tau-on, the nc, the B and so on unlimitely ');
writeln('Each such particle has a unique quantum number n=1,2.... which ');
writeln('is inserted at program start');
writeln('The theoretical formula is not exact, therefore the K-value');
writeln('of the formula m=mo.(K.n)**3 has to be slightly adjused ');
writeln;
writeln('The electron mass is used as reference mass= 1 mass unit ');
writeln('The electron has a mass/energy of 0.51099 Mev, Mega electron volts');
writeln('The electron reference mass = 1 emu ');
writeln('THE PROGRAM IS STOPPED BY PRINTING 0 AS INPUT')
end;


procedure define;
begin
pi:=4*arctan(1);
write('Present standard masses y/n ');readln(q);

end;

procedure write_out;
begin
writeln('Calculated mass= ',m,' emu' );
part_mass;
writeln;

end;

procedure compute;
begin
M:=exp(3*ln(K*n));
end;

begin
present;
define;k:=1;
if q='y' then begin part_mass;end;

while K<>0  do begin
write('Insert quantumnumber 1,2,3,4,5,6,7.... ');readln(n);
write('Insert value near = pi (theoretical value) ');readln(K);
writeln;
compute;
write_out;
end;
readln;
end.



From:	ERA6::OTG          23-SEP-1993 13:29:02.17
To:	ERA5::ERAOTG
CC:	
Subj:	


************************ program for complex particles ************
*******************************************************************

program	pkomplex;

VAR	corr,r0d,z1,m0,mc,mo,xo,xc,yo,yc,x,B1,B2,B3:real;
	mcc,masstabel,y,eo,ec,vc,r,z,k,a,d,fx,fe,fo,fc,korr:real;
	nc,no,n,n1,n12,s,start,nn:integer;
	b,c,o:packed array(.1..3.) of char;
  test:boolean;
CONST	pi=3.1415926;
	e=0.511003414E-3;
	ev=0.511003414;
	me=9.1066E-31;

procedure lines;
begin
for nn:=1 to 25 do begin
writeln;
end;
end;

Procedure explain_terms;
begin
writeln('The program computes the mass of a complex particle');
writeln('according to a new model for elementary particles ');
writeln('consisting of a center mass and one or several orbital masses');
writeln('Input 0 for electron, 0.51099906 Gev or 1 electron masses');
writeln('Input 2 for myon 0.10565839 Gev or 206.77 electron masses');
writeln('Input 3 for K-on 0.493646 Gev or 966 electron masses');
writeln('Input 4 for proton, 0.93827232 Gev or 1836.12 electron masses');
writeln('Input 5 for Tauon, 1.7841 Gev or 3491.4 electron masses');
writeln('Input 6 for nc, 2.980 Gev or 5832 electron masses ');
writeln('Input 7 for B, 5.2776 Gev or 10327 electron masses ');
writeln('Input 8 for z, 8.3 Gev or 16243 electron masses');
writeln;
writeln('These base particles approximately follows the ');
writeln('mathematical relation  M=electron_mass*(Pi*n)**3 ');
writeln;
end;


PROCEDURE calc(s,n:integer;var r,z1,mo,vc,r0d:real);
	begin
	if a=0 then mo:=m0;
	if r=0 then r:=0.00000001;
	mo:=mo*(1+z1);a:=m0/mo;
  vc:=exp(0.5*ln(1-a*a));
	z:=1;
	if n=1 then z:=1/(1+mo/mc);

  xo:=exp(eo*ln(mo));
  xc:=exp(ec*ln(mc));
  yo:=exp(y*ln(mo));
  yc:=exp(x*ln(mc));

	case s of
  1:k:=exp(2*ln(xo));
  2:k:=2*pi*(xo);
  3:k:=exp(2*ln(xc));
	4:k:=2*pi*xc;
	end;

	r0d:=z/(vc*k*n);
	fo:=exp(2*ln(r0d))*exp(4*ln(yo))*(n-1);
	fe:=exp(2*ln(r0d));
	fx:=exp(2*ln(r0d))*exp(2*ln(yo))*exp(2*ln(yc));
	fc:=r0d*mo*exp(2*ln(vc))/z;
	r:=fc-fx-fe-fo;
	end;(*calc*)

procedure information;
begin
if (n=0) then begin write('Write center particle type 2,3,4,5... ');readln(nc);end;
	write('Write orbit particle type 0,2,3,4.... ');readln(no);
  write('Number of orbits 1,2,3...             ');readln(n);writeln;
	corr:=1;
if (no=0) then m0:=1;(*electron mass*)
if (no=1) then m0:=exp(3*ln(pi*1));
if (no=2) then m0:=206.77;(*muon mass*)
if (no=3) then m0:=966.6;(*kaon-mass*)
if (no=4) then m0:=1836.12;(*proton-mass*)
if (no=5) then m0:=1784.2/ev;(*tau-mass*)
if (no=6) then m0:=2978/ev;(*eta-mass*)
if (no>6) then m0:=exp(3*ln(pi*nc));

if (nc=1) then mc:=exp(3*ln(pi*n));
if (nc=2) then mc:=206.77;(*muon-mass*)
if (nc=3) then mc:=966.6;(*kaon-mass*)
if (nc=4) then mc:=1836.12;(*proton-mass*)
if (nc=5) then mc:=1784.2/ev;(*tau-mass*)
if (nc=6) then mc:=2978/ev;(*eta-mass*)
if (nc>6) then mc:=exp(3*ln(pi*nc));
writeln('The computed particle (particles) consist of :');
write('The center mass being ');
if nc=2 then write('a myon ');
if nc=3 then write('a K-on ');
if nc=4 then write('a proton ');
if nc=5 then write('a Tau-on');
if nc=6 then write('a nc-particle');
if nc=7 then write('a B-meson ');
if nc=8 then write('a zeta-particle ');
write(' with ',n,' orbitals of ');
if no=0 then writeln('electrons');
if no=1 then writeln('unknown base particle ');
if no=2 then writeln('my-ons');
if no=3 then writeln('K-ons');
if no=4 then writeln('protons');
Write('The particle(s) is(are)  ');
if ((n=1) or (n=3) or (n=5)) then writeln('neutral uncharged ');
if ((n=2) or (n=4) or (n=6)) then
writeln('positive or negative charged particle');
writeln;
writeln('Computing started, please wait a moment');
end;

procedure result;
begin
test:=false;
if (test=false) then write('    ',s);
if (test=false) then write('              ',mcc,'            ');
if (test=false) then writeln(mcc*ev*1E-3:1:2,' Gev');
if (s<=2) then B1:=b1+mcc;
if (s>2) then B2:=b2+mcc;
b3:=b3+mcc;
mo:=m0;r:=-1;z1:=0.1;
end;

procedure parameters;
begin
write('corr= (theoretical value=1) ');readln(corr);
n:=0;n12:=2;
B1:=0;B2:=0;B3:=0;
eo:=1/3;(*standard value*)
ec:=1/3;(*standard value*)
x:=corr/3;y:=corr/3;
z1:=0.1;
end;

procedure print_medium;
begin
write('mass_medium 1_2  ',B1/n12,'    ');
writeln('  ',1E-3*ev*b1/n12,' Gev');b1:=0;
if (test=false) then write('mass_medium 3_4   ',1E-3*b2/2,'    ');
if (test=false) then writeln('  ',1E-3*ev*b2/2,' Gev');b2:=0;
if (test=false) then write('mass_medium_total ',1E-3*b3/4,'    ');
if (test=false) then writeln(' ',1E-3*ev*b3/4,' Gev');b3:=0;
end;

(************************ START PROGRAM *******************************)
begin
lines;
explain_terms;
parameters;
information;
begin
begin
for s:=1 to 4 do
begin
a:=0;r:=-1;z1:=0.1;
while (z1>0.0000001) do begin
while r<0 do begin
calc(s,n,r,z1,mo,vc,r0d);
end;(*while r-loop*)
mo:=mo*(1-z1);
z1:=z1/1.1;r:=-1;
end;(*while z1-loop*)
mcc:=n*mo+mc;
result;
end;
end;
print_medium;
end;
writeln('ready');
end.

(**** en ****)














