clear
hw4a=cht2dwp(linspace(0,1,10),linspace(0,1,10),.000001);
hw4a.Gr=1000;
hw4a.Pr=2;
hw4a.T(1,:)=1;
hw4a.T(end,:)=0;
hw4a.prelax=1.3026;

while hw4a.issteady(.05)==0
    hw4a.dopsi();
    hw4a.douv();
    hw4a.doT();
    hw4a.dow();
    hw4a.stept();
end