program para real x(100),y(100) real xp(5),yp(5) i_dummy=pgopen('/xwin') xmin=0 xmax=5 ymin=0 ymax=20 call pgenv(xmin,xmax,ymin,ymax,0,0) dx=(xmax-xmin)/(100-1) do i=1,100 x(i)=xmin+(i-1)*dx y(i)=x(i)*x(i) c y(i)=sqrt(x(i)) end do do i=1,5 xp(i)=i yp(i)=i*i end do call pgline(100,x,y) call pgpt(5,xp,yp,9) call pgclos end