java 版本二重积分计算

public abstract class Calculate {
	public Calculate(){
		
	}
	public abstract double function(double x);
	
	public abstract double function(double x,double y);
	
	/**
	 * 牛顿-莱布尼兹积分
	 * */
	public double Integral(double x1,double x2,double eps){
		if(x2==x1){
			try {
				throw new Exception("不好意思,貌似不可以这么积分");
			} catch (Exception e) {
				e.printStackTrace();
			}
		}
		double length=Math.abs(x1-x2);
		int n=(int)(1/eps);
		double steps=length/n;
		double sum=0;
		for(int i=0;i<n;i++){
			double height=function(x1+steps*i);
			double square=steps*height;
			sum=sum+square;
		}
		return sum;
	}
	public double Integral(double x1,double x2,double y1,double y2,double eps){
		double lengthY=Math.abs(y2-y1);
		double lengthX=Math.abs(x2-x1);
		int n=(int)(1/eps);
		double stepsY=lengthY/n;
		double stepsX=lengthX/n;
		double sum=0;
		for(int i=0;i<n;i++){
			for(int j=0;j<n;j++){
				double heightY=function(x1+stepsX*i,y1+stepsY*j);
				double squareY=heightY*stepsY*stepsX;
				sum=sum+squareY;
			}
		}
		return sum;
	}
}

public class test{
	public static void main(String[] args){
		class C extends Calculate{

			@Override
			public double function(double x) {
				// TODO Auto-generated method stub
				return x;
			}

			@Override
			public double function(double x, double y) {
				return E.e(x)*E.e(y);
			}
		}
		C c=new C();
		double r=c.Integral(0, 1,0,1,0.00005);
		System.out.println((E.e(1)-1)*(E.e(1)-1));
		System.out.println(r);
	}
}

评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值