public FixedPoint(int N,double limiteInferior,double limiteSuperior,
Integral.Funcao bf, Integral.Funcao an, FuncaoDupla fp,
Integral.Funcao D, Integral.Funcao funcaoInicial)
{
autoQuad = new AutoQuad(10, limiteInferior, limiteSuperior, precisaoRequerida, 10);
this.funcaoInicial = funcaoInicial;
this.an = an;
this.bf = bf;
this.fp = fp;
this.D = D;
this.N = N;
x0 = limiteInferior;
xf = limiteSuperior;
}
public void AtribuirCampos()
{
metodoDeCoeficientes = MomentosBase.MetodoCoeficientesDosMomentos.Wheeler;
funcaoDeMomentoInicial = FuncaoDeMomentoInicial;
autoQuad = new AutoQuad(10, 0.0, 1.0, precisaoRequerida, 10);
numeroDePontos = 7;
deltaN = 4;
Lmin = 0.0;
Lmax = 1.0;
limiteInferiorDistribuicao = Lmin;
limiteSuperiorDistribuicao = Lmax;
//Parâmetros do Hidrocraqueamento
//Valores dos parâmetros para T = 324 C
alfa = 1.35;
//kMax = 1.35; //1/hora;
kMax = 1.35;
a0 = 6.41;
a1 = 28.15;
delta = 2.6667e-05;
//Valores dos parâmetros para T = 342 C
//alfa = 0.40;
//kMax = 8.08; //1/hora;
//a0 = 6.00;
//a1 = 3.80;
//delta = 9.05e-07;
//Valores dos parâmetros para T = 354 C
//alfa = 0.314;
//kMax = 32.08; //1/hora;
//a0 = 5.82;
//a1 = 3.60;
//delta = 9.05e-07;
funcaoEVezesPolinomio = FuncaoEVezesPolinomio;
}
public void AtribuirCampos()
{
metodoDeCoeficientes = MomentosBase.MetodoCoeficientesDosMomentos.Wheeler;
funcaoDeMomentoInicial = FuncaoDeMomentoInicial;
autoQuad = new AutoQuad(10, 0.0, 1.0, precisaoRequerida, 10);
autoQuadS0 = new AutoQuad2(10,0.0,1.0,precisaoRequerida,10);
numeroDePontos = 5;
deltaN = 2;
numeroDePontosInicial = 5;
deltaNInicial = 2;
//Parâmetros da distribuição inicial
x0 = 6.0;
aP = 99.7373819939209;
alfaGama = 2.11974361925893;
beta = 8.60867323342265;
Lmin = 6.0;
Lmax = 69.5;
limiteInferiorDistribuicao = Lmin;
limiteSuperiorDistribuicao = Lmax;
//Parâmetros do Hidrocraqueamento
//Valores dos parâmetros para T = 324 C
alfa = 0.5;
// alfa = 0.35;
kMax = 2.64; //1/hora;
a0 = 8.14;
a1 = 4.68;
delta = 7.05e-07;
//Valores dos parâmetros para T = 342 C
//alfa = 0.40;
//kMax = 8.08; //1/hora;
//a0 = 6.00;
//a1 = 3.80;
//delta = 9.05e-07;
//Valores dos parâmetros para T = 354 C
//alfa = 0.314;
//kMax = 32.08; //1/hora;
//a0 = 5.82;
//a1 = 3.60;
//delta = 9.05e-07;
funcaoEVezesPolinomio = FuncaoEVezesPolinomio;
}