I'm trying to set up a Matlab M file to plot the results from a circuit.

The circuit consists of a source (14.4kVRMS) a line (Z= .75+j20 ohms) and a load(which varies)

I would like to be able to plot the VA across the load, versus the power factor of the load. Now the rub is that I must keep the load voltage within +-10% of the source voltage. I attached the M-file that I have pulled together so far.

Thanks for your help.