y2[x_]=D[y1[x],x]
     3 + 6 x