function mrp=morph(f1,f2,s)
mrp=s.*f1+(1-s).*f2;
