plotting p-v diagram of stirling engine using schmidt model

Hi, guys. Im having a problem to plotting a P-V diagram of schmidt model for my final year project. Im currntly on my way to complete my bachelor degree. Supposely I need to plot somthing like in the picture below. the given equation as follows:

1) ∆Q_c=mRT ln⁡(1-(p_A cos⁡(θ+ϕ)∆α)/(p_m+p_A sin⁡( θ+ϕ) ))

2) ∆Q_e=mRT ln⁡(1-(p_A cos⁡(θ-ϕ)∆α)/(p_m+p_A sin⁡(θ-ϕ) ))^(-1)

any variable involve may be random and plus any helps are highly appriciated.

[KALYAN ACHARJYA]

What you have tried, any issue?

[saraahhhhh]

Im a newbie, so actually, i don't know to change the equations into matlab coding.