Conduct a Breusch Pagan test for a multiple regression?

Hi,

I want to conduct a Breusch Pagan test for a multiple regression in mathematica. Attached I got the formula for the Breusch Pagan test we used for a single regression. How can I modify this formula so that I get a p-value for the BP test for a multiple regression?

Thanks for your help!
Nora

Please post a full example as code (not an image, that is).