function initial_data(t,N)
global rou T V x dx A prou_pt pV_pt pT_pt;
 dx=3/(N-1);
 x=0:dx:3;
 A=NF_Shape(N);
 rou=zeros(t,N);
 T=zeros(t,N);
 V=zeros(t,N);
 init_condition(N);
 prou_pt=zeros(t,N);
 pV_pt=zeros(t,N);
 pT_pt=zeros(t,N);


function y=NF_Shape(N)
global x;
  y=1+2.2*(x-1.5).^2;
  
function init_condition(N)
global rou T V x;
  rou(1,:)=1-0.3146*x;
  T(1,:)=1-0.2314*x;
  V(1,:)=(0.1+1.09*x).*T(1,:).^(0.5);