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);
