i need to solve I(dw/dt)=MR-MG

i am a new to matlab and i need to solve I(dw/dt)=MR-MG this is what i have written so far but not sure how to solve numerically for w(t)

clear all 
close all 
clc

lambda=[14.06;11.71;10.04;8.79;7.81;7.03;6.39;5.86;5.41;5.02;4.69;4.39;4.13;3.90;3.70;3.51;3.35;3.19;3.06;2.93;2.81];
Cp=[.145;.315;.405;.452;.473;.476;.468;.452;.432;.406;.377;.343;.305;.267;.231;.199;.170;.146;.125;.107;.092];

R=30; %meter radius
rho=1.225; %kg/m^3
Vo=10; %m/s
omega=lambda*Vo/R;
I=7.95*10^5 %kgm^2

Mr=0.5*rho*Vo^3*pi*R^2.*Cp./omega; %aero torque Nm

if omega<2.1;
    Mg=0;
else
    Mg=4.35*10^6*(omega-2.1);
end

Answer 1

[Walter Roberson]

Change

if omega<2.1;
    Mg=0;
else
    Mg=4.35*10^6*(omega-2.1);
end

to

Mg = zeros(size(omega));
omega_boundary = 2.1;
mask = omega >= omega_boundary;
Mg(mask) = 4.35*10^6*(omega(mask) - omega_boundary);