Question

用Python求积分

import math

import sys

from texttable import Texttable

import time

start = time.clock()

def area()

double x1=1;//起始区间

double x2=4.5;//结束区间

double dx=0.0001 ;//步长，步长越小越精确

double y1=0,y2=0,y3=0;//结果

for(double x = x1;x<x2;x=x+dx)

y1=y1+dx*x*x;//细小矩形取左边为高

y2=y2+ dx*(x+0.5*dx)*(x+0.5*dx); //取中间

y3=y3+dx*(x+dx)*(x+dx);

cout<<y1<<" "<<y2<<" "<<y3<<endl;

long clock_t clock(void);

printf(clock_t);
这段代码编译不出来，报错：第七行invlid syntax

Answer 1

科学计算应该要安装numpy，你可以去下Anaconda，全套自动打包好了的。

把代码改了一下：

# -*- coding: utf-8 -*-
import math
import sys
import time
import numpy as np

def area():
    x1=1            #起始区间
    x2=4.5            #结束区间
    dx=0.0001        #步长，步长越小越精确
    y1=0
    y2=0
    y3=0            #结果
    for x in np.arange(x1, x2, dx):
        y1=y1+dx*x*x            #细小矩形取左边为高
        y2=y2+ dx*(x+0.5*dx)*(x+0.5*dx)             #取中间
        y3=y3+dx*(x+dx)*(x+dx)
        print (y1, y2, y3)
        
start = time.clock()        
area()
end = time.clock()
print (end-start)

Answer 2

import sys
import time

start = time.clock()
def area():
    x1, x2, dx=1., 4.5, 0.0001
    y1, y2, y3 = 0., 0., 0.
    x = x1
    while x < x2:
        y1 += dx*x*x
        y2 += dx*(x+0.5*dx)*(x+0.5*dx)
        y3 += dx*(x+dx)*(x+dx)
        x += dx
        print "%.4f  %.4f  %.4f" % (y1, y2, y3)

area()
clock_t = time.clock()
print clock_t - start

Answer 3

python不能用双斜线注释的吧，注释要用 #

追问

我去掉注释也没通啊

追答

哦，分号也不能用，python是不需要用分号来结尾的