simulate The sun shines on the bubbles and turns them into colors

Here’s an example program in Python that simulates the sun shining on bubbles and turning them into colors:

pythonCopy codeimport math
import numpy as np
import matplotlib.pyplot as plt

def bubble_colors(n=1000, radius_mean=1, radius_std=0.1, wavelengths=None, sun_pos=(0, 0)):
    """
    Simulate sun shining on bubbles and turning them into colors.
    :param n: number of bubbles
    :param radius_mean: mean radius of bubbles
    :param radius_std: standard deviation of bubble radius
    :param wavelengths: wavelengths of light to simulate (default: RGB)
    :param sun_pos: position of sun (default: (0, 0))
    :return: (x, y, colors) where x and y are the positions of the bubbles, and colors is a list of RGB tuples
    """
    if wavelengths is None:
        # RGB wavelengths
        wavelengths = [(1, 0, 0), (0, 1, 0), (0, 0, 1)]

    # Generate bubble radii and positions
    radii = np.random.normal(radius_mean, radius_std, size=n)
    thetas = np.random.uniform(0, 2*math.pi, size=n)
    x = radii * np.cos(thetas)
    y = radii * np.sin(thetas)

    # Calculate distance from sun for each bubble
    distances = np.sqrt((x-sun_pos[0])**2 + (y-sun_pos[1])**2)

    # Calculate angle of incidence and angle of reflection for each wavelength
    angles_i = np.arcsin(distances / radii)
    angles_r = angles_i

    # Calculate phase shift for each wavelength
    phase_shifts = [2*math.pi*d/wavelength for d in distances for wavelength in wavelengths]

    # Calculate color of each bubble for each wavelength
    colors = []
    for i in range(n):
        color = np.zeros(3)
        for j, wavelength in enumerate(wavelengths):
            phase = phase_shifts[i*len(wavelengths) + j]
            color += np.array(wavelength) * np.cos(angles_i[i] - angles_r[i] + phase)
        color /= len(wavelengths)
        colors.append(tuple(color))

    return x, y, colors

# Generate bubbles and plot them
x, y, colors = bubble_colors()
plt.scatter(x, y, c=colors)
plt.axis('equal')
plt.show()

This program generates a set of bubbles with random positions and radii, and then calculates the color of each bubble based on the wavelengths of light and the position of the sun. The colors are then plotted using matplotlib.

Note that this is a simplified simulation and does not take into account many factors that can affect the color of bubbles in real life, such as the angle of the sun, the refractive index of water, and the size distribution of bubbles in the water.