<?php 

% CDF & pdf of the power reserve costs

%اين برنامه ابتدا تابع چگالي تجمعي هزينه را محاسبه مي کند که به صورت تابعي

%با پارامتر Csr خواهد شد و سپس با مشتق گيري از آن تابع چگالي احتمال را

%محاسبه مي کند و در انتها آن را رسم مي کند.

syms c P Csr

z=7.0092e-057*((c^44)*exp(-(c/0.2)))*exp(-((P^2)/200));

CDF=int(int(z,c,0,(Csr/P)+2.5),P,0,+inf);

Csr=-5:1:20;

PDF=diff(CDF,Csr);

plot(Csr,PDF)