Using Matlab to show magnetic field of our earth

Hi there, i need to write a Matlab program to show the magnetic field of our earth form given mathematical equation provided by International Geomagnetic Reference(IGRF). Given Tu= du/36525 where du is the number days of Universal time elapsed since JD 2451545.0, Zonal coefficient = 0.00108263 and earth rotation rate= 7.2921151 x 10^-5 rad/s. plz do help. tq

Answer 1

[Arkadiy Turevskiy]

Maybe this could be a starting point?

Arkadiy