integer::I,J,N
real*8::Cr,T
real*8,dimension(:),allocatable ::P,X
N=3 !變量X的個數
Cr=5.0d0 !常量Cr，自己設定
T=4.0d0 !常量T，自己設定
allocate(P(N),X(N))
! =======讀入變量X的值
do I=1,N
write(*,*)"請輸入第 ",I," 個變量的值："
read(*,*)X(I)
enddo
! =======讀入變量X的值
do I=1,N
P(I)=(-4.2d0/Cr**2*X(I)+2.9/Cr)*T
write(*,*)“第 ”,I," 個變量X對應結果：",P(I)
enddo
end