How can I go about calling another function and estimating a solution?

I have attached a matlab function, that returns a solution to a given system of equations with naive Gaussian elimination. I would like to write a script to call the function and estimate the solution for the following system of equations (attached image) with single presicion.

How can I go about doing that?

Answer 1

[Stephan]

Your new script:

A = [3, 123456789, 8, 5; 2, 13,...];
b = [246913625, ...];
[x, det] = single(nvgss(A,b))

Best regards

Stephan