I have a moneynessratio-tenor volatility surface and want to extrapolate the implied volatility for moneynessratios > 150%. The volatility surface was downloaded for different points in time, so I basically have one matrix for every moneyness ratio which holds per column a different (standard) tenor. These matrices were merged together into one big matrix to be able to process the data. This table looks as follows: ![Date 50.7 50.14 .... 60.7 2019-02-01 0.67 0.68 0.97 2019-02-02 0.43 0.26 0.26 2019-02-03 0.69 0.66 0.13] 1 To extrapolate I want to use the quadradic polynom of the form IV = MN² + MN + TN² + TN + alpha with MN = Moneynessratio, TN = Tenor, IV = Implied Volatility, alpha = Intercept. MN is a matrix that contains the moneynessratios which is the same for every day (e.g. 50%) TN is a matrix that contains the tenor which is the same for every day (e.g. 7). If I run now a regression, I regress the 1st column of every explanatory variable on the 1st column of the indepenendt variable and get two betas (for every column of the explanatory variables one beta). After the regression is done, I have a two matrices of betas. How can I now use them to extrapolate for a non-standard moneyness ratio (e.g. 168 %) and a non standard tenor (e.g. 9 days)? Thank you very much for your help!