fork download
copy 
  1. import numpy as np
  2.  
  3.  
  4. def f(x):
  5.     return np.exp(x)
  6.  
  7.  
  8. def trapezoidal_rule(a, b, n):
  9.     h = (b - a) / n
  10.     x = np.linspace(a, b, n+1)
  11.     y = f(x)
  12.     integral = h/2 * (y[0] + 2*np.sum(y[1:n]) + y[n])
  13.     return integral
  14.  
  15.  
  16. a, b = 0, 1
  17. n = 1000  
  18. result = trapezoidal_rule(a, b, n)
  19.  
  20. print(f)
Success #stdin #stdout 0.12s 23672KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
<function f at 0x14c13803ef50>
https://ideone.com/r5fii6
language:
Python (cpython 2.7.16)
created:
13 days ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!