I have script phase portrait but i got an errors

can't plot a phase portrait

[Jan]

Please mention, what the errors are by posting a copy of the complete message. It is easier to fix a problem than to guess, what the problem is.

Providing parameters as extra input arguments to ODE45 is outdated for 20 years now. Use an anonymous function instead:

[t1,y1] = ode45(@(t,x) @bidangfase(t, x, Pg, Ps, pe, beta1, beta2,a1, a2, a4, a5, c1,c2,c4,c5,v,u,d10,d11,d12,d20,d21,d22,d50,d51,d52,d4,tau,mu, beta3, Phi, psi, delta, gamma, Pr), ...
    t0:0.1:tf, X0, options);