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 |
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 |