当前位置:文档之家› Fortran90数值计算chap2c2f9

Fortran90数值计算chap2c2f9

Fortran90数值计算chap2c2f9
Fortran90数值计算chap2c2f9

Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).C2.Alphabetical Listing of Explicit Interfaces The?le supplied as nr.f90contains explicit interfaces for all the Numerical Recipes routines(except those already in the module nrutil).The interfaces are in alphabetical order,by the generic interface name,if one exists,or by the speci?c routine name if there is no generic name.The?le nr.f90is normally invoked via a USE statement within a main program or subroutine that references a Numerical Recipes routine.See§21.1for an example.MODULE nr INTERFACE SUBROUTINE airy(x,ai,bi,aip,bip)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP),INTENT(OUT)::ai,bi,aip,bip END SUBROUTINE airy END INTERFACE INTERFACE SUBROUTINE amebsa(p,y,pb,yb,ftol,func,iter,temptr)USE nrtype INTEGER(I4B),INTENT(INOUT)::iter REAL(SP),INTENT(INOUT)::yb REAL(SP),INTENT(IN)::ftol,temptr REAL(SP),DIMENSION(:),INTENT(INOUT)::y,pb REAL(SP),DIMENSION(:,:),INTENT(INOUT)::p INTERFACE FUNCTION func(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP)::func END FUNCTION func END INTERFACE END SUBROUTINE amebsa END INTERFACE INTERFACE SUBROUTINE amoeba(p,y,ftol,func,iter)USE nrtype INTEGER(I4B),INTENT(OUT)::iter REAL(SP),INTENT(IN)::ftol REAL(SP),DIMENSION(:),INTENT(INOUT)::y

REAL(SP),DIMENSION(:,:),INTENT(INOUT)::p INTERFACE

FUNCTION func(x)

USE nrtype

REAL(SP),DIMENSION(:),INTENT(IN)::x

REAL(SP)::func

END FUNCTION func

END INTERFACE

1384

Appendix C2.Alphabetical Listing of Explicit Interfaces1385 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).END SUBROUTINE amoeba END INTERFACE INTERFACE SUBROUTINE anneal(x,y,iorder)USE nrtype INTEGER(I4B),DIMENSION(:),INTENT(INOUT)::iorder REAL(SP),DIMENSION(:),INTENT(IN)::x,y END SUBROUTINE anneal END INTERFACE INTERFACE SUBROUTINE asolve(b,x,itrnsp)USE nrtype REAL(DP),DIMENSION(:),INTENT(IN)::b REAL(DP),DIMENSION(:),INTENT(OUT)::x INTEGER(I4B),INTENT(IN)::itrnsp END SUBROUTINE asolve END INTERFACE INTERFACE SUBROUTINE atimes(x,r,itrnsp)USE nrtype REAL(DP),DIMENSION(:),INTENT(IN)::x REAL(DP),DIMENSION(:),INTENT(OUT)::r INTEGER(I4B),INTENT(IN)::itrnsp END SUBROUTINE atimes END INTERFACE INTERFACE SUBROUTINE avevar(data,ave,var)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::data REAL(SP),INTENT(OUT)::ave,var END SUBROUTINE avevar END INTERFACE INTERFACE SUBROUTINE balanc(a)USE nrtype REAL(SP),DIMENSION(:,:),INTENT(INOUT)::a END SUBROUTINE balanc END INTERFACE INTERFACE SUBROUTINE banbks(a,m1,m2,al,indx,b)USE nrtype INTEGER(I4B),INTENT(IN)::m1,m2INTEGER(I4B),DIMENSION(:),INTENT(IN)::indx REAL(SP),DIMENSION(:,:),INTENT(IN)::a,al REAL(SP),DIMENSION(:),INTENT(INOUT)::b END SUBROUTINE banbks END INTERFACE INTERFACE SUBROUTINE bandec(a,m1,m2,al,indx,d)USE nrtype INTEGER(I4B),INTENT(IN)::m1,m2INTEGER(I4B),DIMENSION(:),INTENT(OUT)::indx REAL(SP),INTENT(OUT)::d REAL(SP),DIMENSION(:,:),INTENT(INOUT)::a REAL(SP),DIMENSION(:,:),INTENT(OUT)::al END SUBROUTINE bandec

END INTERFACE

INTERFACE

SUBROUTINE banmul(a,m1,m2,x,b)

USE nrtype

INTEGER(I4B),INTENT(IN)::m1,m2

REAL(SP),DIMENSION(:),INTENT(IN)::x

REAL(SP),DIMENSION(:),INTENT(OUT)::b

REAL(SP),DIMENSION(:,:),INTENT(IN)::a

1386Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).END SUBROUTINE banmul END INTERFACE INTERFACE SUBROUTINE bcucof(y,y1,y2,y12,d1,d2,c)USE nrtype REAL(SP),INTENT(IN)::d1,d2REAL(SP),DIMENSION(4),INTENT(IN)::y,y1,y2,y12REAL(SP),DIMENSION(4,4),INTENT(OUT)::c END SUBROUTINE bcucof END INTERFACE INTERFACE SUBROUTINE bcuint(y,y1,y2,y12,x1l,x1u,x2l,x2u,x1,x2,ansy,&ansy1,ansy2)USE nrtype REAL(SP),DIMENSION(4),INTENT(IN)::y,y1,y2,y12REAL(SP),INTENT(IN)::x1l,x1u,x2l,x2u,x1,x2REAL(SP),INTENT(OUT)::ansy,ansy1,ansy2END SUBROUTINE bcuint END INTERFACE INTERFACE beschb SUBROUTINE beschb_s(x,gam1,gam2,gampl,gammi)USE nrtype REAL(DP),INTENT(IN)::x REAL(DP),INTENT(OUT)::gam1,gam2,gampl,gammi END SUBROUTINE beschb_s SUBROUTINE beschb_v(x,gam1,gam2,gampl,gammi)USE nrtype REAL(DP),DIMENSION(:),INTENT(IN)::x REAL(DP),DIMENSION(:),INTENT(OUT)::gam1,gam2,gampl,gammi END SUBROUTINE beschb_v END INTERFACE INTERFACE bessi FUNCTION bessi_s(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),INTENT(IN)::x REAL(SP)::bessi_s END FUNCTION bessi_s FUNCTION bessi_v(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessi_v END FUNCTION bessi_v END INTERFACE INTERFACE bessi0FUNCTION bessi0_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::bessi0_s END FUNCTION bessi0_s FUNCTION bessi0_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessi0_v

END FUNCTION bessi0_v

END INTERFACE

INTERFACE bessi1

FUNCTION bessi1_s(x)

USE nrtype

REAL(SP),INTENT(IN)::x

REAL(SP)::bessi1_s

END FUNCTION bessi1_s

Appendix C2.Alphabetical Listing of Explicit Interfaces1387 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).FUNCTION bessi1_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessi1_v END FUNCTION bessi1_v END INTERFACE INTERFACE SUBROUTINE bessik(x,xnu,ri,rk,rip,rkp)USE nrtype REAL(SP),INTENT(IN)::x,xnu REAL(SP),INTENT(OUT)::ri,rk,rip,rkp END SUBROUTINE bessik END INTERFACE INTERFACE bessj FUNCTION bessj_s(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),INTENT(IN)::x REAL(SP)::bessj_s END FUNCTION bessj_s FUNCTION bessj_v(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessj_v END FUNCTION bessj_v END INTERFACE INTERFACE bessj0FUNCTION bessj0_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::bessj0_s END FUNCTION bessj0_s FUNCTION bessj0_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessj0_v END FUNCTION bessj0_v END INTERFACE INTERFACE bessj1FUNCTION bessj1_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::bessj1_s END FUNCTION bessj1_s FUNCTION bessj1_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessj1_v END FUNCTION bessj1_v END INTERFACE INTERFACE bessjy SUBROUTINE bessjy_s(x,xnu,rj,ry,rjp,ryp)USE nrtype REAL(SP),INTENT(IN)::x,xnu

REAL(SP),INTENT(OUT)::rj,ry,rjp,ryp

END SUBROUTINE bessjy_s

SUBROUTINE bessjy_v(x,xnu,rj,ry,rjp,ryp)

USE nrtype

REAL(SP),INTENT(IN)::xnu

REAL(SP),DIMENSION(:),INTENT(IN)::x

REAL(SP),DIMENSION(:),INTENT(OUT)::rj,rjp,ry,ryp

1388Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).END SUBROUTINE bessjy_v END INTERFACE INTERFACE bessk FUNCTION bessk_s(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),INTENT(IN)::x REAL(SP)::bessk_s END FUNCTION bessk_s FUNCTION bessk_v(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessk_v END FUNCTION bessk_v END INTERFACE INTERFACE bessk0FUNCTION bessk0_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::bessk0_s END FUNCTION bessk0_s FUNCTION bessk0_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessk0_v END FUNCTION bessk0_v END INTERFACE INTERFACE bessk1FUNCTION bessk1_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::bessk1_s END FUNCTION bessk1_s FUNCTION bessk1_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessk1_v END FUNCTION bessk1_v END INTERFACE INTERFACE bessy FUNCTION bessy_s(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),INTENT(IN)::x REAL(SP)::bessy_s END FUNCTION bessy_s FUNCTION bessy_v(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessy_v END FUNCTION bessy_v END INTERFACE

INTERFACE bessy0

FUNCTION bessy0_s(x)

USE nrtype

REAL(SP),INTENT(IN)::x

REAL(SP)::bessy0_s

END FUNCTION bessy0_s

FUNCTION bessy0_v(x)

USE nrtype

Appendix C2.Alphabetical Listing of Explicit Interfaces1389 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessy0_v END FUNCTION bessy0_v END INTERFACE INTERFACE bessy1FUNCTION bessy1_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::bessy1_s END FUNCTION bessy1_s FUNCTION bessy1_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::bessy1_v END FUNCTION bessy1_v END INTERFACE INTERFACE beta FUNCTION beta_s(z,w)USE nrtype REAL(SP),INTENT(IN)::z,w REAL(SP)::beta_s END FUNCTION beta_s FUNCTION beta_v(z,w)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::z,w REAL(SP),DIMENSION(size(z))::beta_v END FUNCTION beta_v END INTERFACE INTERFACE betacf FUNCTION betacf_s(a,b,x)USE nrtype REAL(SP),INTENT(IN)::a,b,x REAL(SP)::betacf_s END FUNCTION betacf_s FUNCTION betacf_v(a,b,x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::a,b,x REAL(SP),DIMENSION(size(x))::betacf_v END FUNCTION betacf_v END INTERFACE INTERFACE betai FUNCTION betai_s(a,b,x)USE nrtype REAL(SP),INTENT(IN)::a,b,x REAL(SP)::betai_s END FUNCTION betai_s FUNCTION betai_v(a,b,x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::a,b,x REAL(SP),DIMENSION(size(a))::betai_v END FUNCTION betai_v END INTERFACE INTERFACE bico FUNCTION bico_s(n,k)

USE nrtype

INTEGER(I4B),INTENT(IN)::n,k

REAL(SP)::bico_s

END FUNCTION bico_s

FUNCTION bico_v(n,k)

USE nrtype

INTEGER(I4B),DIMENSION(:),INTENT(IN)::n,k

REAL(SP),DIMENSION(size(n))::bico_v

1390Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).END FUNCTION bico_v END INTERFACE INTERFACE FUNCTION bnldev(pp,n)USE nrtype REAL(SP),INTENT(IN)::pp INTEGER(I4B),INTENT(IN)::n REAL(SP)::bnldev END FUNCTION bnldev END INTERFACE INTERFACE FUNCTION brent(ax,bx,cx,func,tol,xmin)USE nrtype REAL(SP),INTENT(IN)::ax,bx,cx,tol REAL(SP),INTENT(OUT)::xmin REAL(SP)::brent INTERFACE FUNCTION func(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::func END FUNCTION func END INTERFACE END FUNCTION brent END INTERFACE INTERFACE SUBROUTINE broydn(x,check)USE nrtype REAL(SP),DIMENSION(:),INTENT(INOUT)::x LOGICAL(LGT),INTENT(OUT)::check END SUBROUTINE broydn END INTERFACE INTERFACE SUBROUTINE bsstep(y,dydx,x,htry,eps,yscal,hdid,hnext,derivs)USE nrtype REAL(SP),DIMENSION(:),INTENT(INOUT)::y REAL(SP),DIMENSION(:),INTENT(IN)::dydx,yscal REAL(SP),INTENT(INOUT)::x REAL(SP),INTENT(IN)::htry,eps REAL(SP),INTENT(OUT)::hdid,hnext INTERFACE SUBROUTINE derivs(x,y,dydx)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP),DIMENSION(:),INTENT(IN)::y REAL(SP),DIMENSION(:),INTENT(OUT)::dydx END SUBROUTINE derivs END INTERFACE END SUBROUTINE bsstep END INTERFACE INTERFACE SUBROUTINE caldat(julian,mm,id,iyyy)USE nrtype INTEGER(I4B),INTENT(IN)::julian INTEGER(I4B),INTENT(OUT)::mm,id,iyyy END SUBROUTINE caldat

END INTERFACE

INTERFACE

FUNCTION chder(a,b,c)

USE nrtype

REAL(SP),INTENT(IN)::a,b

REAL(SP),DIMENSION(:),INTENT(IN)::c

REAL(SP),DIMENSION(size(c))::chder

END FUNCTION chder

Appendix C2.Alphabetical Listing of Explicit Interfaces1391 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).END INTERFACE INTERFACE chebev FUNCTION chebev_s(a,b,c,x)USE nrtype REAL(SP),INTENT(IN)::a,b,x REAL(SP),DIMENSION(:),INTENT(IN)::c REAL(SP)::chebev_s END FUNCTION chebev_s FUNCTION chebev_v(a,b,c,x)USE nrtype REAL(SP),INTENT(IN)::a,b REAL(SP),DIMENSION(:),INTENT(IN)::c,x REAL(SP),DIMENSION(size(x))::chebev_v END FUNCTION chebev_v END INTERFACE INTERFACE FUNCTION chebft(a,b,n,func)USE nrtype REAL(SP),INTENT(IN)::a,b INTEGER(I4B),INTENT(IN)::n REAL(SP),DIMENSION(n)::chebft INTERFACE FUNCTION func(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::func END FUNCTION func END INTERFACE END FUNCTION chebft END INTERFACE INTERFACE FUNCTION chebpc(c)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::c REAL(SP),DIMENSION(size(c))::chebpc END FUNCTION chebpc END INTERFACE INTERFACE FUNCTION chint(a,b,c)USE nrtype REAL(SP),INTENT(IN)::a,b REAL(SP),DIMENSION(:),INTENT(IN)::c REAL(SP),DIMENSION(size(c))::chint END FUNCTION chint END INTERFACE INTERFACE SUBROUTINE choldc(a,p)USE nrtype REAL(SP),DIMENSION(:,:),INTENT(INOUT)::a REAL(SP),DIMENSION(:),INTENT(OUT)::p END SUBROUTINE choldc END INTERFACE INTERFACE SUBROUTINE cholsl(a,p,b,x)USE nrtype REAL(SP),DIMENSION(:,:),INTENT(IN)::a

REAL(SP),DIMENSION(:),INTENT(IN)::p,b

REAL(SP),DIMENSION(:),INTENT(INOUT)::x

END SUBROUTINE cholsl

END INTERFACE

INTERFACE

SUBROUTINE chsone(bins,ebins,knstrn,df,chsq,prob)

USE nrtype

1392Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).INTEGER(I4B),INTENT(IN)::knstrn REAL(SP),INTENT(OUT)::df,chsq,prob REAL(SP),DIMENSION(:),INTENT(IN)::bins,ebins END SUBROUTINE chsone END INTERFACE INTERFACE SUBROUTINE chstwo(bins1,bins2,knstrn,df,chsq,prob)USE nrtype INTEGER(I4B),INTENT(IN)::knstrn REAL(SP),INTENT(OUT)::df,chsq,prob REAL(SP),DIMENSION(:),INTENT(IN)::bins1,bins2END SUBROUTINE chstwo END INTERFACE INTERFACE SUBROUTINE cisi(x,ci,si)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP),INTENT(OUT)::ci,si END SUBROUTINE cisi END INTERFACE INTERFACE SUBROUTINE cntab1(nn,chisq,df,prob,cramrv,ccc)USE nrtype INTEGER(I4B),DIMENSION(:,:),INTENT(IN)::nn REAL(SP),INTENT(OUT)::chisq,df,prob,cramrv,ccc END SUBROUTINE cntab1END INTERFACE INTERFACE SUBROUTINE cntab2(nn,h,hx,hy,hygx,hxgy,uygx,uxgy,uxy)USE nrtype INTEGER(I4B),DIMENSION(:,:),INTENT(IN)::nn REAL(SP),INTENT(OUT)::h,hx,hy,hygx,hxgy,uygx,uxgy,uxy END SUBROUTINE cntab2END INTERFACE INTERFACE FUNCTION convlv(data,respns,isign)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::data REAL(SP),DIMENSION(:),INTENT(IN)::respns INTEGER(I4B),INTENT(IN)::isign REAL(SP),DIMENSION(size(data))::convlv END FUNCTION convlv END INTERFACE INTERFACE FUNCTION correl(data1,data2)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::data1,data2REAL(SP),DIMENSION(size(data1))::correl END FUNCTION correl END INTERFACE INTERFACE SUBROUTINE cosft1(y)USE nrtype REAL(SP),DIMENSION(:),INTENT(INOUT)::y END SUBROUTINE cosft1END INTERFACE

INTERFACE

SUBROUTINE cosft2(y,isign)

USE nrtype

REAL(SP),DIMENSION(:),INTENT(INOUT)::y

INTEGER(I4B),INTENT(IN)::isign

END SUBROUTINE cosft2

END INTERFACE

INTERFACE

Appendix C2.Alphabetical Listing of Explicit Interfaces1393 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).SUBROUTINE covsrt(covar,maska)USE nrtype REAL(SP),DIMENSION(:,:),INTENT(INOUT)::covar LOGICAL(LGT),DIMENSION(:),INTENT(IN)::maska END SUBROUTINE covsrt END INTERFACE INTERFACE SUBROUTINE cyclic(a,b,c,alpha,beta,r,x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::a,b,c,r REAL(SP),INTENT(IN)::alpha,beta REAL(SP),DIMENSION(:),INTENT(OUT)::x END SUBROUTINE cyclic END INTERFACE INTERFACE SUBROUTINE daub4(a,isign)USE nrtype REAL(SP),DIMENSION(:),INTENT(INOUT)::a INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE daub4END INTERFACE INTERFACE dawson FUNCTION dawson_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::dawson_s END FUNCTION dawson_s FUNCTION dawson_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::dawson_v END FUNCTION dawson_v END INTERFACE INTERFACE FUNCTION dbrent(ax,bx,cx,func,dbrent_dfunc,tol,xmin)USE nrtype REAL(SP),INTENT(IN)::ax,bx,cx,tol REAL(SP),INTENT(OUT)::xmin REAL(SP)::dbrent INTERFACE FUNCTION func(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::func END FUNCTION func FUNCTION dbrent_dfunc(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::dbrent_dfunc END FUNCTION dbrent_dfunc END INTERFACE END FUNCTION dbrent END INTERFACE INTERFACE SUBROUTINE ddpoly(c,x,pd)

USE nrtype

REAL(SP),INTENT(IN)::x

REAL(SP),DIMENSION(:),INTENT(IN)::c

REAL(SP),DIMENSION(:),INTENT(OUT)::pd

END SUBROUTINE ddpoly

END INTERFACE

INTERFACE

FUNCTION decchk(string,ch)

1394Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).USE nrtype CHARACTER(1),DIMENSION(:),INTENT(IN)::string CHARACTER(1),INTENT(OUT)::ch LOGICAL(LGT)::decchk END FUNCTION decchk END INTERFACE INTERFACE SUBROUTINE dfpmin(p,gtol,iter,fret,func,dfunc)USE nrtype INTEGER(I4B),INTENT(OUT)::iter REAL(SP),INTENT(IN)::gtol REAL(SP),INTENT(OUT)::fret REAL(SP),DIMENSION(:),INTENT(INOUT)::p INTERFACE FUNCTION func(p)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::p REAL(SP)::func END FUNCTION func FUNCTION dfunc(p)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::p REAL(SP),DIMENSION(size(p))::dfunc END FUNCTION dfunc END INTERFACE END SUBROUTINE dfpmin END INTERFACE INTERFACE FUNCTION dfridr(func,x,h,err)USE nrtype REAL(SP),INTENT(IN)::x,h REAL(SP),INTENT(OUT)::err REAL(SP)::dfridr INTERFACE FUNCTION func(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::func END FUNCTION func END INTERFACE END FUNCTION dfridr END INTERFACE INTERFACE SUBROUTINE dftcor(w,delta,a,b,endpts,corre,corim,corfac)USE nrtype REAL(SP),INTENT(IN)::w,delta,a,b REAL(SP),INTENT(OUT)::corre,corim,corfac REAL(SP),DIMENSION(:),INTENT(IN)::endpts END SUBROUTINE dftcor END INTERFACE INTERFACE SUBROUTINE dftint(func,a,b,w,cosint,sinint)USE nrtype REAL(SP),INTENT(IN)::a,b,w REAL(SP),INTENT(OUT)::cosint,sinint INTERFACE

FUNCTION func(x)

USE nrtype

REAL(SP),DIMENSION(:),INTENT(IN)::x

REAL(SP),DIMENSION(size(x))::func

END FUNCTION func

END INTERFACE

END SUBROUTINE dftint

Appendix C2.Alphabetical Listing of Explicit Interfaces1395 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).END INTERFACE INTERFACE SUBROUTINE difeq(k,k1,k2,jsf,is1,isf,indexv,s,y)USE nrtype INTEGER(I4B),INTENT(IN)::is1,isf,jsf,k,k1,k2INTEGER(I4B),DIMENSION(:),INTENT(IN)::indexv REAL(SP),DIMENSION(:,:),INTENT(OUT)::s REAL(SP),DIMENSION(:,:),INTENT(IN)::y END SUBROUTINE difeq END INTERFACE INTERFACE FUNCTION eclass(lista,listb,n)USE nrtype INTEGER(I4B),DIMENSION(:),INTENT(IN)::lista,listb INTEGER(I4B),INTENT(IN)::n INTEGER(I4B),DIMENSION(n)::eclass END FUNCTION eclass END INTERFACE INTERFACE FUNCTION eclazz(equiv,n)USE nrtype INTERFACE FUNCTION equiv(i,j)USE nrtype LOGICAL(LGT)::equiv INTEGER(I4B),INTENT(IN)::i,j END FUNCTION equiv END INTERFACE INTEGER(I4B),INTENT(IN)::n INTEGER(I4B),DIMENSION(n)::eclazz END FUNCTION eclazz END INTERFACE INTERFACE FUNCTION ei(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::ei END FUNCTION ei END INTERFACE INTERFACE SUBROUTINE eigsrt(d,v)USE nrtype REAL(SP),DIMENSION(:),INTENT(INOUT)::d REAL(SP),DIMENSION(:,:),INTENT(INOUT)::v END SUBROUTINE eigsrt END INTERFACE INTERFACE elle FUNCTION elle_s(phi,ak)USE nrtype REAL(SP),INTENT(IN)::phi,ak REAL(SP)::elle_s END FUNCTION elle_s FUNCTION elle_v(phi,ak)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::phi,ak REAL(SP),DIMENSION(size(phi))::elle_v

END FUNCTION elle_v

END INTERFACE

INTERFACE ellf

FUNCTION ellf_s(phi,ak)

USE nrtype

REAL(SP),INTENT(IN)::phi,ak

REAL(SP)::ellf_s

1396Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).END FUNCTION ellf_s FUNCTION ellf_v(phi,ak)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::phi,ak REAL(SP),DIMENSION(size(phi))::ellf_v END FUNCTION ellf_v END INTERFACE INTERFACE ellpi FUNCTION ellpi_s(phi,en,ak)USE nrtype REAL(SP),INTENT(IN)::phi,en,ak REAL(SP)::ellpi_s END FUNCTION ellpi_s FUNCTION ellpi_v(phi,en,ak)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::phi,en,ak REAL(SP),DIMENSION(size(phi))::ellpi_v END FUNCTION ellpi_v END INTERFACE INTERFACE SUBROUTINE elmhes(a)USE nrtype REAL(SP),DIMENSION(:,:),INTENT(INOUT)::a END SUBROUTINE elmhes END INTERFACE INTERFACE erf FUNCTION erf_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::erf_s END FUNCTION erf_s FUNCTION erf_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::erf_v END FUNCTION erf_v END INTERFACE INTERFACE erfc FUNCTION erfc_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::erfc_s END FUNCTION erfc_s FUNCTION erfc_v(x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x REAL(SP),DIMENSION(size(x))::erfc_v END FUNCTION erfc_v END INTERFACE INTERFACE erfcc FUNCTION erfcc_s(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::erfcc_s

END FUNCTION erfcc_s

FUNCTION erfcc_v(x)

USE nrtype

REAL(SP),DIMENSION(:),INTENT(IN)::x

REAL(SP),DIMENSION(size(x))::erfcc_v

END FUNCTION erfcc_v

END INTERFACE

INTERFACE

Appendix C2.Alphabetical Listing of Explicit Interfaces1397 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).SUBROUTINE eulsum(sum,term,jterm)USE nrtype REAL(SP),INTENT(INOUT)::sum REAL(SP),INTENT(IN)::term INTEGER(I4B),INTENT(IN)::jterm END SUBROUTINE eulsum END INTERFACE INTERFACE FUNCTION evlmem(fdt,d,xms)USE nrtype REAL(SP),INTENT(IN)::fdt,xms REAL(SP),DIMENSION(:),INTENT(IN)::d REAL(SP)::evlmem END FUNCTION evlmem END INTERFACE INTERFACE expdev SUBROUTINE expdev_s(harvest)USE nrtype REAL(SP),INTENT(OUT)::harvest END SUBROUTINE expdev_s SUBROUTINE expdev_v(harvest)USE nrtype REAL(SP),DIMENSION(:),INTENT(OUT)::harvest END SUBROUTINE expdev_v END INTERFACE INTERFACE FUNCTION expint(n,x)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP),INTENT(IN)::x REAL(SP)::expint END FUNCTION expint END INTERFACE INTERFACE factln FUNCTION factln_s(n)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP)::factln_s END FUNCTION factln_s FUNCTION factln_v(n)USE nrtype INTEGER(I4B),DIMENSION(:),INTENT(IN)::n REAL(SP),DIMENSION(size(n))::factln_v END FUNCTION factln_v END INTERFACE INTERFACE factrl FUNCTION factrl_s(n)USE nrtype INTEGER(I4B),INTENT(IN)::n REAL(SP)::factrl_s END FUNCTION factrl_s FUNCTION factrl_v(n)USE nrtype INTEGER(I4B),DIMENSION(:),INTENT(IN)::n REAL(SP),DIMENSION(size(n))::factrl_v

END FUNCTION factrl_v

END INTERFACE

INTERFACE

SUBROUTINE fasper(x,y,ofac,hifac,px,py,jmax,prob)

USE nrtype

REAL(SP),DIMENSION(:),INTENT(IN)::x,y

REAL(SP),INTENT(IN)::ofac,hifac

INTEGER(I4B),INTENT(OUT)::jmax

1398Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).REAL(SP),INTENT(OUT)::prob REAL(SP),DIMENSION(:),POINTER::px,py END SUBROUTINE fasper END INTERFACE INTERFACE SUBROUTINE fdjac(x,fvec,df)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::fvec REAL(SP),DIMENSION(:),INTENT(INOUT)::x REAL(SP),DIMENSION(:,:),INTENT(OUT)::df END SUBROUTINE fdjac END INTERFACE INTERFACE SUBROUTINE fgauss(x,a,y,dyda)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x,a REAL(SP),DIMENSION(:),INTENT(OUT)::y REAL(SP),DIMENSION(:,:),INTENT(OUT)::dyda END SUBROUTINE fgauss END INTERFACE INTERFACE SUBROUTINE fit(x,y,a,b,siga,sigb,chi2,q,sig)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x,y REAL(SP),INTENT(OUT)::a,b,siga,sigb,chi2,q REAL(SP),DIMENSION(:),OPTIONAL,INTENT(IN)::sig END SUBROUTINE fit END INTERFACE INTERFACE SUBROUTINE fitexy(x,y,sigx,sigy,a,b,siga,sigb,chi2,q)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::x,y,sigx,sigy REAL(SP),INTENT(OUT)::a,b,siga,sigb,chi2,q END SUBROUTINE fitexy END INTERFACE INTERFACE SUBROUTINE fixrts(d)USE nrtype REAL(SP),DIMENSION(:),INTENT(INOUT)::d END SUBROUTINE fixrts END INTERFACE INTERFACE FUNCTION fleg(x,n)USE nrtype REAL(SP),INTENT(IN)::x INTEGER(I4B),INTENT(IN)::n REAL(SP),DIMENSION(n)::fleg END FUNCTION fleg END INTERFACE INTERFACE SUBROUTINE flmoon(n,nph,jd,frac)USE nrtype INTEGER(I4B),INTENT(IN)::n,nph INTEGER(I4B),INTENT(OUT)::jd REAL(SP),INTENT(OUT)::frac END SUBROUTINE flmoon

END INTERFACE

INTERFACE four1

SUBROUTINE four1_dp(data,isign)

USE nrtype

COMPLEX(DPC),DIMENSION(:),INTENT(INOUT)::data

INTEGER(I4B),INTENT(IN)::isign

END SUBROUTINE four1_dp

Appendix C2.Alphabetical Listing of Explicit Interfaces1399 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).SUBROUTINE four1_sp(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE four1_sp END INTERFACE INTERFACE SUBROUTINE four1_alt(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE four1_alt END INTERFACE INTERFACE SUBROUTINE four1_gather(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE four1_gather END INTERFACE INTERFACE SUBROUTINE four2(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:,:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE four2END INTERFACE INTERFACE SUBROUTINE four2_alt(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:,:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE four2_alt END INTERFACE INTERFACE SUBROUTINE four3(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:,:,:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE four3END INTERFACE INTERFACE SUBROUTINE four3_alt(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:,:,:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE four3_alt END INTERFACE INTERFACE SUBROUTINE fourcol(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:,:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE fourcol END INTERFACE INTERFACE

SUBROUTINE fourcol_3d(data,isign)

USE nrtype

COMPLEX(SPC),DIMENSION(:,:,:),INTENT(INOUT)::data

INTEGER(I4B),INTENT(IN)::isign

END SUBROUTINE fourcol_3d

END INTERFACE

INTERFACE

SUBROUTINE fourn_gather(data,nn,isign)

1400Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).USE nrtype COMPLEX(SPC),DIMENSION(:),INTENT(INOUT)::data INTEGER(I4B),DIMENSION(:),INTENT(IN)::nn INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE fourn_gather END INTERFACE INTERFACE fourrow SUBROUTINE fourrow_dp(data,isign)USE nrtype COMPLEX(DPC),DIMENSION(:,:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE fourrow_dp SUBROUTINE fourrow_sp(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:,:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE fourrow_sp END INTERFACE INTERFACE SUBROUTINE fourrow_3d(data,isign)USE nrtype COMPLEX(SPC),DIMENSION(:,:,:),INTENT(INOUT)::data INTEGER(I4B),INTENT(IN)::isign END SUBROUTINE fourrow_3d END INTERFACE INTERFACE FUNCTION fpoly(x,n)USE nrtype REAL(SP),INTENT(IN)::x INTEGER(I4B),INTENT(IN)::n REAL(SP),DIMENSION(n)::fpoly END FUNCTION fpoly END INTERFACE INTERFACE SUBROUTINE fred2(a,b,t,f,w,g,ak)USE nrtype REAL(SP),INTENT(IN)::a,b REAL(SP),DIMENSION(:),INTENT(OUT)::t,f,w INTERFACE FUNCTION g(t)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::t REAL(SP),DIMENSION(size(t))::g END FUNCTION g FUNCTION ak(t,s)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::t,s REAL(SP),DIMENSION(size(t),size(s))::ak END FUNCTION ak END INTERFACE END SUBROUTINE fred2END INTERFACE INTERFACE FUNCTION fredin(x,a,b,t,f,w,g,ak)USE nrtype

REAL(SP),INTENT(IN)::a,b

REAL(SP),DIMENSION(:),INTENT(IN)::x,t,f,w

REAL(SP),DIMENSION(size(x))::fredin

INTERFACE

FUNCTION g(t)

USE nrtype

REAL(SP),DIMENSION(:),INTENT(IN)::t

REAL(SP),DIMENSION(size(t))::g

Appendix C2.Alphabetical Listing of Explicit Interfaces1401 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).END FUNCTION g FUNCTION ak(t,s)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::t,s REAL(SP),DIMENSION(size(t),size(s))::ak END FUNCTION ak END INTERFACE END FUNCTION fredin END INTERFACE INTERFACE SUBROUTINE frenel(x,s,c)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP),INTENT(OUT)::s,c END SUBROUTINE frenel END INTERFACE INTERFACE SUBROUTINE frprmn(p,ftol,iter,fret)USE nrtype INTEGER(I4B),INTENT(OUT)::iter REAL(SP),INTENT(IN)::ftol REAL(SP),INTENT(OUT)::fret REAL(SP),DIMENSION(:),INTENT(INOUT)::p END SUBROUTINE frprmn END INTERFACE INTERFACE SUBROUTINE ftest(data1,data2,f,prob)USE nrtype REAL(SP),INTENT(OUT)::f,prob REAL(SP),DIMENSION(:),INTENT(IN)::data1,data2END SUBROUTINE ftest END INTERFACE INTERFACE FUNCTION gamdev(ia)USE nrtype INTEGER(I4B),INTENT(IN)::ia REAL(SP)::gamdev END FUNCTION gamdev END INTERFACE INTERFACE gammln FUNCTION gammln_s(xx)USE nrtype REAL(SP),INTENT(IN)::xx REAL(SP)::gammln_s END FUNCTION gammln_s FUNCTION gammln_v(xx)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::xx REAL(SP),DIMENSION(size(xx))::gammln_v END FUNCTION gammln_v END INTERFACE INTERFACE gammp FUNCTION gammp_s(a,x)USE nrtype REAL(SP),INTENT(IN)::a,x

REAL(SP)::gammp_s

END FUNCTION gammp_s

FUNCTION gammp_v(a,x)

USE nrtype

REAL(SP),DIMENSION(:),INTENT(IN)::a,x

REAL(SP),DIMENSION(size(a))::gammp_v

END FUNCTION gammp_v

END INTERFACE

1402Appendix C2.Alphabetical Listing of Explicit Interfaces Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).INTERFACE gammq FUNCTION gammq_s(a,x)USE nrtype REAL(SP),INTENT(IN)::a,x REAL(SP)::gammq_s END FUNCTION gammq_s FUNCTION gammq_v(a,x)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::a,x REAL(SP),DIMENSION(size(a))::gammq_v END FUNCTION gammq_v END INTERFACE INTERFACE gasdev SUBROUTINE gasdev_s(harvest)USE nrtype REAL(SP),INTENT(OUT)::harvest END SUBROUTINE gasdev_s SUBROUTINE gasdev_v(harvest)USE nrtype REAL(SP),DIMENSION(:),INTENT(OUT)::harvest END SUBROUTINE gasdev_v END INTERFACE INTERFACE SUBROUTINE gaucof(a,b,amu0,x,w)USE nrtype REAL(SP),INTENT(IN)::amu0REAL(SP),DIMENSION(:),INTENT(INOUT)::a,b REAL(SP),DIMENSION(:),INTENT(OUT)::x,w END SUBROUTINE gaucof END INTERFACE INTERFACE SUBROUTINE gauher(x,w)USE nrtype REAL(SP),DIMENSION(:),INTENT(OUT)::x,w END SUBROUTINE gauher END INTERFACE INTERFACE SUBROUTINE gaujac(x,w,alf,bet)USE nrtype REAL(SP),INTENT(IN)::alf,bet REAL(SP),DIMENSION(:),INTENT(OUT)::x,w END SUBROUTINE gaujac END INTERFACE INTERFACE SUBROUTINE gaulag(x,w,alf)USE nrtype REAL(SP),INTENT(IN)::alf REAL(SP),DIMENSION(:),INTENT(OUT)::x,w END SUBROUTINE gaulag END INTERFACE INTERFACE SUBROUTINE gauleg(x1,x2,x,w)USE nrtype REAL(SP),INTENT(IN)::x1,x2REAL(SP),DIMENSION(:),INTENT(OUT)::x,w

END SUBROUTINE gauleg

END INTERFACE

INTERFACE

SUBROUTINE gaussj(a,b)

USE nrtype

REAL(SP),DIMENSION(:,:),INTENT(INOUT)::a,b

END SUBROUTINE gaussj

END INTERFACE

Appendix C2.Alphabetical Listing of Explicit Interfaces1403 Sample page from NUMERICAL RECIPES IN FORTRAN 90: THE Art of PARALLEL Scientific Computing (ISBN 0-521-57439-0)Copyright (C) 1986-1996 by Cambridge University Press.Programs Copyright (C) 1986-1996 by Numerical Recipes Software. Permission is granted for internet users to make one paper copy for their own personal use. Further reproduction, or any copying of machine-readable files (including this one) to any server computer, is strictly prohibited. To order Numerical Recipes books or CDROMs, visit website https://www.doczj.com/doc/a516262195.html, or call 1-800-872-7423 (North America only),or send email to directcustserv@https://www.doczj.com/doc/a516262195.html, (outside North America).INTERFACE gcf FUNCTION gcf_s(a,x,gln)USE nrtype REAL(SP),INTENT(IN)::a,x REAL(SP),OPTIONAL,INTENT(OUT)::gln REAL(SP)::gcf_s END FUNCTION gcf_s FUNCTION gcf_v(a,x,gln)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::a,x REAL(SP),DIMENSION(:),OPTIONAL,INTENT(OUT)::gln REAL(SP),DIMENSION(size(a))::gcf_v END FUNCTION gcf_v END INTERFACE INTERFACE FUNCTION golden(ax,bx,cx,func,tol,xmin)USE nrtype REAL(SP),INTENT(IN)::ax,bx,cx,tol REAL(SP),INTENT(OUT)::xmin REAL(SP)::golden INTERFACE FUNCTION func(x)USE nrtype REAL(SP),INTENT(IN)::x REAL(SP)::func END FUNCTION func END INTERFACE END FUNCTION golden END INTERFACE INTERFACE gser FUNCTION gser_s(a,x,gln)USE nrtype REAL(SP),INTENT(IN)::a,x REAL(SP),OPTIONAL,INTENT(OUT)::gln REAL(SP)::gser_s END FUNCTION gser_s FUNCTION gser_v(a,x,gln)USE nrtype REAL(SP),DIMENSION(:),INTENT(IN)::a,x REAL(SP),DIMENSION(:),OPTIONAL,INTENT(OUT)::gln REAL(SP),DIMENSION(size(a))::gser_v END FUNCTION gser_v END INTERFACE INTERFACE SUBROUTINE hqr(a,wr,wi)USE nrtype REAL(SP),DIMENSION(:),INTENT(OUT)::wr,wi REAL(SP),DIMENSION(:,:),INTENT(INOUT)::a END SUBROUTINE hqr END INTERFACE INTERFACE SUBROUTINE hunt(xx,x,jlo)USE nrtype INTEGER(I4B),INTENT(INOUT)::jlo REAL(SP),INTENT(IN)::x

REAL(SP),DIMENSION(:),INTENT(IN)::xx

END SUBROUTINE hunt

END INTERFACE

INTERFACE

SUBROUTINE hypdrv(s,ry,rdyds)

USE nrtype

REAL(SP),INTENT(IN)::s

REAL(SP),DIMENSION(:),INTENT(IN)::ry

数值计算方法课程设计

重庆邮电大学 数学与应用数学 专业 《数值计算方法》课程设计 姓名: 李金徽 王莹 刘姝楠 班级: 1131001 1131002 1131002 学号: 2010213542 2010213570 2010213571 设计时间: 2012-6-4 指导教师: 朱伟

一、课程设计目的 在科学计算与工程设计中,我们常会遇到求解线性方程组的问题,对于系数矩阵为低阶稠密矩阵的线性方程组,可以用直接法进行消元,而对于系数矩阵为大型稀疏矩阵的情况,直接法就显得比较繁琐,而迭代法比较适用。比较常用的迭代法有Jacobi 迭代与Gauss - seidel 迭代。本文基于两种方法设计算法,并比较他们的优劣。 二、课程设计内容 给出Jacobi 迭代法和Gauss-Seidel 迭代法求解线性方程组的算法思想和MATLAB 程序实现,并对比分析这两种算法的优劣。 三、问题的分析(含涉及的理论知识、算法等) Jacobi 迭代法 方程组迭代法的基本思想和求根的迭代法思想类似,即对于线性 方程组Ax = b( 其中n n n R b R R A ∈?∈,),即方程组 )1(2211222221211 1212111?? ???? ?=+?++??=+?++=+?++n n nn n n n n n n b x a x a x a b x a x a x a b x a x a x a 将系数矩阵A 写为 )2(000000 21122 12122 11U L D a a a a a a a a a A n n n n nn --≡??? ?? ? ? ??---- ??????? ??----??????? ??= 若选取D M =,则U L A M N +=-=,方程组)1(转化为等价方程组 b x U L Dx ++=)(

数值计算方法学习心得

数值计算方法学习心得 ------一个代码的方法是很重要,一个算法的思想也很重要,但 在我看来,更重要的是解决问题的方法,就像爱因斯坦说的内容比 思维本身更重要。 我上去讲的那次其实做了挺充分的准备,程序的运行,pdf文档,算法公式的推导,程序伪代码,不过有一点缺陷的地方,很多细节 没有讲的很清楚吧,下来之后也是更清楚了这个问题。 然后一学期下来,总的来说,看其他同学的分享,我也学习到 许多东西,并非只是代码的方法,更多的是章胜同学的口才,攀忠 的排版,小冯的深入挖掘…都是对我而言比算法更加值得珍惜的东西,又骄傲地回想一下,曾同为一个项目组的我们也更加感到做项 目对自己发展的巨大帮助了。 同时从这些次的实验中我发现以前学到的很多知识都非常有用。 比如说,以前做项目的时候,项目导师一直要求对于要上传的 文件尽量用pdf格式,不管是ppt还是文档,这便算是对产权的一种 保护。 再比如代码分享,最基础的要求便是——其他人拿到你的代码 也能运行出来,其次是代码分享的规范性,像我们可以用轻量级Ubuntu Pastebin,以前做过一小段时间acm,集训队里对于代码的分享都是推荐用这个,像数值计算实验我觉得用这个也差不多了,其 次项目级代码还是推荐github(被微软收购了),它的又是可能更 多在于个人代码平台的搭建,当然像readme文档及必要的一些数据 集放在上面都更方便一些。

然后在实验中,发现debug能力的重要性,对于代码错误点的 正确分析,以及一些与他人交流的“正规”途径,讨论算法可能出 错的地方以及要注意的细节等,比如acm比赛都是以三人为一小组,讨论过后,讲了一遍会发现自己对算法理解更加深刻。 然后学习算法,做项目做算法一般的正常流程是看论文,尽量 看英文文献,一般就是第一手资料,然后根据论文对算法的描述, 就是如同课上的流程一样,对算法进一步理解,然后进行复现,最 后就是尝试自己改进。比如知网查询牛顿法相关论文,会找到大量 可以参考的文献。 最后的最后,想说一下,计算机专业的同学看这个数值分析, 不一定行云流水,但肯定不至于看不懂写不出来,所以我们还是要 提高自己的核心竞争力,就是利用我们的优势,对于这种算法方面 的编程,至少比他们用的更加熟练,至少面对一个问题,我们能思 考出对应问题的最佳算法是哪一个更合适解决问题。 附记: 对课程的一些小建议: 1. debug的能力不容忽视,比如给一个关于代码实现已知错误的代码给同学们,让同学们自己思考一下,然后分享各自的debug方法,一步一步的去修改代码,最后集全班的力量完成代码的debug,这往往更能提升同学们的代码能力。 2. 课堂上的效率其实是有点低的,可能会给学生带来一些负反馈,降低学习热情。 3. 总的来说还是从这门课程中学到许多东西。 数值分析学习心得体会

电大网院统考大学英语B模拟试题9套

电大网院统考大学英语B模拟试题9套

电大网院统考大学英语B模拟试题9套 大学英语B9 一、交际英语 1、-Our class won the girls' relay race in our school's sports meeting.- _____________________ A:That's right. B:Thanks a lot. C:Congratulations! D:It's a pleasure. 答案: C 2、- I heard you got a full mark in the English exam. Congratulations!- _________ A:Oh, no, no. B:No, I didn't do very well. C:Thanks. D:Don't praise me. 答案: C 3、- Could I ask you a rather personal question? - Of course, _____. A:good idea B:that's right C:never mind D:go ahead 答案: D 4、- Marilyn, I'm afraid I have to be leaving now.- ____________ A:That sounds wonderful. B:Oh, so early? C:Not at all. D:Good luck! 答案: B 5、- Excuse me, how much is the jacket?- It's 499 Yuan. A:Oh, no. That's OK! B:How do you like it? C:Which do you prefer? D:Would you like to try it on? 答案: D

现代数值计算方法习题答

现代数值计算方法习题答案 习 题 一 1、解:根据绝对误差限不超过末位数的半个单位,相对误差限为绝对误差限除以 有效数字本身,有效数字的位数根据有效数字的定义来求.因此 49×10 -2 :E = 0.005; r E = 0.0102; 2位有效数字. 0.0490 :E = 0.00005;r E = 0.00102; 3位有效数字. 490.00 :E = 0.005; r E = 0.0000102;5位有效数字. 2、解: 7 22 = 3.1428 …… , π = 3.1415 …… , 取它们的相同部分3.14,故有3位有效数字. E = 3.1428 - 3.1415 = 0.0013 ;r E = 14 .3E = 14 .30013.0 = 0.00041. 3、解:101的近似值的首位非0数字1α = 1,因此有 |)(*x E r |) 1(10 1 21--??=n < = 2 1× 10 -4 , 解之得n > = 5,所以 n = 5 . 4、证:) ()(1)()(1)(* 1 1* * 1 1 * * x x x n x E x n x E n n n -= ≈ -- )(11)()(1) ()(* * * * * 1 1 ** * * x E n x x x n x x x x n x x E x E r n n n n n r = -= -≈ = - 5、解:(1)因为=20 4.4721…… , 又=)(*x E |*x x -| = |47.420-| = 0.0021 < 0.01, 所以 =*x 4.47. (2)20的近似值的首位非0数字1α = 4,因此有 |)(*x E r |) 1(10 4 21--??= n < = 0.01 , 解之得n > = 3 .所以,=*x 4.47. 6、解:设正方形的边长为x ,则其面积为2x y =,由题设知x 的近似值为*x = 10 cm . 记*y 为y 的近似值,则

数值计算方法课程报告

课程报告 课程名称______《数值计算》 __ 学生学院_____机电工程学院___ 专业班级_____微电子(1)班____ 学号________ 学生姓名_______________ 指导教师_____ ________ XXXX年XX月XX日

姓 名: 线 学 号 : 订 装专 业:学院: 广东工业大学考试试卷( A ) 课程名称: 数值计算试卷满分100 分考试时间: 2015 年 12 月 26 日(第 17 周星期六) 题号一二三四五六七八九十总分 评卷得分 评卷签名 复核得分 复核签名 “数值计算”考试要求 “数值计算”考试以开卷形式进行。在“数值计算”课程考试日(2015 年12 月 19 日,第 12 周星期五)考试时间,在考试教室领取试题,在 2015 年12 月 26 日(第 17 周星期六)进行答辩。不参加答辩者将取消考试成绩。 “数值计算”考试结果要求独立在计算机上完成,可使用Matlab或 C 程序编程实现。考试结果将以报告书形式提交,内容包括对问题描述、计算程序以及算例、计算结果、分析组成。计算程序要求具有通用性,能够处理异常情况,可以输入问题、算法参数、算例及初始值,在计算过程中显示当前计算状态、计算完成后显示计算结果。上述内容将作为试卷成绩的主要评定依据。特别提醒,不得使用教师在讲课和实验时的范例作为考试结果。报告书具体格式参考毕业设计手册。 以考生学号命名的文件夹存放程序及报告书电子版,以班级为单位刻录在一张光盘中,与打印版报告书一起由班长和学习委员一起上交任课教师。 数值计算课程总成绩将由试卷成绩(70%)、平时成绩(30%)组成。

2011高考英语模拟题第7套试题

抚顺协作校英语高考专题 2010高考英语模拟题第7套试题 英语试卷 第I卷(三部分,共115分) 第一部分:听力理解(共两节,满分30分) 第一节(共5小题,每小题1.5分,满分7.5分) 听下面5段对话。每段对话后有一个小题,从题中所给的A、B、C三个选项中选出最佳选项,并标在试卷的相应位置。听完每段对话后,你都有10秒钟的时间来回答有关小题和阅读下一小题。每段对话仅读一遍。 1.How much will the woman pay if she buys two pounds of tomatoes ? A.$0.80 B.$2.20 C.$0.30 2.How long does it take the woman to drive home when it isn’t rush hour ? A.Twenty minutes . B.Twenty-five minutes . C.Fifty minutes . 3.What does the man mean about Betty ? A.She doesn’t like going shopping . B.She went shopping yesterday . C.She prefers shopping to studying . 4.What do you know from the conversation ? A.The woman had to work overtime. B.Henry failed to meet the woman . C.The woman had a traffic accident . 5.What does the woman mean ? A.She doesn’t beli eve the man . B.They are not going to land . C.She isn’t afraid . 第二节(共15小题;每小题1.5分,满分22.5分) 听下面5段对话或独白。每段对话或独白后有几个小题,从题中所给的A,B,C三个选项中选出最佳选项,并标在试卷的相应位置。听每段对话或独白前,你将有时间阅读各个小题,每小题5秒钟;听完后,每小题将给出5秒钟的作答时间。每段对话或独白读两遍。听第6段材料,回答第6至8三个小题。 6.What can the woman be ? A.She is a waitress at a hotel . B.She is the secretary of Mr . Li . C.She works at the Blackwood Hotel.

现代数值分析

研硕16《化工数值方法及Matlab应用》试题 班级姓名成绩 1.(15分)数值计算方法的主要研究对象有哪些?其常用基本算法主要包括哪三个方面?举例说明Matlab在解决化工数值计算问题方面有什么样实用价值?答:(1)数值计算方法的主要研究对象为非线性方程求根,插值法、曲线拟合、数值积分、常微分方程(组)、初值问题求解、线性和非线性方程组求解。(2)基本算法包括①离散化方法:用差商代替导数、差分代替微分等,将连续的数学问题转化为离散问题。②逼近方法:用简单函数的值近似代替求解困难或形式未知的复杂函数的值。③迭代法:用一个固定公式反复计算,对较为粗糙的根的近似值进行加工直到满足精度要求的方法。 (3)Matlab在解决化工数值计算问题的实用价值有:数值计算和符号计算功能;图形功能;MATLAB语言;功能性和学科性工具箱。 2.(10分)数值计算中的“曲线拟合”,一般有哪些方法?请至少指出四种,并简述各自的基本特点。 答:(1)拉格朗日插值:,优点在于不要求数据点事等间隔的,缺点是数据点不易过多,当数据比较多时,差值函数有偏离原函数的风险; (2)牛顿插值法:它不仅克服了“增加一个节点时整个计算工作必须重新开始”的缺点,而且可以节省乘、除法运算次数。同时,在牛顿插值多项式中用到的差分与差商等概念,又与数值计算的其他方面有着密切的关系。

(3)牛顿迭代法:牛顿迭代法是求方程根的重要方法之一,其最大优点是在方程f(x) = 0的单根附近具有平方收敛,而且该法还可以用来求方程的重根、复根,此时线性收敛,但是可通过一些方法变成超线性收敛。 (4)区间二分法:优点:算法简单,容易理解,且总是收敛的。缺点:收敛速度太慢,浪费时间,二分法不能求复根跟偶数重根。 (5)最小二乘法:通过最小化误差的平方和寻找数据的最佳函数匹配。利用最小二乘法可以简便地求得未知的数据,并使得这些求得的数据与实际数据之间误差的平方和为最小。 3. (15分)在298K 下,化学反应 2OF 2=O 2+2F 2 的平衡常数为0.410 atm ,如在298K 下将OF 2 通入容器,当t=0 时为1 atm ,问最后总压是多少?取计算精度为10-3 。 解:首先写出求解问题的数学方程式。 假设气体是理想气体,由反应的化学计量式可知, 22222F O OF += 设氧的分压为p ,平衡时有p 21- p p 2。 平衡时,有()410.02142 3=-p p 整理得 0410.064.1640.1423=-+-p p p 函数关系为 ()0410.064.1640.1423=-+-=p p p p f

太原理工大学数值计算方法实验报告

本科实验报告 课程名称:计算机数值方法 实验项目:方程求根、线性方程组的直接解 法、线性方程组的迭代解法、代数插值和最 小二乘拟合多项式 实验地点:行勉楼 专业班级: ******** 学号: ********* 学生姓名: ******** 指导教师:李誌,崔冬华 2016年 4 月 8 日

y = x*x*x + 4 * x*x - 10; return y; } float Calculate(float a,float b) { c = (a + b) / 2; n++; if (GetY(c) == 0 || ((b - a) / 2) < 0.000005) { cout << c <<"为方程的解"<< endl; return 0; } if (GetY(a)*GetY(c) < 0) { return Calculate(a,c); } if (GetY(c)*GetY(b)< 0) { return Calculate(c,b); } } }; int main() { cout << "方程组为:f(x)=x^3+4x^2-10=0" << endl; float a, b; Text text; text.Getab(); a = text.a; b = text.b; text.Calculate(a, b); return 0; } 2.割线法: // 方程求根(割线法).cpp : 定义控制台应用程序的入口点。// #include "stdafx.h" #include"iostream"

心得体会 使用不同的方法,可以不同程度的求得方程的解,通过二分法计算的程序实现更加了解二分法的特点,二分法过程简单,程序容易实现,但该方法收敛比较慢一般用于求根的初始近似值,不同的方法速度不同。面对一个复杂的问题,要学会简化处理步骤,分步骤一点一点的循序处理,只有这样,才能高效的解决一个复杂问题。

2020届全国各地高考英语模拟试卷:书面表达汇总(9套,含范文)

全国各地高考英语书面表达汇总9套 新课标I卷 假定你是李华,暑假想去一家外贸公司兼职,已写好申请书和个人简历(resume)。给外教MS Jenkins 写信,请她帮你修改所附材料的文字和格式(format)。 注意: 1. 词数100左右; 2. 可以适当增加细节,已使行文连贯。 示例1 Dear Ms Jenkins, I'm Li Hua from your English writing class last term. I'm writing to ask for your help. I'm applying for a part-time job at a foreign company in my city during the summer vacation, and I have just completed my application letter and resume. However, I am not quite sure of the lan guage and the format I've used. I know you have a very busy schedule, but I'd be very grateful if you could take some time to go through them and make necessary changes. Please find my application letter and resume in the attachment. Thank you for your kindness! Yours, Li Hua 示例2 Dear Ms Jenkins, I am Li Hua, I am writing to tell you something about my plan for the coming summer vacation and I also want you to do me a favor. In order to get some practical experience, I am planning to take a part-time job in a foreign capital company. I have already finished my job application and personal resume. But this is the first time that I have written an application and the personals resume, so I don’t even know if there are something to pay attention to. So, I’m writing you the letter , hoping you can give me some help. I will be very grateful if you can help me. Looking forward to your reply. And I’d be really thankful. Yours, Li Hua 新课标Ⅱ卷 假定你是李华,你校摄影俱乐部(photography club)将举办国际中学摄影展。请给你的英国朋友Peter写封信。请他提供作品。信的内容包括:

数值分析实验报告总结

数值分析实验报告总结 随着电子计算机的普及与发展,科学计算已成为现代科 学的重要组成部分,因而数值计算方法的内容也愈来愈广泛和丰富。通过本学期的学习,主要掌握了一些数值方法的基本原理、具体算法,并通过编程在计算机上来实现这些算法。 算法算法是指由基本算术运算及运算顺序的规定构成的完 整的解题步骤。算法可以使用框图、算法语言、数学语言、自然语言来进行描述。具有的特征:正确性、有穷性、适用范围广、运算工作量少、使用资源少、逻辑结构简单、便于实现、计算结果可靠。 误差 计算机的计算结果通常是近似的,因此算法必有误差, 并且应能估计误差。误差是指近似值与真正值之差。绝对误差是指近似值与真正值之差或差的绝对值;相对误差:是指近似值与真正值之比或比的绝对值。误差来源见表 第三章泛函分析泛函分析概要 泛函分析是研究“函数的函数”、函数空间和它们之间 变换的一门较新的数学分支,隶属分析数学。它以各种学科

如果 a 是相容范数,且任何满足 为具体背景,在集合的基础上,把客观世界中的研究对象抽 范数 范数,是具有“长度”概念的函数。在线性代数、泛函 分析及相关的数学领域,泛函是一个函数,其为矢量空间内 的所有矢量赋予非零的正长度或大小。这里以 Cn 空间为例, Rn 空间类似。最常用的范数就是 P-范数。那么 当P 取1, 2 ,s 的时候分别是以下几种最简单的情形: 其中2-范数就是通常意义下的距离。 对于这些范数有以下不等式: 1 < n1/2 另外,若p 和q 是赫德尔共轭指标,即 1/p+1/q=1 么有赫德尔不等式: II = ||xH*y| 当p=q=2时就是柯西-许瓦兹不等式 般来讲矩阵范数除了正定性,齐次性和三角不等式之 矩阵范数通常也称为相容范数。 象为元素和空间。女口:距离空间,赋范线性空间, 内积空间。 1-范数: 1= x1 + x2 +?+ xn 2-范数: x 2=1/2 8 -范数: 8 =max oo ,那 外,还规定其必须满足相容性: 所以

大学英语模拟考试试题

郑州大学远程教育学院入学测试机考 专升本大学英语模拟题 1、One of my teeth is so ______ that it is going to be missing soon.() A.lose B.loose C.loss D.lost 标准答案:B 2、—How about taking a walk —Oh, I think it’s ______ cold for a walk. () A.very much B.too much C.much too D.so much 标准答案:C 3、You are supposed to write your composition every ______ line.() A.one B.a C.other D.another 标准答案:C

4、After living for years in a big city, they found it difficult to settle ______ in a town.() A.for B.at C.up D.down 标准答案:D 5、The reason for my absence was ______ I had fallen ill.() A.why B.because C.for D.that 标准答案:D 6、Only with your help, such difficulties in the field. A. can we overcome B. we can overcome C. should we overcome D. We should overcome 标准答案:A 7、The ski resorts are usually crowded. There are many people skiing. A. enjoy B. that want C. who enjoy D. want

MATLAB与数值分析课程总结

MATLAB与数值分析课程总结 姓名:董建伟 学号:2015020904027 一:MATLAB部分 1.处理矩阵-容易 矩阵的创建 (1)直接创建注意 a中括号里可以用空格或者逗号将矩阵元素分开 b矩阵元素可以是任何MATLAB表达式,如实数复数等 c可以调用赋值过的任何变量,变量名不要重复,否则会被覆盖 (2)用MATLAB函数创建矩阵如:a空阵[] b rand/randn——随机矩阵 c eye——单位矩阵 d zeros ——0矩阵 e ones——1矩阵 f magic——产生n阶幻方矩阵等 向量的生成 (1)用冒号生成向量 (2)使用linspace和logspace分别生成线性等分向量和对 数等分向量 矩阵的标识和引用 (1)向量标识 (2)“0 1”逻辑向量或矩阵标识 (3)全下标,单下标,逻辑矩阵方式引用 字符串数组 (1)字符串按行向量进行储存 (2)所有字符串用单引号括起来 (3)直接进行创建 矩阵运算 (1)注意与数组点乘,除与直接乘除的区别,数组为乘方对应元素的幂

(2)左右除时斜杠底部靠近谁谁是分母 (3)其他运算如,inv矩阵求逆,det行列式的值, eig特征值,diag 对角矩阵 2.绘图-轻松 plot-绘制二维曲线 (1)plot(x)绘制以x为纵坐标的二维曲线 plot(x,y) 绘制以x为横坐标,y为纵坐标的二维曲线 x,y为向量或矩阵 (2)plot(x1,y1,x2,y2,。。。。。。)绘制多条曲线,不同字母代替不同颜色:b蓝色,y黄色,r红色,g绿色 (3)hold on后面的pl ot图像叠加在一起 hold off解除hold on命令,plot将先冲去窗口已有图形(4)在hold后面加上figure,可以绘制多幅图形 (5)subplot在同一窗口画多个子图 三维图形的绘制 (1)plot3(x,y,z,’s’) s是指定线型,色彩,数据点形的字 符串 (2)[X,Y]=meshgrid(x,y)生成平面网格点 (3)mesh(x,y,z,c)生成三维网格点,c为颜色矩阵 (4)三维表面处理mesh命令对网格着色,surf对网格片着色 (5)contour绘制二维等高线 (6)axis([x1,xu,y1,yu])定义x,y的显示范围 3.编程-简洁 (1)变量命名时可以由字母,数字,下划线,但是不得包含空格和标点 (2)最常用的数据类型只有双精度型和字符型,其他数据类型只在特殊条件下使用 (3)为得到高效代码,尽量提高代码的向量化程度,避免使用循环结构

2017高职考英语试卷(模拟试题)

2017 年云南省高等职业技术教育招生考试试题(模拟试题) A. are B. be C. were D. is 英语15. I don’t know I should say. 本试题满分100 分,考试时间120 分钟。考生必须在答题卡上答题,在试题纸、草稿纸 A. what B. where C. be D. are 上答题无效。16. he had cleaned the classroom, he left school. I、语音(本大题共 5 小题,每小题 1 分,共 5 分) A. Before B. When C. After D. If 从A、B、C、D 四个选项中找出其画线部分与所给单词的画线部分读音相同的 17. She asked me whether I up with the others already. 选项,并用2B 铅笔把答题卡上对应题目的答案标号涂黑。 A. had held B. was holding C. had caught D. was catching 1. cat 8. This book is interesting than that one. A. take B. made C. have D. same A. much more B. more much C. many more D. more many 2. school 19. You don't have to hurry and there is time. A. foot B. food C. stood D. wood A. few B. a few C. little D. a little 3. hear 20. It's a n u rse ’s duty after patients. A. dear B. hear C. pear D. earth A. to look B. to looking C. looks D. look 4. short 21. Let me congratulate you on the entrance exam. A. word B. work C. worst D. worn A. you pass B. to pass C. passing D. pass 5. honor 22.--W hy isn ’t he in? A. house B. hospital C. hour D. hurry --He to the hospital. Ⅱ、词汇与语法(本大题共25 小题,每小题 1 分,共25 分) A. must go B. can go C. must be D. many have gone 从A、B、C、D 四个选项中选出一个最佳答案,并用2B 铅笔把答题卡上对应题目的答23. Any was the girl who they about all the time. 案标号涂黑。 A. were talked B. had been talking C. had been talked D. had talked 6. There are on the table. 24. There are days In a year. A. two glasses of water B. two glass of milks A. three- hundred-sixth-five B. three hundreds and sixty-five C. a glass of milk D. two glasses of milks C. three hundreds and sixty- fives D. three hundred and sixty-five 7. I am used to TV in the evening. 25. How many English words you the end of last term? A. watch B. watched C. watches D. watching A. do, learn B. are, learned C. had, learned D. are, learning 8. Go and before you do your homework. 26. He told me he forgot all of the . A. has a sleep B. have a sleep C. had a sleep D. slept A. telephone numbers B. telephones number 9. She has never told a lie, . C. telephones number D. telephone number A. has she B. doesn't she C. hasn ’t she D. isn't she 27. September 10th is Day. 10. There is in the word “h o n s e t”. A. Teacher's B. Teachers’ C. Teacher D. Teachers A. a B. the C. an D. X 28. China is famous her Great Wall. 11. Mr. Black arrived in Shanghai 12:45 August20,2007. A. as B. to C. for D. of A. at, on B. to C. at, In D. in, on 29. My mother is ill this week, so . 12. When I passed Mary's house, she the flowers. A. is my father B. does my father C. has my father D. will my father. A. is watering B. was watering C. waters D. watered 30. The glass by my little dog just now. 13.He here if it tomorrow A. broken B. is broken C. is broken D. was broken Ⅲ.完形填空(本大题共10 小题,每小题 1 分,共10 分)A. will stay, will rain B. will stay, rain 阅读下列短文,掌握其大意,然后从每题所给的四个选项中选出最佳选项,并用2B 铅 C. stays, rains D. stayed, rained 14. Neither the girls nor the teacher going to visit Nanjing on Sunday. 笔在答题卡上将该项涂黑。

现代数值分析复习题

复习题(一) 一、填空题: 1、求方程0.5x2 101x 1 0的根,要求结果至少具有6位有效数字。已知 V10203 101.0099,贝卩两个根为x1 _____________________________ , X2 ________________________________ .(要有计算过程和结果) 4 1 0 A A 1 4 1 2、0 1 4,则A的LU分解为。 1 2 A 3、 3 5,贝卩(A) ____________ ,A __________ . 4、已知f(1)「Q f(2)「2,f(3) =3,则用抛物线(辛卜生)公式计算求 3 得1 f(x)dx -------------------- ,用三点式求得f (1) ________________ . 5、f(1) 1,f(2) 2,f(3) 1,则过这三点的二次插值多项式中x2的系数 为_____ ,拉格朗日插值多项式为 _________________________ . 二、单项选择题: 1、Jacobi迭代法解方程组Ax b的必要条件是( ). A. A的各阶顺序主子式不为零 B. (A) 1 C a ii 0,i 1,2, ,n D|| A 1 2、设f(x) 3x99 5x 7,均差f[1,2,22, ,299]=(). D. 3

4、三点的高斯求积公式的代数精度为 ( ). A.3 B. -3 C. 5 D.0 2 2 3 A 0 5 1 3、设 0 0 7 ,则 (A )为( ). A. 2 B. 5 C. 7

分别用拉格朗日插值法和牛顿插值法求 f (x )的三次插值多项式P 3(x ),并 求f (2)的近似值(保留四位小数). 4、 取步长h 0.2,用预估-校正法解常微分方程初值问题 y 2x 3y y (0) 1 (0 x 1) 5、 已知 A. 2 B.5 C. 3 D. 4 5、幕法的收敛速度与特征值的分布 A.有关 B.不一定 C. 无关 三、计算题: 1、用高斯-塞德尔方法解方程组 4X ! 2X 2 X 3 11 X 1 4X 2 2X 3 18 2X ! X 2 5X 3 22 (°) /c c c\T ,取 x (°,°,°),迭 四次(要求按五位有效数字计算 ). 1 2、求A 、B 使求积公式 1 f (X )dX A[f( 1) f (1)] 1 B [f (2)f (2)] 的代数精 度尽量高,并求其代数精度;利用此公式求 I 21dx 1 x (保留四位小 数)。 3、已知

数值分析心得体会

数值分析心得体会 篇一:学习数值分析的经验 数值分析实验的经验、感受、收获、建议班级:计算131 学号:XX014302 姓名:曾欢欢 数值分析实验主要就是学习MATLAB的使用以及对数值分析类容的应用,可以使学生更加理解和记忆数值分析学得类容,也巩固了MATLAB的学习,有利于以后这个软件我们的使用。在做实验中,我们需要具备较好的编程能力、明白MATLAB软件的使用以及掌握数值分析的思想,才能让我们独立自主的完成该作业,如果是上述能力有限的同学,需要借助MATLAB的书以及网络来完成实验。数值分析实验对于我来说还是有一定难度,所以我课下先复习了MATLAB的使用方法以及编写程序的基本类容,借助互联网和同学老师资源完成了数值分析得实验的内容。在实验书写中,我复习了各种知识,所以我认为这门课程是有必要且是有用处的,特别是需要处理大量实验数据的人员,很有必要深入了解学习它,这样在以后的工作学习里面就减少了很多计算问题也提高了实验结果的精确度。 学习数值分析的经验、感受、收获、建议数值分析的内容包括插值与逼近,数值微分与数值积分,非线性方程与线性方程组的数值解法,矩阵的特征值与特征向量计算,常微分方程数值解等。

首先我们必须明白数值分析的用途。通常所学的其他数学类学科都是由公式定理开始,从研究他们的定义,性质再到证明与应用。但实际上,尤其是工程,物理,化学等其它具体的学科。往往我们拿到 手的只是通过实验得到的数据。如果是验证性试验,需要代回到公式 进行分析,验证。但往往更多面对的是研究性或试探性试验,无具体 公式定理可代。那就必须通过插值,拟合等计算方法进行数据处理以得到一个相对可用的一般公式。还有许多计算公式理论上非常复杂,在工程中不实用,所以必须根据实际情况把它转化成多项式近似表 示。学习数值分析,不应盲目记公式,因为公事通常很长且很乏味。其次,应从公式所面临的问题以及用途出发。比如插值方法,就 是就是把实验所得的数据看成是公式的解,由这些解反推出一个近似公式,可以具有局部一般性。再比如说拟合,在插值的基础上考虑实 验误差,通过拟合能将误差尽可能缩小,之后目的也是得到一个具有 一定条件下的一般性的公式。。建议学习本门课程要结合知识与实际,比如在物理实验里面很多

数值计算方法总结计划复习总结提纲.docx

数值计算方法复习提纲 第一章数值计算中的误差分析 1 2.了解误差 ( 绝对误差、相对误差 ) 3.掌握算法及其稳定性,设计算法遵循的原则。 1、误差的来源 模型误差 观测误差 截断误差 舍入误差 2误差与有效数字 绝对误差E(x)=x-x * 绝对误差限x*x x* 相对误差E r (x) ( x x* ) / x ( x x* ) / x* 有效数字 x*0.a1 a2 ....a n10 m 若x x*110m n ,称x*有n位有效数字。 2 有效数字与误差关系 ( 1)m 一定时,有效数字n 越多,绝对误差限越小; ( 2)x*有 n 位有效数字,则相对误差限为E r (x)1 10 (n 1)。 2a1 选择算法应遵循的原则 1、选用数值稳定的算法,控制误差传播; 例 I n 11n x dx e x e I 0 1 1 I n1nI n1 e △ x n n! △x0 2、简化计算步骤,减少运算次数; 3、避免两个相近数相减,和接近零的数作分母;避免

第二章线性方程组的数值解法 1.了解 Gauss 消元法、主元消元法基本思想及算法; 2.掌握矩阵的三角分解,并利用三角分解求解方程组; (Doolittle 分解; Crout分解; Cholesky分解;追赶法) 3.掌握迭代法的基本思想,Jacobi 迭代法与 Gauss-Seidel 4.掌握向量与矩阵的范数及其性质, 迭代法的收敛性及其判定。 本章主要解决线性方程组求解问题,假设n 行 n 列线性方程组有唯一解,如何得到其解? a 11x 1 a 12 x 2... a 1n x n b1 a 21x 1 a 22 x 2... a 2n x n b2 ... a n1x 1 a n 2 x 2... a nn x n b n 两类方法,第一是直接解法,得到其精确解; 第二是迭代解法,得到其近似解。 一、Gauss消去法 1、顺序G auss 消去法 记方程组为: a11(1) x1a12(1) x2... a1(1n) x n b1(1) a21(1) x1a22(1) x2... a2(1n) x n b2(1) ... a n(11) x1a n(12) x2... a nn(1) x n b n(1) 消元过程: 经n-1步消元,化为上三角方程组 a11(1) x1b1(1) a 21(2) x1a22(2 ) x2b2( 2 ) ... a n(1n) x1a n(n2) x2...a nn(n ) x n b n( n ) 第k步 若a kk(k)0 ( k 1)( k) a ik(k )(k )( k 1)( k )a ik(k )( k) a ij a ij a kk(k ) a kj b i b i a kk(k )b k k 1,...n 1 i, j k 1,....,n 回代过程:

英语听力模拟试题共2套

英语听力模拟试题(五) 第一节 听下面5段对话。每段对话后有一个小题,从题中所给的A、B、C三个选项中选出最佳选项,并标在试卷的相应位置。听完每段对话后,你都有10秒钟的时间来回答有关小题和阅读下一小题。每段对话仅读一遍。 1. What does the man advise the woman to do? A. Throw the bread away. B. Eat the bread. C. Take the bread back to the supermarket. 2. What is Peter going to do this afternoon? A. Play basketball. B. Ride a bicycle. C. Go swimming. 3. How much is the man going to pay for the oranges? A. 40 pence. B. 50 pence. C. 60 pence. 4. What does the man think of his classmate? A. He’s stupid. B. He’s strange. C. He’s shameless. 5. What do they decide to do? A. The man will bring some food back for dinner. B. The woman will fill the refrigerator before supper. C. They will eat out. 第二节 听下面5段对话或独白。每段对话或独白后有几个小题,从题中所给的A、B、C三个选项中选出最佳选项,并标在试卷的相应位置。听每段对话或独白前,你将有时间阅读各个小题,每小题5秒钟;听完后,各小题将给出5秒钟的作答时间。每段对话或独白读两遍。 听第6段材料,回答第6至7题。 6. What will the weather be like in the afternoon? A. A bit colder. B. Very warm. C. Less cold. 7. What is the woman most probably going to do this afternoon? A. Go skating. B. Go out by car. C. Make a snowman. 听第7段材料,回答第8至9题。 8. When will the man go to see the doctor? A. On Tuesday. B. On Wednesday. C. On Thursday. 9. What’s wrong with the man? A. He was hit by a ladder. B. He broke his leg. C. He hurt his foot. 听第8段材料,回答第10至12题。 10. In which season might the conversation take place? A. Autumn. B. Winter. C. Spring. 11. Which skirt is the popular one of the year? A. The flowery skirt. B. The green skirt. C. The tan skirt. 12. What is the original price of the clothes? A. 240 dollars. B. 300 dollars. C. 1,200 dollars.

相关主题
文本预览
相关文档 最新文档