program p1

complex *16 complex1
complex *16 complex2
complex complex3

complex1=(1,2)
print *,'c1=',complex1
complex2=(5,10)

complex1=complex1+complex2
print *,'c1=',complex1

complex1=complex1*complex2
print *,'c1=',complex1

complex3=(25,10)
complex1=complex1/complex3
print *,'c1=',complex1

print *,'cos=',cos(complex1)
print *,'sin=',sin(complex1)

end program p1

 c1= (1.00000000000000,2.00000000000000)
 c1= (6.00000000000000,12.0000000000000)
 c1= (-90.0000000000000,120.000000000000)
 c1= (-1.44827586206897,5.37931034482759)
 cos= (13.2527369180607,107.621174445038)
 sin= (-107.625750880667,13.2521733884507)