program demo_plottuse M_datapac,only:plott,labelimplicit noneinteger::iinteger,parameter::dp=kind(0.0d0)real(kind=dp),allocatable::x(:),y(:)call label('plott')y=[(real(i)/10.0,i=1,30)]x=y**3.78-6*y**2.52+9*y**1.26call plott(x,y,size(x))end program demo_plott