在Fortran用数组打印杨辉三角

在Fortran中，打印杨辉三角通常会涉及到循环结构和二维数组的处理。杨辉三角是一个数学术语，它的一行由连续的自然数构成，每个数字等于其上方两个数字之和。下面是一个简单的示例程序，展示如何使用Fortran打印前几行的杨辉三角：

program pythagorean_triangle
    implicit none
    integer :: i, j, n
    integer, dimension(:,:) :: triangle

    ! 初始化数组和行数
    n = 5 ! 可以调整这个值以打印更多的行
    do i = 0, n
        triangle(i+1, 1:i+1) = [(i+1-j, j=0:i), elemental]
    end do

    ! 打印三角形
    do i = 1, n
        do j = 1, i
            write(*, "(I4)", advance="no") triangle(i,j)
        end do
        write(*, *) ! 换行
    end do

end program pythagorean_triangle

在这个程序中，`triangle`是一个二维数组，我们通过嵌套循环填充数组，并使用`write(*, "(I4)")`输出整数。`elemental`关键字用于向量化操作，使得每次循环都能独立计算一行。