call('python C:\OSGeo4W64\bin\gdal_calc.py -A {0} -B {1} -C {2} --outfile={3} --calc="(A*B+(10*C)"'.format(data1,data2,data3,myout))

data = gdal.Open(mydata1)
data1 = data.ReadAsArray()

data2 = gdal.Open(mydata2)
data3 = data2.ReadAsArray()

data4 = gdal.Open(mydata3)
data5 = data4.ReadAsArray()

data_uri1 = array(data1, dtype=float)
data_uri2 = array(data3, dtype=float)
data_uri3 = array(data5, dtype=float)

calc = (data_uri1*data_uri2+(10*data_uri3))
geo = data.GetGeoTransform()
proj = data.GetProjection()
shape = data1.shape

driver = gdal.GetDriverByName("GTiff")
dst_ds = driver.Create(tdi_path, shape[1], shape[0], 1, gdal.GDT_Float32)
dst_ds.SetGeoTransform(geo)
dst_ds.SetProjection(proj)
dst_ds.GetRasterBand(1).WriteArray(calc)
dst_ds = None