As @imi-miri points out, you can import stan but I realized that the API is quite different and I had to go back and forth with the AIs to even get some test code to work.

This runs:

import stan

model_code = """
data {
  int<lower=0> N;
  array[N] real y;  // NEW SYNTAX for arrays
}
parameters {
  real mu;
}
model {
  y ~ normal(mu, 1);
}
"""

# Compile and fit the model
posterior = stan.build(model_code, data={'N': 10, 'y': [1.2, 2.4, 1.5, 0.9, 3.2, 1.8, 2.7, 3.1, 1.6, 2.0]})
fit = posterior.sample(num_chains=4, num_samples=1000)

print(fit)