pde3 := diff(P(t, r, s), t)+k.(a-r).(diff(P(t, r, s), r))+.5*sigma^2*r^(-2*b)*(diff(P(t, r, s), r, r))-2*r*P(t, r, s) = 0