Abbreviation of evaluated conflicts with built-in function eval about typst-physics HOT 2 OPEN

#import "@preview/cetz:0.2.0": canvas, plot, draw

#let E-hermitepsi = ((0, "1/(calc.sqrt(calc.sqrt(calc.pi)))*(1 * f)"), (1, "1/(calc.sqrt(calc.sqrt(4*calc.pi)))*(2 * y * f)"), (2, "1/(calc.sqrt(calc.sqrt(64*calc.pi)))*(4 * y * y - 2) * f"), (3, "1/(calc.sqrt(calc.sqrt(2304*calc.pi)))*(8 * y * y * y - 12 * y)*f"))
#figure(
  grid(
    columns: 2,     // 2 means 2 auto-sized columns
    gutter: 2mm,    // space between columns
    canvas(length: 1.2cm, {
      import draw: *
      plot.plot(
        size: (5, 10),
        x-tick-step: 2,
        x-label: [$y$],
        // x-ticks: ( (0, $0$), (1, $a$)),
        y-label: [$V(x)$, #text([$E$], fill:blue), #text([$psi(y)$], fill:red)],
        y-tick-step: none,
        y-ticks: ((0.5, $v=0$), (1.5, $v=1$), (2.5, $v=2$), (3.5, $v=3$)),
        x-min: -4,
        x-max: 4,
        y-min: 0,
        // y-max: 12,
        {  
          plot.add(style: (stroke: (paint:black)),
            domain: (-3, 3), samples: 100, y => 0.5*y*y
          )
          for (v, hermitepsi) in E-hermitepsi {
            let E = v + 1/2
            let psi(y) = eval(hermitepsi, scope: (y:y, f: calc.exp(-y*y/2)))
              plot.add(style: (stroke: (paint:blue)),
                domain: (-4, 4), samples: 10, y => E
              )
              plot.add(style: (stroke: (paint:red)),
               domain: (-4, 4), samples: 100, y => psi(y) + E
              )
          }
        }
      )
      line(style: (stroke: (thickness: 5pt)), (2.5, 0), (2.5, 10))  
    }),
    canvas(length: 1.2cm, {
      import draw: *
      plot.plot(
        size: (5, 10),
        x-tick-step: 2,
        x-label: [$y$],
        y-label: [$V(x)$, #text([$E$], fill:blue), #text([$psi^2(y)$], fill:red)],
        y-tick-step: none,
        y-ticks: ( (0.5, $v=0$), (1.5, $v=1$),(2.5, $v=2$),(3.5, $v=3$)),
        x-min: -4,
        x-max: 4,
        y-min: 0,
        // y-max: 12,
        {  
          plot.add(style: (stroke: (paint:black)),
            domain: (-3, 3), samples: 100, y => 0.5*y*y
          )
          for (v, hermitepsi) in E-hermitepsi {
            let E = v + 1/2
            let psi(y) = eval(hermitepsi, scope: (y:y, f: calc.exp(-y*y/2), E: E))
            plot.add(style: (stroke: (paint:blue)),
              domain: (-4, 4), samples: 10, y => E
            )
            plot.add(style: (stroke: (paint:red)),
             domain: (-4, 4), samples: 100, y => psi(y) * psi(y) +  E
            )
          }
        }
      )
      line(style: (stroke: (thickness: 5pt)), (2.5, 0), (2.5, 10))
    })
  )
)