(* File: Phi-matrix-rglr   Author: Hongyou Wu   Date: 2000 Dec 02 *)
(*                                                                *)
(* Input: aa, bb, f, q, w, \[Lambda].                             *)
(*                                                                *)
(* Output: pZZ[i,j,t_,\[Lambda]], ph[i,j,t_,\[Lambda]]            *)


Clear[mZZ,pZZ,vZZ,yZZ,zZZ,sZZ];

mZZ=(aa+bb)/2;

sZZ=NDSolve[{yZZ'[vZZ]==zZZ[vZZ]/f[vZZ],
             zZZ'[vZZ]==(q[vZZ]-\[Lambda]*w[vZZ])yZZ[vZZ],
             yZZ[mZZ]==1,
             zZZ[mZZ]==0},
            {yZZ,zZZ},{vZZ,aa,bb}];

pZZ[1,1,vZZ_,\[Lambda]]=(yZZ[vZZ]/.sZZ)[[1]];

pZZ[2,1,vZZ_,\[Lambda]]=(zZZ[vZZ]/.sZZ)[[1]];

sZZ=NDSolve[{yZZ'[vZZ]==zZZ[vZZ]/f[vZZ],
             zZZ'[vZZ]==(q[vZZ]-\[Lambda]*w[vZZ])yZZ[vZZ],
             yZZ[mZZ]==0,
             zZZ[mZZ]==1},
            {yZZ,zZZ},{vZZ,aa,bb}];

pZZ[1,2,vZZ_,\[Lambda]]=(yZZ[vZZ]/.sZZ)[[1]];

pZZ[2,2,vZZ_,\[Lambda]]=(zZZ[vZZ]/.sZZ)[[1]];

ph[1,1,vZZ_,\[Lambda]]=pZZ[1,1,vZZ,\[Lambda]]pZZ[2,2,aa,\[Lambda]]-
                       pZZ[1,2,vZZ,\[Lambda]]pZZ[2,1,aa,\[Lambda]];

ph[1,2,vZZ_,\[Lambda]]=pZZ[1,2,vZZ,\[Lambda]]pZZ[1,1,aa,\[Lambda]]-
                       pZZ[1,1,vZZ,\[Lambda]]pZZ[1,2,aa,\[Lambda]];

ph[2,1,vZZ_,\[Lambda]]=pZZ[2,1,vZZ,\[Lambda]]pZZ[2,2,aa,\[Lambda]]-
                       pZZ[2,2,vZZ,\[Lambda]]pZZ[2,1,aa,\[Lambda]];

ph[2,2,vZZ_,\[Lambda]]=pZZ[2,2,vZZ,\[Lambda]]pZZ[1,1,aa,\[Lambda]]-
                       pZZ[2,1,vZZ,\[Lambda]]pZZ[1,2,aa,\[Lambda]];

Clear[mZZ,yZZ,zZZ,sZZ];