I got something a bit different from the other person who responded to you.

function force=Gmmr2Array(mass1, mass2) {

    var G = 6.67408e-11;

    var m1 = mass1;

    var m2 = mass2;

    var r = 1;

    var force = G \* m1 \* m2 / (r \* r);

    return force;

}

function [force, torque]=pointMatrixGravity(array1,array2) {

    var G = 6.67408e-11;

    var m1 = array1(1);

    var m2 = array2(1);

    var r = norm(array1(2:4)-array2(2:4));

    var force = G \* m1 \* m2 / (r \* r);

    var torque = cross(array1(2:4)-array2(2:4), force);

    return [force, torque];

}