【Java开源代码栏目提醒】:本文主要为网学会员提供d3r8.java,希望对需要d3r8.java网友有所帮助,学习一下!
import java.text.*;
public class d3r8
{
public static void main (String[] args)
{
//program d3r8
//driver for routine gauleg
int i,npoint = 10;
double x1 = 0.0;
double x2 = 1.0;
double x3 = 10.0;
double x[] = new double[11];
double w[] = new double[11];
d3r8F g = new d3r8F();
DecimalFormat form = new DecimalFormat("0.000000");
g.gauleg(x1, x2, x, w, npoint);
System.out.println("# x(i) w(i)");
for (i = 1; i <= npoint; i++)
{
System.out.print(i + " ");
System.out.print(form.format(x[i]) + " ");
System.out.println(form.format(w[i]));
}
//demonstrate the use of gauleg for an integral
g.gauleg(x1, x3, x, w, npoint);
double xx = 0.0;
for (i = 1; i <= npoint; i++)
{
xx = xx + w[i] * g.func(x[i]);
}
System.out.print("integral from gauleg: ");
System.out.println(form.format(xx));
System.out.print("actual value: ");
System.out.println(form.format((1.0 - (1.0 + x3) * Math.exp(-x3))));
}
}