The soil-water coupled analysis program of the invention computes a global L matrix and a modified global L matrix regarding a volume change rate of a soil skeleton over time, a global H matrix regarding a water permeability of soil, a global M matrix regarding a mass, and a global K matrix regarding a tangent stiffness of the soil skeleton, based on input settings of soils, such as clay, intermediate soil, and sand, to respective elements of a soil foundation, input settings of a solid soil model, and input settings of analysis conditions (steps S140 to S165). The soil-water coupled analysis program formulates a global tangent stiffness equation (simultaneous linear equations) using all these computed matrixes and determines an unknown ‘jerk field’ and a ‘pore water pressure field’ under given boundary conditions, for example, a given deformation condition and a given stress rate condition (step S170). This enables highly-accurate dynamic and static analyses in soil foundations of various soils from sand to intermediate soils and clay.