Package com.nr.ode

Examples of com.nr.ode.Odeint.integrate()

41
42
43
44
45
46
47
48
49
50
51
    y[0]=1.0; y[1]=1.0; y[2]=0.0;
    Output out=new Output(20);
    rhs_StepperRoss d = new rhs_StepperRoss();
    StepperRoss s = new StepperRoss();
    Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
    ode.integrate();

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
//      System.out.printf(setprecision(12) << yout[i]);
    }
View Full Code Here
40
41
42
43
44
45
46
47
48
49
50
    y[0]=1.0; y[1]=1.0; y[2]=0.0;
    Output out=new Output(20);
    rhs_StepperSie d = new rhs_StepperSie();
    StepperSie s= new StepperSie();
    Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
    ode.integrate();

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
//      System.out.printf(setprecision(12) << yout[i]);
    }
View Full Code Here
50
51
52
53
54
55
56
57
58
59
60
    Output out=new Output(20);
    x1=0.0;x2=1.0;
    rhs_StepperStoerm d = new rhs_StepperStoerm();
    StepperStoerm s1= new StepperStoerm();
    Odeint ode1 = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s1);
    ode1.integrate();
    yexp[0]=bess.airy_ai(1.0);
    yexp[1]=bess.airy_bi(1.0);

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
View Full Code Here
78
79
80
81
82
83
84
85
86
87
88
    dydx[1]=0.0;

    x1=0.0;x2=-1.0;
    StepperStoerm s2 = new StepperStoerm();
    Odeint ode2 = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s2);
    ode2.integrate();
    yexp[0]=bess.airy_ai(-1.0);
    yexp[1]=bess.airy_bi(-1.0);

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
View Full Code Here
45
46
47
48
49
50
51
52
53
54
55
    }
    Output out = new Output(20);
    rhs_StepperDopr5 d = new rhs_StepperDopr5();
    StepperDopr5 s = new StepperDopr5();
    Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
    ode.integrate();

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
      System.out.printf("%f  %f\n", yout[i],yexp[i]);
    }
View Full Code Here
68
69
70
71
72
73
74
75
76
77
78
    yy[3]=y[1].im();
    Hypderiv d = new Hypderiv(a,b,c,z0,dz);
    Output out = new Output();
    StepperBS s = new StepperBS();
    Odeint ode =new Odeint(yy,0.0,1.0,atol,rtol,0.1,0.0,out,d,s);
    ode.integrate();
    y[0]=new Complex(yy[0],yy[1]);
    return y[0];
  }
 
  public static Complex hypgeo(final Complex a, final Complex b,final Complex c,
View Full Code Here
55
56
57
58
59
60
61
62
63
64
65
    h1=(x2-x1)/100.0;
    y=buildVector(load1(x1,v));
    Output out = new Output();
    StepperDopr853 s1 = new StepperDopr853();
    Odeint integ1 = new Odeint(y,x1,xf,atol,rtol,h1,hmin,out,d,s1);
    integ1.integrate();
    f1=buildVector(score(xf,y));
    y=buildVector(load2(x2,v2));
   
    StepperDopr853 s2 = new StepperDopr853();
    Odeint integ2 = new Odeint(y,x2,xf,atol,rtol,h1,hmin,out,d,s2);
View Full Code Here
61
62
63
64
65
66
67
68
69
70
71
    f1=buildVector(score(xf,y));
    y=buildVector(load2(x2,v2));
   
    StepperDopr853 s2 = new StepperDopr853();
    Odeint integ2 = new Odeint(y,x2,xf,atol,rtol,h1,hmin,out,d,s2);
    integ2.integrate();
    f2=buildVector(score(xf,y));
    for (int i=0;i<nvar;i++) f1[i] -= f2[i];
    return f1;
  }
 
View Full Code Here
40
41
42
43
44
45
46
47
48
49
50
    h1=(x2-x1)/100.0;
    y=load(x1,v);
    Output out = new Output();
    StepperDopr853 s = new StepperDopr853();
    Odeint integ = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d, s);
    integ.integrate();
    return score(x2,y);
  }
 
  public abstract double[] load(final double x, final double[] v);
 
View Full Code Here
45
46
47
48
49
50
51
52
53
54
55
    }
    Output out=new Output(20);
    rhs_StepperDopr853 d = new rhs_StepperDopr853();
    StepperDopr853 s = new StepperDopr853();
    Odeint ode = new Odeint(y,x1,x2,atol,rtol,h1,hmin,out,d,s);
    ode.integrate();

    for (i=0;i<nvar;i++) {
      yout[i]=out.ysave[i][out.count-1];
      System.out.printf("%f  %f\n", yout[i],yexp[i]);
    }
View Full Code Here
TOP
Copyright © 2018 www.massapi.com. All rights reserved.
All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.