Thanks for response ...
But this method is highly error generating... As you can see that initially it produces the surface
Where as after interpolation it produces the Surface
Which shohs a huge error in Z axis...
You can use Interpolation, after converting your triples {x,y,z} into the structure {{x,y},z}:
f = Interpolation[root1 /. {x_Real, y_, z_} :> {{x, y}, z}]; Plot3D[f[x, y], {x, .1, .4}, {y, .1, 9}]