       subroutine plasma(flam,rho,mue,cv,cvp,ca,cap,n,qplas)

       implicit double precision (a-h,o-z)

       dimension a(0:2),b(3)
       double precision mue
       xi=(((rho/mue)/1d9)**(0.33333333333d0))/flam
       a(0)=2.320d-7
       a(1)=8.449d-8
       a(2)=1.787d-8
       b(1)=2.581d-2
       b(2)=1.734d-2
       b(3)=6.990d-4
       c=0.56457d0
       call qcalc(a,b,c,xi,flam,f)
       qplas=(cv**2+n*cvp**2)*((rho/mue)**3)*f

       return
       end

************************************************************************

