Lab 06/05/25

# Define variables
depth <- 300        # in meters
density <- 1025     # in kg/m³
gravity <- 9.8      # in m/s²

# Calculate raw pressure
pressure <- density * gravity * depth

# Subtract 100,000 to account for air compression inside the sub
adjusted_pressure <- pressure - 100000

# Take the square root of adjusted pressure
sqrt_adjusted <- sqrt(adjusted_pressure)

# Raise that result to the power of 1.1
hull_stress <- sqrt_adjusted^1.1

# Print all results
print(paste("Raw pressure:", pressure, "Pascals"))

[1] "Raw pressure: 3013500 Pascals"

print(paste("Adjusted pressure:", adjusted_pressure, "Pascals"))

[1] "Adjusted pressure: 2913500 Pascals"

print(paste("Square root of adjusted pressure:", sqrt_adjusted))

[1] "Square root of adjusted pressure: 1706.89777081113"

print(paste("Final hull stress estimate:", hull_stress))

[1] "Final hull stress estimate: 3592.76045889325"