Robotics Toolbox question

I was recently working towards coding an inverse kinematics method, but due to the complexity of the robots "7dof redundant manipulator", i'm working on, i wanted to see other ways to derive joint angles from a series of x,y,z of each joint "motion capture". i'm a self taught in robotics.

I wanted to know, and correct if i'm wrong please, that there are two options, the anlytical, which tries to use transformation matricies entries to deduce the angles, and the one that employs the jacobian "pseudo-inverse rather". The first approach looks very tedious to me because T08 "from end to base trans" is very large when using the DH parameters. so maybe this is the way to do it, i don't know. or the jacobian "computationally demanding".

I have then two questions, which one is suitable for this number of DoF ? does the robotic toolbox have an inverse kinematics that i can set accordingly from DH parameters and xyz data, or just the one available from puma 6dof.

Note: the data file i have, has a stream of 4 joints coordinates in 3D.

please, suggestions is very appreciated