# HG changeset patch # User A.R. Burgers # Date 1614894133 -3600 # Thu Mar 04 22:42:13 2021 +0100 # Node ID 676db18cb83f00cedf03279ddeabca84821492a1 # Parent 88307452af1a9698dd2a3e5c84541d3043a6cc58 first step towards fixing b60113 * tricdf.m: fix diff -r 88307452af1a -r 676db18cb83f inst/tricdf.m --- a/inst/tricdf.m Mon Mar 01 19:53:49 2021 +0100 +++ b/inst/tricdf.m Thu Mar 04 22:42:13 2021 +0100 @@ -63,7 +63,7 @@ cdf(k_temp) += full_area; k_temp = k & (a < x) & (x < c); - area = (x(k_temp) - a).^2 * h; + area = (x(k_temp) - a).^2 * h / (2 * ( c - a)); cdf(k_temp) += area; k_temp = k & (b <= x); @@ -71,7 +71,7 @@ cdf(k_temp) += full_area; k_temp = k & (c < x) & (x < b); - area = (b-x(k_temp)).^2 * h; + area = (b-x(k_temp)).^2 * h / (2 * (b - c)); cdf(k_temp) += full_area - area; else h = 2 ./ (b-a); @@ -81,7 +81,7 @@ cdf(k_temp) += full_area; k_temp = k & (a <= x) & (x < c); - area = (x(k_temp) - a(k_temp)).^2 .* h(k_temp); + area = (x(k_temp) - a(k_temp)).^2 .* h(k_temp) ./ (2 * (c(k_temp) - a(k_temp))); cdf(k_temp) += area; k_temp = k & (b <= x); @@ -89,7 +89,7 @@ cdf(k_temp) += full_area; k_temp = k & (c <= x) & (x < b); - area = (b(k_temp)-x(k_temp)).^2 .* h(k_temp); + area = (b(k_temp)-x(k_temp)).^2 .* h(k_temp) ./ (2 * (b(k_temp) - c(k_temp))); cdf(k_temp) += full_area - area; endif