labNames <- c('xLab','yLabl')

plot(c(1:10),xlab=expression(paste(labName[1], x^2)),ylab=expression(paste(labName[2], y^2)))

labNames <- c('xLab','yLab')

xlab <- bquote(.(labNames[1]) ~ x^2)

ylab <- bquote(.(labNames[2]) ~ y^2)

plot(c(1:10), xlab = xlab, ylab = ylab)

labNames <- c('xLab','yLab')

plot(c(1:10),

    xlab=substitute(paste(nn, x^2), list(nn=labNames[1])),

    ylab=substitute(paste(nn, y^2), list(nn=labNames[2])))

par(mar=c(5, 4.3, 4, 2) + 0.1)

plot(c(1:10), xlab = expression(xLab ~ x^2 ~ m^-2),

    ylab = expression(yLab ~ y^2 ~ m^-2),

    main="Plot 1")

plot(c(1:10), xlab = expression(xLab * x^2 * m^-2),

    ylab = expression(yLab * y^2 * m^-2),

    main="Plot 2")

plot(c(1:10), xlab = expression(xLab ~ x^2 * m^-2),

    ylab = expression(yLab ~ y^2 * m^-2),

    main="Plot 3")

plot(c(1:10), xlab = expression('Your x label' ~ mu[3] * phi),

    ylab = expression("Temperature (" * degree * C *")"))