/[lmdze]/trunk/dyn3d/coefpoly.f
ViewVC logotype

Diff of /trunk/dyn3d/coefpoly.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 90 by guez, Wed Mar 5 14:57:53 2014 UTC revision 91 by guez, Wed Mar 26 17:18:58 2014 UTC
# Line 18  SUBROUTINE coefpoly(xf1, xf2, xprim1, xp Line 18  SUBROUTINE coefpoly(xf1, xf2, xprim1, xp
18    
19    ! On en revient a resoudre un systeme de 4 equat.a 4 inconnues a0,a1,a2,a3    ! On en revient a resoudre un systeme de 4 equat.a 4 inconnues a0,a1,a2,a3
20    
21    DOUBLE PRECISION xf1, xf2, xprim1, xprim2, xtild1, xtild2, xi    DOUBLE PRECISION xf1, xf2, xprim1, xprim2, xtild1, xtild2
   DOUBLE PRECISION xfout, xprim  
22    DOUBLE PRECISION a1, a2, a3, a0, xtil1car, xtil2car, derr, x1x2car    DOUBLE PRECISION a1, a2, a3, a0, xtil1car, xtil2car, derr, x1x2car
23    
24    xtil1car = xtild1*xtild1    xtil1car = xtild1*xtild1
# Line 35  SUBROUTINE coefpoly(xf1, xf2, xprim1, xp Line 34  SUBROUTINE coefpoly(xf1, xf2, xprim1, xp
34    a1 = xprim1 - 3.*a3*xtil1car - 2.*a2*xtild1    a1 = xprim1 - 3.*a3*xtil1car - 2.*a2*xtild1
35    a0 = xf1 - a3*xtild1*xtil1car - a2*xtil1car - a1*xtild1    a0 = xf1 - a3*xtild1*xtil1car - a2*xtil1car - a1*xtild1
36    
   RETURN  
37  END SUBROUTINE coefpoly  END SUBROUTINE coefpoly
Legend:
Removed from v.90  
changed lines
  Added in v.91
  ViewVC Help
Powered by ViewVC 1.1.21