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