0% found this document useful (0 votes)
4 views3 pages

Codingan Kalkulus

The document contains a Python script that calculates Riemann sums for various functions using different sampling methods. It includes six problems, each demonstrating how to compute the Riemann sum for specified intervals and sample points, with outputs showing the values of the function, interval widths, and areas. The results for each problem are printed, including the final Riemann sum for each calculation.

Uploaded by

nazrielalfian024
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
4 views3 pages

Codingan Kalkulus

The document contains a Python script that calculates Riemann sums for various functions using different sampling methods. It includes six problems, each demonstrating how to compute the Riemann sum for specified intervals and sample points, with outputs showing the values of the function, interval widths, and areas. The results for each problem are printed, including the final Riemann sum for each calculation.

Uploaded by

nazrielalfian024
Copyright
© All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd

3/4/26, 9:59 PM Untitled0.

ipynb - Colab

import numpy as np

def f1(x):
return -x**2 + 4*x

# Subinterval dan titik sampel


intervals_1 = [(1, 2.5), (2.5, 3.5), (3.5, 5)]
sample_points_1 = [2, 3, 4.5]

riemann_sum_1 = 0

print("SOAL 1")
print("-------------------------------------------------")
print("xi\tf(xi)\tΔx\tf(xi)Δx")

for i in range(len(intervals_1)):
a, b = intervals_1[i]
dx = b - a
xi = sample_points_1[i]
fx = f1(xi)
area = fx * dx
riemann_sum_1 += area

print(f"{xi:.2f}\t{fx:.4f}\t{dx:.2f}\t{area:.4f}")

print("-------------------------------------------------")
print("Riemann Sum =", round(riemann_sum_1,4))

def f2(x):
return x**2 - 4*x + 3

intervals_2 = [(0,0.5), (0.5,2), (2,4)]


sample_points_2 = [0.7, 1.7, 3.5]

riemann_sum_2 = 0

print("\nSOAL 2")
print("-------------------------------------------------")
print("xi\tf(xi)\tΔx\tf(xi)Δx")

for i in range(len(intervals_2)):
a, b = intervals_2[i]
dx = b - a
xi = sample_points_2[i]
fx = f2(xi)
area = fx * dx
riemann_sum_2 += area

print(f"{xi:.2f}\t{fx:.4f}\t{dx:.2f}\t{area:.4f}")

print("-------------------------------------------------")
print("Riemann Sum =", round(riemann_sum_2,4))

def f3(x):
return x - 1

partition_3 = [3, 3.75, 4.25, 5.5, 6, 7]


sample_points_3 = [3, 4, 4.75, 6, 6.5]

riemann_sum_3 = 0

print("\nSOAL 3")
print("-------------------------------------------------")
print("xi\tf(xi)\tΔx\tf(xi)Δx")

for i in range(len(sample_points_3)):
dx = partition_3[i+1] - partition_3[i]
xi = sample_points_3[i]
fx = f3(xi)
area = fx * dx
riemann_sum_3 += area

print(f"{xi:.2f}\t{fx:.4f}\t{dx:.2f}\t{area:.4f}")

print("-------------------------------------------------")
i t("Ri S " d( i 3 4))
[Link] 1/3
3/4/26, 9:59 PM [Link] - Colab
print("Riemann Sum =", round(riemann_sum_3,4))

def f4(x):
return -x/2 + 3

partition_4 = [-3, -1.3, 0, 0.9, 2]


sample_points_4 = [-2, -0.5, 0, 2]

riemann_sum_4 = 0

print("\nSOAL 4")
print("-------------------------------------------------")
print("xi\tf(xi)\tΔx\tf(xi)Δx")

for i in range(len(sample_points_4)):
dx = partition_4[i+1] - partition_4[i]
xi = sample_points_4[i]
fx = f4(xi)
area = fx * dx
riemann_sum_4 += area

print(f"{xi:.2f}\t{fx:.4f}\t{dx:.2f}\t{area:.4f}")

print("-------------------------------------------------")
print("Riemann Sum =", round(riemann_sum_4,4))

def f5(x):
return x**2/2 + x

a, b = -2, 2
n = 8
dx = (b - a) / n

riemann_sum_5 = 0

print("\nSOAL 5 (Midpoint)")
print("-------------------------------------------------")
print("xi\tf(xi)\tΔx\tf(xi)Δx")

for i in range(n):
x_left = a + i*dx
x_right = x_left + dx
xi = (x_left + x_right)/2
fx = f5(xi)
area = fx * dx
riemann_sum_5 += area

print(f"{xi:.4f}\t{fx:.4f}\t{dx:.4f}\t{area:.4f}")

print("-------------------------------------------------")
print("Riemann Sum =", round(riemann_sum_5,4))

def f6(x):
return 4*x**3 + 1

a, b = 0, 3
n = 6
dx = (b - a) / n

riemann_sum_6 = 0

print("\nSOAL 6 (Right Endpoint)")


print("-------------------------------------------------")
print("xi\tf(xi)\tΔx\tf(xi)Δx")

for i in range(1, n+1):


xi = a + i*dx
fx = f6(xi)
area = fx * dx
riemann_sum_6 += area

print(f"{xi:.4f}\t{fx:.4f}\t{dx:.4f}\t{area:.4f}")

print("-------------------------------------------------")
print("Riemann Sum =", round(riemann_sum_6,4))

[Link] 2/3
3/4/26, 9:59 PM [Link] - Colab

SOAL 1
-------------------------------------------------
xi f(xi) Δx f(xi)Δx
2.00 4.0000 1.50 6.0000
3.00 3.0000 1.00 3.0000
4.50 -2.2500 1.50 -3.3750
-------------------------------------------------
Riemann Sum = 5.625

SOAL 2
-------------------------------------------------
xi f(xi) Δx f(xi)Δx
0.70 0.6900 0.50 0.3450
1.70 -0.9100 1.50 -1.3650
3.50 1.2500 2.00 2.5000
-------------------------------------------------
Riemann Sum = 1.48

SOAL 3
-------------------------------------------------
xi f(xi) Δx f(xi)Δx
3.00 2.0000 0.75 1.5000
4.00 3.0000 0.50 1.5000
4.75 3.7500 1.25 4.6875
6.00 5.0000 0.50 2.5000
6.50 5.5000 1.00 5.5000
-------------------------------------------------
Riemann Sum = 15.6875

SOAL 4
-------------------------------------------------
xi f(xi) Δx f(xi)Δx
-2.00 4.0000 1.70 6.8000
-0.50 3.2500 1.30 4.2250
0.00 3.0000 0.90 2.7000
2.00 2.0000 1.10 2.2000
-------------------------------------------------
Riemann Sum = 15.925

SOAL 5 (Midpoint)
-------------------------------------------------
xi f(xi) Δx f(xi)Δx
-1.7500 -0.2188 0.5000 -0.1094
-1.2500 -0.4688 0.5000 -0.2344
-0.7500 -0.4688 0.5000 -0.2344
-0.2500 -0.2188 0.5000 -0.1094
0.2500 0.2812 0.5000 0.1406
0.7500 1.0312 0.5000 0.5156
1.2500 2.0312 0.5000 1.0156
1.7500 3.2812 0.5000 1.6406
-------------------------------------------------
Riemann Sum = 2.625

SOAL 6 (Right Endpoint)


-------------------------------------------------
xi f(xi) Δx f(xi)Δx
0.5000 1.5000 0.5000 0.7500
1.0000 5.0000 0.5000 2.5000

[Link] 3/3

You might also like