File talk:Judd's (r,g) diagram.svg

Source

r=[0	0	0	0	0	0.001	0.005	0.023	0.095	0.136	0.177	0.238	0.287	0.33	0.368	0.399	0.438	0.487	0.55	0.639	0.719	0.801	0.862	0.906	0.951	0.975	0.985	1;
0.001	0.002	0.002	0.002	0.003	0.005	0.014	0.037	0.116	0.157	0.196	0.251	0.294	0.331	0.363	0.389	0.423	0.465	0.52	0.6	0.669	0.743	0.799	0.839	0.882	0.904	0.912	0.928;
0.003	0.004	0.004	0.004	0.006	0.011	0.023	0.053	0.137	0.178	0.214	0.263	0.3	0.331	0.358	0.381	0.409	0.446	0.493	0.561	0.623	0.689	0.739	0.776	0.813	0.834	0.843	0.856;
0.005	0.006	0.006	0.007	0.01	0.017	0.035	0.072	0.159	0.198	0.231	0.274	0.305	0.331	0.354	0.373	0.397	0.427	0.468	0.526	0.58	0.637	0.681	0.713	0.747	0.766	0.774	0.784;
0.008	0.009	0.01	0.011	0.015	0.026	0.049	0.092	0.182	0.218	0.248	0.284	0.31	0.331	0.351	0.366	0.386	0.411	0.444	0.493	0.539	0.588	0.626	0.654	0.683	0.699	0.706	0.716;
0.012	0.013	0.014	0.016	0.022	0.037	0.067	0.117	0.205	0.238	0.264	0.294	0.315	0.332	0.347	0.359	0.375	0.395	0.422	0.462	0.5	0.54	0.572	0.596	0.621	0.635	0.641	0.65;
0.018	0.02	0.021	0.024	0.032	0.052	0.09	0.145	0.229	0.258	0.279	0.303	0.319	0.332	0.344	0.353	0.366	0.381	0.402	0.433	0.463	0.495	0.521	0.54	0.56	0.572	0.576	0.584;
0.028	0.03	0.032	0.036	0.047	0.073	0.121	0.179	0.254	0.277	0.293	0.311	0.323	0.333	0.341	0.348	0.357	0.368	0.383	0.406	0.428	0.452	0.472	0.486	0.501	0.51	0.513	0.52;
0.045	0.048	0.051	0.057	0.073	0.109	0.163	0.22	0.279	0.296	0.307	0.319	0.327	0.333	0.338	0.343	0.348	0.356	0.365	0.381	0.395	0.411	0.424	0.434	0.444	0.45	0.452	0.456;
0.087	0.092	0.096	0.106	0.129	0.174	0.227	0.27	0.306	0.315	0.321	0.326	0.33	0.333	0.336	0.338	0.34	0.344	0.349	0.356	0.363	0.371	0.378	0.382	0.388	0.391	0.392	0.394;
0.142	0.148	0.154	0.165	0.19	0.232	0.273	0.3	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32	0.32];
g=[0	0.002	0.003	0.007	0.017	0.044	0.105	0.203	0.397	0.492	0.562	0.625	0.638	0.624	0.604	0.584	0.551	0.506	0.446	0.359	0.28	0.199	0.138	0.094	0.049	0.025	0.015	0;
0.001	0.004	0.005	0.009	0.02	0.047	0.111	0.209	0.391	0.475	0.534	0.585	0.594	0.582	0.564	0.548	0.52	0.482	0.431	0.355	0.287	0.216	0.161	0.122	0.081	0.059	0.051	0.036;
0.003	0.006	0.007	0.011	0.023	0.052	0.117	0.216	0.386	0.458	0.508	0.548	0.555	0.544	0.53	0.515	0.492	0.459	0.416	0.352	0.293	0.231	0.183	0.148	0.113	0.093	0.084	0.072;
0.005	0.008	0.009	0.014	0.026	0.058	0.125	0.223	0.38	0.442	0.483	0.514	0.519	0.509	0.497	0.485	0.466	0.439	0.403	0.349	0.299	0.246	0.205	0.175	0.143	0.125	0.118	0.108;
0.008	0.011	0.013	0.017	0.031	0.066	0.136	0.232	0.374	0.426	0.458	0.483	0.486	0.478	0.467	0.458	0.442	0.42	0.391	0.347	0.305	0.26	0.225	0.199	0.172	0.158	0.151	0.142;
0.012	0.015	0.017	0.023	0.037	0.075	0.148	0.242	0.368	0.41	0.435	0.453	0.455	0.449	0.44	0.433	0.42	0.403	0.38	0.344	0.31	0.274	0.245	0.224	0.201	0.188	0.183	0.175;
0.018	0.021	0.024	0.031	0.047	0.088	0.164	0.254	0.361	0.394	0.413	0.426	0.427	0.422	0.415	0.41	0.4	0.387	0.369	0.342	0.315	0.287	0.264	0.247	0.229	0.219	0.215	0.208;
0.028	0.032	0.034	0.042	0.061	0.107	0.185	0.268	0.355	0.379	0.392	0.401	0.401	0.397	0.393	0.388	0.381	0.372	0.359	0.339	0.32	0.299	0.282	0.27	0.256	0.249	0.246	0.24;
0.045	0.05	0.053	0.063	0.086	0.138	0.215	0.285	0.348	0.363	0.371	0.377	0.377	0.374	0.372	0.369	0.364	0.358	0.35	0.337	0.325	0.311	0.3	0.292	0.282	0.277	0.275	0.272;
0.087	0.093	0.098	0.111	0.139	0.194	0.26	0.307	0.341	0.348	0.352	0.355	0.354	0.353	0.352	0.35	0.348	0.345	0.341	0.335	0.329	0.322	0.317	0.313	0.308	0.306	0.305	0.303;
0.142	0.149	0.155	0.169	0.197	0.245	0.292	0.319	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337	0.337];
b=[1	0.998	0.997	0.993	0.983	0.955	0.89	0.774	0.508	0.372	0.261	0.137	0.075	0.046	0.028	0.017	0.011	0.007	0.004	0.002	0.001	0	0	0	0	0	0	0;
0.998	0.994	0.993	0.989	0.978	0.948	0.875	0.754	0.493	0.368	0.27	0.164	0.112	0.087	0.073	0.063	0.057	0.053	0.049	0.045	0.044	0.041	0.04	0.039	0.037	0.037	0.037	0.036;
0.994	0.99	0.989	0.985	0.971	0.937	0.86	0.731	0.477	0.364	0.278	0.189	0.145	0.125	0.112	0.104	0.099	0.095	0.091	0.087	0.084	0.08	0.078	0.076	0.074	0.073	0.073	0.072;
0.99	0.986	0.985	0.979	0.964	0.925	0.84	0.705	0.461	0.36	0.286	0.212	0.176	0.16	0.149	0.142	0.137	0.134	0.129	0.125	0.121	0.117	0.114	0.112	0.11	0.109	0.108	0.108;
0.984	0.98	0.977	0.972	0.954	0.908	0.815	0.676	0.444	0.356	0.294	0.233	0.204	0.191	0.182	0.176	0.172	0.169	0.165	0.16	0.156	0.152	0.149	0.147	0.145	0.143	0.143	0.142
0.976	0.972	0.969	0.961	0.941	0.888	0.785	0.641	0.427	0.352	0.301	0.253	0.23	0.219	0.213	0.208	0.205	0.202	0.198	0.194	0.19	0.186	0.183	0.18	0.178	0.177	0.176	0.175;
0.964	0.959	0.955	0.945	0.921	0.86	0.746	0.601	0.41	0.348	0.308	0.271	0.254	0.246	0.241	0.237	0.234	0.232	0.229	0.225	0.222	0.218	0.215	0.213	0.211	0.209	0.209	0.208;
0.944	0.938	0.934	0.922	0.892	0.82	0.694	0.553	0.391	0.344	0.315	0.288	0.276	0.27	0.266	0.264	0.262	0.26	0.258	0.255	0.252	0.249	0.246	0.244	0.243	0.241	0.241	0.24;
0.91	0.902	0.896	0.88	0.841	0.753	0.622	0.495	0.373	0.341	0.322	0.304	0.296	0.293	0.29	0.288	0.288	0.286	0.285	0.282	0.28	0.278	0.276	0.275	0.274	0.273	0.273	0.272;
0.826	0.815	0.806	0.783	0.732	0.632	0.513	0.423	0.353	0.337	0.327	0.319	0.316	0.314	0.312	0.312	0.312	0.311	0.31	0.309	0.308	0.306	0.305	0.305	0.304	0.303	0.303	0.303;
0.716	0.703	0.691	0.666	0.613	0.525	0.435	0.381	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343	0.343];

pr=[0.38	0.37	0.358	0.346	0.339	0.334	0.329	0.325	0.32	0.314	0.312;
0.426	0.405	0.38	0.357	0.344	0.334	0.325	0.318	0.309	0.296	0.294;
0.474	0.439	0.401	0.367	0.349	0.334	0.322	0.312	0.299	0.281	0.277;
0.522	0.472	0.42	0.376	0.353	0.335	0.319	0.306	0.29	0.267	0.262;
0.57	0.505	0.438	0.385	0.357	0.335	0.316	0.301	0.281	0.255	0.249;
0.694	0.582	0.478	0.402	0.364	0.336	0.311	0.291	0.264	0.228	0.221;
0.826	0.654	0.513	0.417	0.37	0.336	0.306	0.282	0.25	0.207	0.198;
0.96	0.723	0.543	0.429	0.375	0.336	0.303	0.275	0.239	0.19	0.179;
1.106	0.788	0.569	0.439	0.379	0.336	0.3	0.27	0.229	0.175	0.163;
0.386	0.386	0.386	0.386	0.386	0.337	0.295	0.261	0.215	0.152	0.138;
0.391	0.391	0.391	0.391	0.391	0.337	0.292	0.254	0.203	0.135	0.12;
0.245	0.245	0.245	0.245	0.245	0.245	0.245	0.245	0.188	0.111	0.094;
0.174	0.174	0.174	0.174	0.174	0.174	0.174	0.174	0.174	0.088	0.069;
0.153	0.153	0.153	0.153	0.153	0.153	0.153	0.153	0.153	0.056	0.033;
0.035	0.035	0.035	0.035	0.035	0.035	0.035	0.035	0.035	0.035	0.012;
0.024	0.024	0.024	0.024	0.024	0.024	0.024	0.024	0.024	0.024	0];
pg=[0.31	0.304	0.298	0.295	0.295	0.297	0.3	0.303	0.306	0.311	0.312;
0.287	0.276	0.265	0.26	0.262	0.266	0.271	0.276	0.282	0.292	0.294;
0.263	0.248	0.235	0.229	0.232	0.238	0.245	0.252	0.261	0.275	0.277;
0.239	0.222	0.207	0.202	0.206	0.213	0.222	0.23	0.242	0.259	0.262;
0.215	0.195	0.18	0.176	0.181	0.191	0.202	0.211	0.225	0.245	0.249;
0.153	0.133	0.121	0.122	0.131	0.144	0.158	0.171	0.189	0.216	0.221;
0.087	0.075	0.07	0.078	0.091	0.107	0.124	0.138	0.161	0.192	0.198;
0.02	0.02	0.027	0.041	0.058	0.077	0.096	0.113	0.137	0.172	0.179;
-0.053	-0.032	-0.012	0.01	0.03	0.051	0.072	0.091	0.117	0.155	0.163;;
-0.014	-0.014	-0.014	-0.014	-0.014	0.011	0.035	0.057	0.087	0.13	0.138;
-0.046	-0.046	-0.046	-0.046	-0.046	-0.018	0.008	0.031	0.063	0.111	0.12;
-0.004	-0.004	-0.004	-0.004	-0.004	-0.004	-0.004	-0.004	0.032	0.084	0.094;
0.001	0.001	0.001	0.001	0.001	0.001	0.001	0.001	0.001	0.057	0.069;
-0.041	-0.041	-0.041	-0.041	-0.041	-0.041	-0.041	-0.041	-0.041	0.021	0.033;
-0.002	-0.002	-0.002	-0.002	-0.002	-0.002	-0.002	-0.002	-0.002	-0.002	0.012;
-0.014	-0.014	-0.014	-0.014	-0.014	-0.014	-0.014	-0.014	-0.014	-0.014	0];
pb=[0.31	0.326	0.344	0.359	0.366	0.369	0.371	0.372	0.374	0.375	0.376;
0.287	0.319	0.355	0.383	0.394	0.4	0.404	0.406	0.409	0.412	0.412;
0.263	0.313	0.364	0.404	0.419	0.428	0.433	0.436	0.44	0.444	0.446;
0.239	0.306	0.373	0.422	0.441	0.452	0.459	0.464	0.468	0.474	0.476;
0.215	0.3	0.382	0.439	0.462	0.474	0.482	0.488	0.494	0.5	0.502;
0.153	0.285	0.401	0.476	0.505	0.52	0.158	0.171	0.547	0.556	0.558;
0.087	0.271	0.417	0.505	0.539	0.557	0.57	0.58	0.589	0.601	0.604;
0.02	0.257	0.43	0.53	0.567	0.587	0.601	0.612	0.624	0.638	0.642;
-0.053	0.244	0.443	0.551	0.591	0.613	0.628	0.639	0.654	0.67	0.674;
0.628	0.628	0.628	0.628	0.628	0.652	0.67	0.682	0.698	0.718	0.724;
0.655	0.655	0.655	0.655	0.655	0.681	0.7	0.715	0.734	0.754	0.76;
0.759	0.759	0.759	0.759	0.759	0.759	0.759	0.759	0.78	0.805	0.812;
0.825	0.825	0.825	0.825	0.825	0.825	0.825	0.825	0.825	0.855	0.862;
0.888	0.888	0.888	0.888	0.888	0.888	0.888	0.888	0.888	0.923	0.934;
0.967	0.967	0.967	0.967	0.967	0.967	0.967	0.967	0.967	0.967	0.976;
0.99	0.99	0.99	0.99	0.99	0.99	0.99	0.99	0.99	0.99	1];

L1=[410:10:490 495 500:10:670];
L2=[492.5 495 500:10:570 571.8];
plot(g',r','k',pg',pr','m'); hold;
plot([repmat(1/3,1,length(L1)); g(1,:)],[repmat(1/3,1,length(L1)); r(1,:)],':','Color',[0.5 0.5 0.5]);
intersection=(pr(1,:)-pg(1,:))./(1-3*pg(1,:));
plot(repmat([1/3; 0],1,length(L2)),[repmat(1/3,1,length(L2)); intersection],':','Color',[0.75 0.75 0.75]);
text(g(1,:),r(1,:),num2str(L1'));
axis('equal',[0 0.7 0 1]);