~ubuntu-branches/ubuntu/karmic/scilab/karmic

« back to all changes in this revision

Viewing changes to routines/scicos/fsv.f

  • Committer: Bazaar Package Importer
  • Author(s): Torsten Werner
  • Date: 2002-03-21 16:57:43 UTC
  • Revision ID: james.westby@ubuntu.com-20020321165743-e9mv12c1tb1plztg
Tags: upstream-2.6
ImportĀ upstreamĀ versionĀ 2.6

Show diffs side-by-side

added added

removed removed

Lines of Context:
 
1
      subroutine fsv(flag,nevprt,t,xd,x,nx,z,nz,tvec,
 
2
     $        ntvec,rpar,nrpar,ipar,nipar
 
3
     $        ,u1,nu1,y1,ny1)
 
4
c
 
5
c
 
6
      double precision t,xd(*),x(*),z(*),tvec(*)
 
7
      integer flag,nevprt,nx,nz,ntvec,nrpar,ipar(*)
 
8
      double precision rpar(*),u1(*),y1(*)
 
9
c
 
10
c
 
11
      doubleprecision a,a0,b0,g
 
12
      a=u1(1)
 
13
      y1(1)=0.d0
 
14
      if(a.gt.1.0d0) then
 
15
         return
 
16
      endif
 
17
      g=1.4d0
 
18
      a0=2.0d0/g
 
19
      b0=(g+1)/g
 
20
      if(a.lt.0.528d0) then
 
21
         a=(2/(g+1.0d0))**(g/(g-1.0d0))  
 
22
      endif
 
23
      y1(1)=sqrt(2.0d0*g*(a**a0-a**b0)/(g-1.0d0))
 
24
   
 
25
 
 
26
c
 
27
c
 
28
      end