c****************************************************************

	program masses
	
c	calcul du facteur B=Mb/(Ma+Mb)


c	Auteur: P. Morel, Departement J.D. Cassini, O.C.A., Observatoire de Nice

c---------------------------------------------------------

	use ctephy
	
	implicit none
	
	double precision :: radian,UA=1.495978706d13
	real :: para,a,i,Ka,P,e,Ma,Mb,cte1=6.283,B,sum
	
	character (len=15) :: orbite,parallaxe,vitesse
		
c----------------------------------------------------------

	radian=pi/180.d0 
	
	orbite='Soderhjelm' ; parallaxe='Hipparcos' ; vitesse='Scarfe et al.'
	para=0.0926 ; a=1.33 ; i=131. ; Ka=4.01 ; P=34.45 ; e=0.46
	
	i=i*radian ; Ka=Ka*1.d5/UA*secon6*1.d-6
	B=para*Ka*P*sqrt(1.-e**2)/cte1/a/abs(sin(i))
	
	sum=(a/para)**3/P**2 ; Mb=B*sum ; Ma=sum-Mb
	
	print*
	write(6,1)orbite(:len_trim(orbite)),parallaxe(:len_trim(parallaxe)),
	1	vitesse(:len_trim(vitesse))
1	format(1x,'orbite: ',a,', parallaxe: ',a,', vitesse: ',a)
	write(6,2)a,i,P,e
2	format(1x,'a=',1pd10.3,', i=',1pd10.3,', P=',1pd10.3,', e=',1pd10.3)
	write(6,3)para,Ka
3	format(1x,'pi=',1pd10.3,' Ka=',1pd10.3)		
	write(6,4)B,sum,Ma,Mb
4	format(1x,'B=',1pd10.3,', Ma+Mb=',1pd10.3,', Ma=',1pd10.3,
	1	', Mb=',1pd10.3)

	print*		
	orbite='BaizeII' ; parallaxe='Lippincott' ; vitesse='Scarfe et al.'
	para=0.098 ; a=1.355 ; i=132.9 ; Ka=4.01 ; P=34.487 ; e=0.460
	
	i=i*radian ; Ka=Ka*1.d5/UA*secon6*1.d-6
	B=para*Ka*P*sqrt(1.-e**2)/cte1/a/abs(sin(i))
	sum=(a/para)**3/P**2 ; Mb=B*sum ; Ma=sum-Mb
	

	write(6,1)orbite(:len_trim(orbite)),parallaxe(:len_trim(parallaxe)),
	1	vitesse(:len_trim(vitesse))
	write(6,2)a,i,P,e
	write(6,3)para,Ka
	write(6,4)B,sum,Ma,Mb

	print*		
	orbite='BaizeII' ; parallaxe='Lippincott(rel)' ; vitesse='Scarfe et al.'
	para=0.0914 ; a=1.355 ; i=132.9 ; Ka=4.01 ; P=34.487 ; e=0.460
	
	i=i*radian ; Ka=Ka*1.d5/UA*secon6*1.d-6
	B=para*Ka*P*sqrt(1.-e**2)/cte1/a/abs(sin(i))
	sum=(a/para)**3/P**2 ; Mb=B*sum ; Ma=sum-Mb
	

	write(6,1)orbite(:len_trim(orbite)),parallaxe(:len_trim(parallaxe)),
	1	vitesse(:len_trim(vitesse))
	write(6,2)a,i,P,e
	write(6,3)para,Ka
	write(6,4)B,sum,Ma,Mb

	end